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ABSTRACT 


Two  quasi -static,  steady  state  flow  problems  of  an  incom¬ 
pressible,  isotropic,  elastic-perfectly  plastic  solid  are  considered  in 
this  thesis;  namely,  the  plane  radial  flow  through  a  rough  converging 
channel  and  the  axially  symmetric  radial  flow  through  a  rough  conver¬ 
ging  conical  channel.  The  classical  PRANDTL-REUSS  constitutive  equations 
are  rewritten  in  a  form  which  is  invariant  under  arbitrary  changes  of 
the  frame  of  reference.  This  invariant  form  is  necessary  for  the 
problems  considered  since  the  classical  form  is  applicable  only  to 
problems  which  involve  small  deformations  and  for  which  the  use  of  the 
infinitesimal  strain  tensor  is  justified. 

The  relationship  between  the  form  invariant  PRANDTL-REUSS 
equations  and  the  constitutive  equations  of  a  hypo-elastic  solid  of 
grade  two  is  discussed.  Also,  two  steady  state  flow  problems  of  a  hypo- 
elastic  solid  of  grade  two  are  considered:  the  plane  radial  flow  through 
a  smooth  converging  channel  with  and  without  consideration  of  inertia 
effects . 

A  frame  indifferent  definition  of  stress  rate  is  of  critical 
importance  in  modern  theories  of  hypo-elasticity  and  elasto-plasticity. 
Consequently,  a  discussion  on  the  definitions  of  stress  rate  is  contained 
in  this  thesis. 

An  investigation  is  made  of  the  nature  of  the  governing  equations 
for  the  plane  strain  flow  of  an  elastic-perfectly  plastic  solid.  The 
results  of  this  investigation  are  applied  in  part  to  the  plane  radial 
flow  problem  of  the  elastic-perfectly  plastic  solid. 
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CHAPTER  I 


CONSTITUTIVE  EQUATIONS  IN  THE  FLOW  THEORY  OF 
ELASTIC-PERFECTLY  PLASTIC  MEDIA  FOR  INFINITESIMAL  STRAIN 

This  chapter  contains  a  review  of  the  classical  elastic- 
perfectly  plastic  theory  of  PRANDTL-REUSS.  This  theory  is  applicable 
only  to  problems  involving  small  deformations  since  approximations  are 
made  for  the  actual  strains  and  stress  rates.  In  this  chapter  all  dis¬ 
placement  components  u^i  =  1,2,3)  are  referred  to  a  fixed  rectangular 
CARTESIAN  frame  X ^. ( i  =  1,2,3)  and  are  assumed  to  be  single-valued  con¬ 
tinuous  functions  of  the  spatial  coordinates  x.  or  the  material  coordi¬ 
nates  a.  and  time  t.  The  review  begins  with  definitions  characterizing 
small  deformations  and  the  usual  CARTESIAN  tensor  notation  is  used 
throughout. 

1.1  INFINITESIMAL  STRAIN  TENSORS  AND  STRAIN  RATE  TENSORS 

The  classical  PRANDTL-REUSS  theory  is  valid  for  small  deforma¬ 
tions,  that  is,  for  deformations  that  involve  small  displacements,  small 
strains  and  small  rotations.  The  precise  meaning  of  the  term  small  is 
discussed  in  TRUESDELL  and  TOUPIN  [1]. 

GREEN'S  strain  tensor  is  given  by 
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can  be  neglected  in  comparison  with  E..  .,  then 


ALMANSI'S  strain  tensor  is  given  by 
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and  for  small  deformations 
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The  displacements  u.  =  x.  -  a.  if  small  imply  that 
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where 
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J  * 

defines  the  infinitesimal  strain  tensor  and  no  distinction  is  made  between 
spatial  and  material  variables. 

The  material  rate  of  change  of  GREEN'S  strain  tensor  is  given 

by 


DE^l  9^^ 

Dt  st  9a .  3a. 

J 


where 


1  3vc  3vt 
dst  =  2  + 


(1.1.1) 


are  components  of  the  rate  of  deformation  tensor.  The  rate  of  defor¬ 
mation  tensor  is  of  importance  in  constitutive  equations  since  it  is  a 
tensor  that  is  independent  of  the  rigid  body  rotation  of  a  material. 

In  plasticity  literature  d  t  are  usually  called  the  components  of  the 
strain  rate  tensor  and  henceforth  in  this  thesis  this  terminology  will 
be  used. 

The  infinitesimal  strain  increment  during  the  infinitesimal  time 
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increment  dt  may  be  defined  as 


dt 


and 
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where  du^  is  the  displacement  increment  and  x.  are  spatial  variables. 
Another  definition  of  finite  strain  is 


/  d.  dt 


where  the  integration  is  taken  over  the  strain  path. 


1.2  ELASTIC  STRESS-STRAIN  RELATIONS  FOR  SMALL  STRAINS 

The  stress-strain  relation  for  the  classical  HOOKEAN  solid  is 
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where  a.  .  are  components  of  the  stress  tensor  and  A  and  y  are  the  LAMf! 

•  J 

constants.  This  stress-strain  relation  can  be  expressed  in  terms  of  the 
stress  deviation  components  s..  and  the  strain  deviation  components 

'  J 


(1.2.1a) 
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and 


P  =  -  K  ekk  ,  (1.2.1b) 

1  2 

where  p  =  -  y  a ^  and  K  =  A  +  j  y  is  the  bulk  modulus.  For  an  in¬ 
compressible  HOOKEAN  solid, 

Ekk  =  0  and  £ij  =  eij  • 

Equations  (1.2.1a,b),  differentiated  with  respect  to  time, 

gives 


...  =  2u  f . . ,  f..  =  d.. 
U  U  1J  ij 


‘  3  dkk 


1J 


and 


(1.2.2a,b) 


P  =  -  K  dkk  • 


The  time  derivative  of  stress  is  s. .  =  since  the  deformation  is 

3sii 

small.  For  a  HOOKEAN  solid,  the  convected  terms  v.  ^  ^  may  be  neglected, 

k  axk 

but  for  small  deformation  of  an  elastic-plastic  solid  these  terms  may 
not  be  neglected  if  the  rate  of  work-hardening  is  small  compared  with  E, 
YOUNG'S  modulus  of  the  material. 


1.3  PERFECTLY  PLASTIC  STRESS-STRAIN  RATE  RELATIONS  FOR  SMALL  DEFORMATIONS 

Plasticity  theory  of  the  incremental  or  flow  type  relates  the 
plastic  strain  increment  or  strain  rate  components  and  not  the  plastic 
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strain  components  with  the  current  stress  components.  ST.  VENANT  [2] 
proposed  that  during  two  dimensional  plastic  deformation  of  a  solid  a 
co-axial  relationship  exists  between  the  strain-rate  tensor  and  the 


stress  tensor.  Independently,  LEVY  [3]  and  VON  MISES  [4]  proposed  a 
similar  relationship  for  three  dimensional  plastic  deformation. 


In  the  following  derivations  of  the  constitutive  equations  for 


the  perfectly  plastic  theory  of  small  deformations  the  approach  used 
involves  the  concept  of  a  yield  function  and  its  associated  flow  rule 
based  on  the  existence  of  a  potential  function  taken  to  be  the  yield 
function.  This  method  is  due  to  VON  MISES  [5]  and  MELAN  [6]. 

According  to  the  concept  of  a  yield  function  there  exists  a 
function  f  of  the  stress  such  that 


f(aij)  =  C,  C  >  0, 


(1.3.1) 


and  where  C  is  a  material  parameter  dependent  on  the  strain  history. 
Relationship  (1.3.1)  is  called  a  yield  criterion.  The  condition 


f(a.  .)  <  C  identifies  the  elastic  domain  and  f(a. .)  >_  C  the  plastic 
I J  i  J 


domain.  For  a  non-hardening  material,  C  is  a  material  constant  and 


discussion  is  confined  strictly  to  non-hardening  materials. 

The  conditions  of  isotropy  and  the  independence  of  plastic 
deformation  on  hydrostatic  pressure  require  f  to  be  of  the  form 


f(0£,  J(j)  =  c 


(1.3.2) 


where 
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•  • 
ij 


and 


J 


*jk  ki 


are  respectively  the  second  and  third  stress  deviation  invariants. 
Absence  of  a  BAUSCHINGER  effect  [fta  —  )  =  f ( j ) ]  requires  f  to  be  an 
even  function  in  d3  • 

Two  important  yield  criteria  are  the  VON  MISES  yield  criterion 
and  the  TRESCA  yield  criterion.  The  VON  MISES  yield  criterion  is 


(1.3.3) 


which  written  in  terms  of  the  principal  stress  components  (i  =  1,2,3) 
is 


(a1  -  a2)2  +  (a2  -  a3)2  +  (a3  -  a^2  =  6k2 

where  k  is  the  yield  stress  in  simple  shear.  The  TRESCA  yield  criterion, 
although  it  can  be  expressed  in  the  form  of  equation  (1.3.2),  is  written 

4 

more  conveniently  as 


(1.3.4) 


where  Y  is  the  uniaxial  yield  stress  and  cmx  and  amin  are  the  maximum 
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and  minimum  principal  stresses. 

Depicted  in  a  principal  stress  space  with  the  principal  stress 
a-j ,  G2»  and  a ^  taken  as  rectangular  CARTESIAN  coordinates,  each  criterion 
represents  a  surface  called  a  yield  surface.  The  VON  MISES  yield  surface 
is  a  circular  cylinder  and  the  TRESCA  yield  surface  is  a  regular  hexagonal 
cylinder  both  infinite  in  extent  and  having  their  axes  pass  through  the 
origin  with  direction  numbers  (1,1,1).  Since  plastic  deformation  is  in¬ 
dependent  of  the  hydrostatic  pressure,  any  section  parallel  to  the  plane 
c?i  +  +  cjg  =  0  may  be  taken  as  representative  of  the  stress  domain. 

It  is  customary  to  take  the  orthogonal  projection  of  a  section  parallel 
to  =  constant  onto  the  o- 1  -  o ^  plane.  Figure  1  shows  such  a  pro¬ 
jection,  the  boundary  of  which  is  called  a  yield  locus,  and  illustrates 
the  relationship  taken  between  the  VON  MISES  yield  locus  and  the  TRESCA 
yield  locus.  The  TRESCA  yield  locus  is  a  linearization  of  the  VON  MISES 
yield  locus. 

Convexity  of  yield  surfaces  is  fundamental  in  classical  plas¬ 
ticity  theory.  A  yield  surface  f(o..)  =  C  is  convex  if 

'  vJ 


df 


f(aij>  -  f(aij)  -  (aij  '  °ij)  W~. 


(1.3.5) 


for  any  two  arbitrary  stress  states  at.  and  a.,  producing  plastic  de- 

'  J  *  vJ 

formation.  Inequality  in  relationship  (1.3.5)  holds  only  for  strictly 
convex  surfaces.  The  convexity  of  the  VON  MISES  and  TRESCA  yield  surfaces 
are  evident  from  a  geometrical  standpoint.  However,  convexity  of  any 
general  yield  surface  cannot  be  shown  using  assumptions  of  perfect  plas¬ 
ticity  alone  unless  some  additional  hypothesis  such  as  the  principle  of 
maximum  work  [7]  is  made. 
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IGURE  1  -  THE  VON  MISES  AND  TRESCA  YIELD  LOCI 
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The  hypothesis  of  the  existence  of  a  plastic  potential  states 
that  there  exists  a  function  g  of  the  stress  such  that  the  plastic  strain 
rate  components  are  given  by  the  flow-rule 


A  , 


where  X  is  a  non-negative  scalar  invariant  dependent  only  on  coordinates 
x.  and  time  t.  The  dependence  of  g  is  upon  the  stress  deviation  in¬ 
variants  and  in  order  that  the  co-axial  relationship  between  the 
plastic  strain-rate  tensor  and  the  stress  deviation  tensor  be  satisfied. 

A  further  hypothesis  is  made  that  g  is  identical  to  f  [8].  Thus  in 
principal  stress  space,  the  principal  plastic  strain  rate  vector  Gct(d!^) 
associated  with  the  principal  stress  vector  0(0^.)  has  the  same  direction 
as  the  outward  normal,  if  defined,  to  the  yield  surface  f(cr  . .)  =  C  at 

vJ 

the  point  a.(i  =  1,2,3)  representing  the  plastic  stress  state.  The 
dj^  (i  =  1,2,3)  are  the  principal  components  of  the  plastic  strain  rate 
tensor  and  the  factor  G  is  used  to  give  dimensions  of  stress  to  the 
plastic  strain  rate  vector  ?(d!^)  in  principal  stress  space.  The  plastic 
strain  rate  components  for  perfectly  plastic  materials  are  thus  deter¬ 
mined  by 


d!p>  . 

1J 


9f 


9Qij 


X 


(1.3.6) 


Where  X  =  0  if  f  <  C  and  also  if  f  =  C  and  f  <  0  , 

X  >  0  if  f  =  C  and  f  =  0. 

An  upper  dot  associated  with  a  quantity  denotes  differentiation  with 
respect  to  time  or  any  other  monotonical ly  increasing  parameter  cor- 


YftosVoq  io5*  2Jn9noqmoo  nr  &Tf2 
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related  with  progressive  deformation.  The  flow  rule  (1.3.6)  corresponding 
to  the  particular  yield  function  f  is  called  the  associated  flow  rule  for 
f  .  The  most  general  associated  flow  rule  for  f  in  the  form  given  by 
equation  (1.3.2)  is 


alttt 


3f 


(s, 


3J«  i j  3J '  v  ik  jk 


1  ,, 

"  3  J2 


V] 


Components  of  the  plastic  strain  rate  tensor  are  determined 

up  to  an  arbitrary  non-negative  scalar  factor  at  a  plastic  stress  state 
3  f 

provided  r -  are  uniquely  defined  there.  This  is  possible  for  all 

i  j 

points  on  the  yield  locus  if  the  corresponding  yield  surface  is  a  regular 
surface.  The  VON  MISES  yield  surface  is  regular;  hence  the  plastic 
strain  rate  components  are  derivable  using  the  flow  rule  (1.3.6)  to  give 


d-j1  ■  X  sij  ,  (1.3.7) 

the  LEVY-MISES  flow  relation.  If,  however,  the  yield  surface  is  a  singu¬ 
lar  surface,  the  corresponding  associated  flow  rule  is  given  by  the 
KOITER  generalization  of  the  plastic  potential  [9].  According  to  this 
generalization,  the  plastic  strain  rate  components  are  determined  by 


=  i 


3f. 


3a.  . 

r  ij 


(1.3.8) 


Xr  are  non-negative  scalar  invariants  dependent  only  on  coordinates  x. 
and  time  t  and 


Xr  =  0  if  fr  <  C  and  also  if  fr  =  C,  ff  <  0 
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1  >  0  if  f  =  C  and  f  =  0. 
r  —  r  r 


The  summation  is  taken  over  all  surfaces  f  =  C  whose  join  forms  the 
yield  surface  in  the  neighbourhood  of  the  plastic  stress  point  considered. 

x.  /  D  \ 

In  principal  stress  space,  the  direction  of  Gd ( d^  ')  at  a  singular  point 
must  be  within  the  region  spanned  by  the  unique  normals  drawn  outwards 
to  the  faces  intersecting  at  the  particular  singular  point  considered. 

There  is  only  one  plastic  regime  available  as  the  stress  point 
traverses  the  VON  MISES  yield  locus.  However  the  TRESCA  yield  surface  is 
singular  and  as  the  stress  point  traverses  the  TRESCA  yield  locus  different 
plastic  regimes  are  available.  If  the  stress  point  lies  on  a  flat  of  the 
yield  locus  then, from  the  concept  of  the  plastic  potential, 


d<P>  =  X 
max 


=  -  d 


(P) 

’min 


dint  ‘  0 


where  d^  ,  d(^2  and  d^  are  the  maximum,  intermediate  and  minimum 
principal  components  of  the  plastic  strain  rate  and  X  >_  0  .  If  the 
stress  point  lies  on  a  corner,  for  example  point  B  of  Figure  1,  it 
follows  from  KOITER'S  generalization  of  the  plastic  potential  that 

djP)  =  A1  ,  d£P)  =  X2>  d<P)  =  -  (A,  +  A2) 


where  A,  >0,  A2  >  0.  For  point  C 

djP)  =  -  Ar  d£P)  =  A1  +  A2,  d*R)  =  -  A2 

where  X-j  >_  0, 


J  o  too  9  n  .t  ,Tm  t  .  i  a u  o  f  b  F  9  r  x, 
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>• 

The  form  of  the  scalar  invariant  A  in  the  LEVY-MISES  flow 
relation  (1.3.7)  is  determined  by  squaring  this  equation  to  give 


2 

2A*J 


i 


2  ‘ 


Since  A  d!^ 

I  J 


d(p)  =  v 
ai j  “  l2 


the  second  invariant  of  d 


(P) 

'ij 


it  follows  that 


A  is  thus  a  function  of  the  invariants  of  the  plastic  strain  rate  tensor 
and  stress  deviation  tensor.  Alternatively,  multiplying  both  sides  of 
equation  (1.3.7)  by  s..  gives 

*  J 


2AJ 


2 


from  which 


d(P> 
=  _Ll_ 


'Am  Am 


s..  d(p) 
_1L  a  _U_ 
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2k‘ 


The  LEVY-MISES  flow  relation  for  perfectly  plastic  theory  is  now  rewritten 
as 


d(P)  s 

j(P)  Am  Am  „ 

d .  .  =  - n -  s  .  . 

U  2k2  U 


(1.3.9a) 


or 


(1.3.9b) 


Each  relation  (1.3.9a  or  b)  is  homogeneous  of  order  one  in  the  plastic 
strain  rates.  Thus  perfectly  plastic  flow  is  inviscid. 

1 .4  ST.  VENANT-LEVY-MISES  ELASTIC-PERFECTLY  PLASTIC  THEORY 

The  classical  theory  describing  the  response  of  an  elastic- 

perfectly  plastic  material  in  which  no  account  is  made  for  elastic 

effects  during  plastic  deformation  is  called  the  ST.  VENANT-LEVY-MISES 

theory.  If  the  yielding  of  the  material  is  governed  by  the  VON  MISES 

yield  criterion,  the  theory  is  referred  to  as  the  LEVY-MISES  theory  [10] 

>  >• 

For  analysis  employing  the  ST.  VENANT-LEVY-MISES  theory,  two 
sets  of  constitutive  equations  are  required;  one  for  the  region  de¬ 
forming  elastically  and  the  other  for  the  region  deforming  plastically. 
For  the  elastic  region,  the  constitutive  equations  are  given  by 


£ij  >  Ekk  =  0  • 


for  incompressible  deformations  and  by 


sij  =  2w  £ij  •  p  =  '  Kekk  ’ 

for  compressible  deformations.  For  the  plastic  region,  the  constitutive 
equations  depend  upon  the  yield  criterion  used.  With  the  VON  MISES 
yield  criterion  (1.3.3)  written  as 
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f  =  s. .  s..  =  2k  , 


0.4.1) 


the  constitutive  equations  are  the  LEVY-MISES  relation  (1.3.7),  namely, 


dj'l-n.. 

TJ  TJ 


d.4.2) 


where 


X  >  0  for  f  =  2k  and  f  =  0. 


With  the  TRESCA  yield  criterion  (1.3.4),  the  constitutive  equations  are 


d'p)  -  i  dfr 

U  k 


X 


3  r  3aij  T 


(1.4.3) 


where 


• 

X„  >  0  for  f  =  2k  and  f  =  0. 
r  —  r  r 

In  any  specific  el astic-perfectly  plastic  deformation  problem 
where  the  equations  (1.4.3)  are  used,  the  plastic  regime  applicable  must 
first  be  determined.  However,  this  is  not  always  possible  and  an  additional 
assumption  is  made  to  make  the  problem  either  statically  or  kinematically 
determinate.  This  amounts  to  specifying  the  singular  point  or  flat  used 

*  a 

on  the  TRESCA  yield  locus.  One  such  assumption  is  the  HAAR-KARMAN  hypo¬ 
thesis  [11]  which  states  that  in  some  axially  symmetric  perfectly  plastic 
problems,  the  circumferential  stress  is  equal  to  one  of  the  principal 
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stresses  in  the  meridional  plane*  This  situation  does  not  occur  for 
the  LEVY-MISES  theory  since  there  is  only  one  plastic  regime  available. 

A 

However,  the  LEVY-MISES  theory  has  one  undesirable  property  in  that  the 
governing  equations  for  axially  symmetric  flow  problems  are  non-hyperbol ic 
in  character  thus  adding  difficulties  to  the  mathematical  analysis  [12]. 

A  feature  of  the  ST.  VENANT-LEVY-MISES  theory  which  causes 
difficulty  is  the  determination  of  the  elastic-plastic  boundary. 

If  in  the  deformation  of  an  elastic-perfectly  plastic  material, 
elastic  deformations  are  taken  as  vanishing  identically  so  that  the  de¬ 
formation  of  the  material  is  just  that  due  to  the  perfectly  plastic  de¬ 
formation,  one  is  then  led  to  the  'rigid-plastic'  theory.  This  theory 
is  applicable  strictly  to  a  hypothetical  'rigid-plastic'  material  which 
remains  rigid  for  stress  states  below  yield  and  whose  elastic  moduli  are 
indefinitely  large.  The  constitutive  equations  are  given  by  equation 
(1.4.2)  or  equation  (1.4.3)  with  d.  ,'  replaced  by  d.  ..  This  theory  has 

*  J  *  J 

been  applied  to  the  problems  involving  unconstrained  plastic  flow  of  an 
elastic-perfectly  plastic  solid  where  the  elastic  deformations  are  small 
compared  to  the  plastic  deformations  and  so  may  be  justifiably  neglected 
provided  the  analysis  is  not  in  the  immediate  vicinity  of  the  transition 
zone  between  elastic  and  plastic  zones. 

1.5  CLASSICAL  ELASTIC-PERFECTLY  PLASTIC  THEORY  OF  PRANDTL-REUSS 

The  incorporation  of  the  elastic  strains  in  the  stress-rate 
of  strain  relations  for  a  perfectly  plastic  solid  was  proposed  by 
PRANDTL  [13]  for  two  dimensional  flow  problems  and  by  REUSS  [14]  for 
three  dimensional  flow  problems.  The  definitions  of  strain  and  stress 
rate  used  are  applicable  only  for  small  deformations. 


. 
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During  the  continued  loading  of  an  element  of  an  elastic- 
perfectly  plastic  material  deforming  plastically,  the  total  displace¬ 
ment  of  each  material  point  of  the  element  may  be  represented  by 


+  u 


(P) 


( 6 ) 

where  u)  '  is  the  displacement  vector  corresponding  to  the  elastic 

( p) 

deformation,  and  u)  '  is  the  displacement  vector  corresponding  to  the 

plastic  deformation  remaining  after  unloading.  The  decomposition  of 

( 6 ) 

u.  into  elastic  and  plastic  components  is  not  unique  since  u)  '  and 
ujP>  are  determined  up  to  a  rigid  body  displacement  of  the  whole  material. 
The  same  remarks  apply  to  the  velocity  vector  given  by 


v 


i 


+  v 


(P) 


(1.5.1) 


where  v?e^  and  vjP^  are  the  velocity  vectors  corresponding  to  the  elastic 
and  plastic  deformations  respectively.  Consequently* the  strain  rate 
tensor  defined  by  equation  (1.1.1)  is  used  in  formulating  the  classical 
PRANDTL-REUSS  equations  since  it  is  a  tensor  independent  of  the  rigid 
body  rotation  of  the  material.  From  equation  (1.5.1)  it  follows  that 


d.  .  =  d!®'  +  d(^ 

ij  ij  ij 


(1.5.2) 


where  d!®^  and  d..  are  respectively  the  uniquely  determined  components 
1  J  •  J 

of  the  elastic  and  plastic  strain  rate  tensors  defined  by 

,  v  ,  8v(e^  8v(e) 

d(e)  =  1  ( — ] —  +  — J —  ) 
ij  “  2  v3x,  3x.  ; 

J 
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and 


d(p)  =  1  (!!i 

ij  “  2  9x . 


(P) 


In  further  discussion.a  solid  which  yields  plastically  according  to 
the  VON  MISES  yield  criterion  will  be  called  a  MISES  solid  and  one 
which  yields  plastically  according  to  the  TRESCA  yield  criterion  will 
be  called  a  TRESCA  solid. 

For  an  elastic-perfectly  plastic  MISES  solid  which  is  elasti¬ 
cally  compressible  and  which  is  deforming  plastically,  the  strain  rate 
tensor  is  obtained  from  equation  (1.5.2)  rewritten  as 


:(e) 

ij 


1 

3 


iff  ^  6  .  . 

kk  ij 


+  d 


(P) 

ij 


(1.5.3) 


(e) 

where  f. .  are  components  of  the  elastic  strain  rate  deviation  tensor. 
■  J 

From  equations  (1.2.2a,b)  it  then  follows  that 


1_  m 


Ds. 


2p  Dt 


1_  D£ 

"  3K  Dt 


+  X  s 


ij 


(1.5.4a) 


This  together  with 


d 


1 

K  Dt 


(1.5.4b) 


constitute  the  constitutive  equations  for  the  MISES  solid.  The  operator 

_  n  s  ^ 

jrr  is  the  material  derivative  defined  by  ^  (  )  =  ^  (  )  +  vk  (  ). 

ut  Ak 

However,  differentiation  need  not  be  with  respect  to  time  t  but  to  any 
other  monotonically  increasing  parameter  correlated  with  progressive  de- 
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formation.  This  is  possible  since  plastic  deformation  is  inviscid. 

(e) 

If  d\  '  =  0,  it  follows  from  equations  ( 1 . 5 . 4a , b )  that 


1 


Ds .  . 
11 


2p  Dt  +  A  sij 


and 


(1 .5.5a ,b) 


Equations  (1.5.5a,b)  are  the  constitutive  equations  for  an  incompressible 
elastic-perfectly  plastic  MISES  solid. 

An  expression  for  X  in  the  case  X  f  0  can  be  obtained  by 
multiplying  each  side  of  equation  (1.5.4a)  or  equation  (1.5.5a)  by  s. . 

*  J 

and  using  the  result  on  differentiating  materially  the  VON  MISES  yield 
criterion  (1.4.1)  that 


Df  .  2s  !!ii  =  0 

Dt  “1  j  Dt  u* 


This,  along  with  the  VON  MISES  yield  criterion,  yields 


d .  .  s  . . 

A  0  9 

2k 


(1.5.6) 


For  an  elastic-perfectly  plastic  TRESCA  solid  which  is  elasti¬ 
cally  compressible  and  deforming  plastically,  the  constitutive  equations 


are 
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d 


•  • 


ij 


1  Dsi,i 

2  y  Dt 


!_  2fi. 

3K  Dt 


and 


(1.5.7a,b) 


d  =  .  1 
kk  K  Dt  ' 


The  role  of  Ar  is  as  given  in  relationship  (1.3.8).  Equations  (1.5.7a,b) 

follow  analogously  from  equations  (1.5.3)  as  did  equations  (1.5. 4a ,b )  on 

(P) 

simply  replacing  d..  by  its  form  given  in  equation  (1.3.8). 

'  J 

If  d^  =  d^  =  0,  then  from  equations  (1.5.7a,b)  one  obtains 


l  Dsii 

aij  “  2y  Dt  +  l  V  3 s i j 


and  (1 . 5 . 8a ,b) 


These  are  the  constitutive  equations  for  an  incompressible  elastic- 
perfectly  plastic  TRESCA  solid.  The  role  of  is  as  given  in  relation 
(1.3.8). 

From  the  constitutive  equations  of  the  classical  PRANDTL-REUSS 
equations  considered  in  the  foregoing  discussion,  those  for  the  rigid 
plastic  theory  are  derivable  on  letting  y  be  indefinitely  large,  pro¬ 
vided  that  all  stress  derivatives  and  velocities  are  finite. 

3s .  . 

On  first  neglecting  convected  terms  of  type  vk  ,  inte- 
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gration  with  respect  to  time  t  of  equations  (1.5. 4a 9b)  or  (1.5.7a,b) 
with  X,  X^  =  0  yields 

£.  .J-s 

ij  2u  ij 

and 


e 


These  are  the  constitutive  equations  of  a  compressible  HOOKEAN  solid 
for  infinitesmal  strains  under  the  usual  initial  conditions  of  a  stress 
free-strain  free  natural  state.  Similarly  from  equations  (1.5.5a,b) 
or  (1.5.8a,b),  one  obtains 


2p  Sij 


and 


the  constitutive  equations  of  an  incompressible  HOOKEAN  solid. 

The  constitutive  equations  of  the  classical  PRANDTL-REUSS  theory 

>  * 

are  similar  to  those  of  ST.  VENANT-LEVY-MISES  theory  with  the  main  exception 
that  elastic  deformations  are  accounted  for  in  determining  the  deformations 
occurring  during  plastic  flow.  There  are  still  two  sets  of  constitutive 
equations  needed;  one  valid  for  the  elastic  region  and  the  other  in  the 
plastic  region.  For  loading  systems  passing  from  elastic  to  plastic 
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states,  a  determination  of  the  elastic-plastic  boundary  is  still  re¬ 
quired.  However,  the  constitutive  equations  have  been  made  more  compli¬ 
cated  by  the  introduction  of  stress  rates. 

Attempts  at  formulating  an  elastic-perfectly  plastic  theory 
which  would  lead  to  a  gradual  transition  from  elastic  to  plastic  states 
with  the  use  of  just  a  single  set  of  constitutive  equations  valid  for 
any  stress  state  not  violating  the  yield  criterion,  were  made  by 
PRAGER  [15]  and  THOMAS  [16].  THOMAS'S  work  is  a  generalization  and 
extension  of  an  earlier  work  [17]  on  axiomatizing  a  perfectly  plastic 
theory  valid  for  infinitesimal  strain.  This  latter  work  is  presently 
considered. 

1 .6  PERFECTLY  PLASTIC  THEORY  OF  T.Y.  THOMAS 

T.Y.  THOMAS  [17]  proposed  a  theory  aimed  at  revealing  the 
possible  forms  that  the  constitutive  equations  of  a  perfectly  plastic 
material  may  take  on  using  established  characteristics  of  this  material 
under  plastic  deformation.  To  formulate  this  new  mathematical  science, 
THOMAS  proposed: 

Axiom  1.  Perfectly  plastic  deformation  occurs  without  volume  change, 
that  is, 


d 


Axiom  2.  The  stress  deviation  tensor  is  an  isotropic  tensor  function 
of  the  plastic  strain  rate  deviation  tensor,  that  is. 


(p) 


ij  v  m  ' 


9 


(1.6.1) 
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Axiom  3.  Relationship  (1.6.1)  does  not  establish  a  one-to-one  corres¬ 
pondence  between  the  stress  deviation  tensor  and  the  plastic  strain 
rate  deviation  tensor. 

With  these  axioms,  THOMAS  obtained  the  constitutive  equations 
for  the  LEVY-MISES  theory  and  constitutive  equations  for  a  perfectly 
plastic  TRESCA  solid.  This  latter  set  of  constitutive  equations  differ 
from  those  given  previously  in  SECTION  1.3.  The  perfectly  plastic 
theory  proposed  by  THOMAS  is  not  consistent  with  the  perfectly  plastic 
theory  based  on  the  concept  of  a  plastic  potential  identified  as  the 
yield  function. 

Axioms  1  and  2  imply  that 


sij  =  <I>ij(I2,I3^  (1.6.2) 

where  I2  =  d!^  d(^  and  I3  =  d!^  dj^  d^  are  the  two  non-zero  basic 
invariants  of  the  plastic  strain  rate  tensor.  Since  .  are  components 

*  J 

of  a  symmetric  second  order  tensor  and  dependent  on  invariants  and 

I  and  since  s..  =0,  equation  (1.6.2)  can  be  expressed  [18]  as 
J 


ij 


=  a 


d(p> 

U 


+  a0  (d 


(P)  ^(P)  1 


2  v  i  k  ukj 


"  3  *2 


6.  .) 

ij' 


(1.6.3) 


where  a.  =  a.  (I9,lJ,  i  =  1,2.  Multiplication  of  equation  (1.6.3)  by 

1  1  C.  ■J 

J  s1j  and  1  s j k  ski  alternately  yields 

a  2 

J2  =  — I2  +  aia2^3  +  T2  a2^^2^  (1.6. 4. a) 


and 
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a. 


j;  = 


a 


2  T  2 
ol2I2 


+  1 


a*j  ot  ^ 


yh 


a. 


(I 


1  T  3\ 
T2  l2  h 


(1.6.4b) 


By  axiom  3,  the  JACOBIAN  A  of  the  equations  (1.6.4a,b)  is  zero.  This 
results  in  a  partial  differential  equation  involving  the  unknown 
and  c^.  Guided  by  the  existing  forms  of  the  previously  established 
constitutive  equations  of  the  LEVY-MISES  theory,  THOMAS  considered  the 
specific  case  where  =  0  and  a-j  >  0.  The  quasi-linear  partial  dif¬ 
ferential  equation 


9a-.  9ou 

2I2  31^  +  3I3  +  al  =  0 


so  formed  has  the  general  solution 


a 


1 


) 


for  arbitrary  differentiable  function  F.  Thus  a  general  form  for  the 
constitutive  equations  of  perfect  plasticity  is 


(1.6.5) 


Squaring  equation  (1.6.5),  one  obtains  the  relation 


(1.6.6) 
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which  plays  the  role  of  a  yield  criterion  since  it  is  a  restriction  on 
the  stress  deviation  invariants.  For  F  a  constant  function  defined  by 


F ( J J 3 )  =  k  a  material  constant, 


equation  (1.6.6)  becomes 


sij  sij  =  2k  * 

the  VON  MISES  yield  criterion  with  the  appropriate  interpretation  of  k. 
Also  equation  (1.6.5)  becomes 


Z(P)  d  P 
,(P)  /d£m  % 
aij  s-’-* 


/2k 


ij 


d(P>  s 

£m  s£m 


2k 


7 


s  •  •  , 

ij 


the  constitutive  equations  for  the  LEVY-MISES  theory. 

Alternative  representations  of  equation  (1.6.5)  and  equation 

(1.6.6)  are  available  if  equation  (1.6.5)  is  expressed  in  terms  of  si 

and  d(P^  ,  the  principal  components  of  the  stress  deviation  tensor  and 

(P) 

plastic  strain  rate  tensor  respectively.  If  6)  '  are  ordered  such 
that  djP^  >  d^  >  d^  ,  then  equation  (1.6.5)  can  be  written  as 


G(d'P>/d<P>) 


si 


TFT 


-  d 


1  'h(P) 

(FT  d1 


(1.6.7) 


for  arbitrary  differentiable  function  G  >  0.  From  this  one  gets 
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s 


(1.6.8) 


a  relation  analogous  to  relation  (1.6.6). 
material  constant, equation  (1.6.8)  becomes 


Choosing  G(-^-)  =  Y  ,  Y  a 
S1 


the  TRESCA  yield  criterion  for  plastic  regime  AB  in  Figure  1  for  Y 
interpreted  as  the  yield  stress  in  simple  shear.  The  corresponding 
constitutive  equations  are 


s 


(1.6.9) 


as  follows  from  equation  (1.6.7).  These,  however,  differ  from  the 
constitutive  equations  for  plastic  regime  AB  using  the  concept  of  the 
plastic  potential  since  then  d£  '  =  0,  whereas  from  equation  (1.6.9), 
d^  t  0  necessarily. 


Y 


CHAPTER  II 


CONSTITUTIVE  EQUATIONS  OF  THE  ELASTIC-PERFECTLY 

PLASTIC  THEORY  OF  PRANDTL-REUSS  FOR  FINITE  STRAIN 

The  constitutive  equations  for  the  finite  elastic-perfectly 
plastic  theory  of  PRANDTL-REUSS  are  developed  in  this  chapter.  This 
necessarily  involves  a  preliminary  discussion  on  invariant  forms  of 
constitutive  equations;  a  topic  of  recent  interest  in  non-linear  field 
theory  of  mechanics.  The  development  of  this  topic  has  been  varied 
both  in  approach  and  in  notation. 

As  seen  in  Chapter  I,  the  constitutive  equations  of  the  classi¬ 
cal  PRANDTL-REUSS  theory  are  suitable  only  for  problems  involving 
infinitesimal  strains  and  rotations.  A  generalization  of  these  equations 
valid  with  respect  to  a  fixed  reference  system  and  applicable  for  finite 
strains  has  been  achieved  by  THOMAS  [18].  Primarily  it  is  the  presence 
of  the  material  derivative  of  the  stress  tensor  that  invalidates  the 
use  of  the  classical  PRANDTL-REUSS  equations  in  problems  involving 
finite  deformation.  The  reason  for  this  is  that  these  constitutive 
equations  violate  the  'principle  of  material  objectivity',  the  invariance 
principle  of  constitutive  equations  under  arbitrary  changes  of  the  frame 
of  reference  [19], 

The  local  motion  of  a  deforming  material  element  surrounding 
a  point  P  can  be  resolved  into  a  translation,  an  instantaneous  rigid 
body  rotation  about  an  axis  through  P,  a  spherical  dilatation  and  a  pure 
shear  deviation.  The  components  of  this  local  motion  Involving  rotation 
and  translation  do  not  produce  deformation  of  the  element,  hence  con¬ 
stitutive  equations  can  best  be  formulated  with  respect  to  a  rectangular 
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CARTESIAN  coordinate  reference  system  at  P  moving  such  that  the  local 
rotation  and  translation  components  vanish.  Such  a  coordinate  reference 
system  has  been  termed  a  'kinematically  preferred  coordinate  system1  by 
THOMAS  [20]  and  a  1 co-rotational  reference  system'  in  TRUESDELL  and 
TOUPIN  [21].  A  more  general  approach  was  employed  by  OLDROYD  [22] 
whereby  the  constitutive  equations  for  homogeneous  materials  are 
formulated  with  respect  to  a  'convected'  coordinate  system  which  moves 
and  deforms  with  the  material.  NOLL  [23]  in  his  earlier  work  on  in¬ 
variant  forms  of  constitutive  equations  uses  the  term  'principle  of 
isotropy  of  space'  and  RIVLIN  and  ERICKSEN  [24]  uses  the  term  'form  in¬ 
variance  under  rotation  of  the  physical  system'.  Since  physical  pro¬ 
perties  of  material  are  independent  of  the  coordinate  reference  system 
used,  the  object  of  the  above  works  was  to  determine,  within  the  frame¬ 
work  of  classical  mechanics,  quantities  associated  with  the  material 
which  are  independent  of  the  rigid  body  motions  of  this  material.  More 
generally,  OLDROYD  sought  invariant  forms  which  were  independent  of  any 
arbitrary  motion  of  the  material.  It  is  to  be  observed  that  the  less 
restrictive  'principle  of  material  objectivity'  mentioned  earlier  re¬ 
quires  that  the  constitutive  equations  be  invariant  under  any  arbitrary 
change  of  the  frame  of  reference  and  hence  of  the  observer. 

In  SECTION  2.1  of  this  chapter,  physical  quantities  are  de¬ 
termined  which  are  invariant  under  arbitrary  rigid  body  motions  con¬ 
sidered  as  elements  of  the  proper  orthogonal  transformation  group. 

Such  quantities  are  called  objective  and  attention  will  be  confined 
only  to  those  which  are  needed  in  the  formulation  of  the  constitutive 
equations  of  the  PRANDTL-REUSS  theory  for  finite  strain.  Whether  or  not 
the  rigid  body  motions  are  subjected  to  the  material  body  and  the  applied 
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forces  or  to  the  coordinate  reference  system  is  of  no  consequence  under 
this  approach.  This  would  not  be  true,  however,  if  the  rigid  body 
motions  were  considered  as  elements  of  the  full  orthogonal  transfor¬ 
mation  group.  Reflection  would  then  be  a  permissible  operation  but 
this  is  physically  meaningless  if  applied  to  material  bodies.  All  co¬ 
ordinate  reference  systems  used  are  rectangular  CARTESIAN  systems  and 
time  is  the  absolute  time  of  classical  mechanics.  Thus  the  approach 
taken  here  is  the  same  as  that  considered  in  THOMAS  [20]  and  ERINGEN 
[25]  but  differs  from  that  in  GREEN  and  ADKINS  [26]  and  NOLL  [27]. 

2.1  OBJECTIVE  TENSORS  AND  STRESS  RATES 

Consider  a  deforming  material  body  moving  through  space  and 
let  x.j  and  >T^(i  =  1,2,3)  be  the  spatial  coordinates  of  the  same  ma¬ 
terial  point  with  respect  to  two  coordinate  reference  frames  S  and  J. 

Then  the  two  motions  -j ,  ^2*  ^3*  ^  ^  =  ^-j  ( C  ^ »  t )  and  x^  (C-j,  ^2*  ^3*^ 
x.  (£.. ,  t)  are  said  to  be  objectively  equivalent  if  and  only  if 

x.  (£.,  t)  =  aij.(t)  x.  (£.,t)  +  bi(t) 


where 


(2.1 .la,b) 


aij(t)  aik^  '  aji^  aki^  '  Sjk  ’  laij(t)l  +  1 

and  T  =  t  -  a,  a  being  a  constant.  The  absolute  times  t  and  t  used 
with  the  S  and  S  systems  respectively  can  have  different  fiducial  time 
origins.  The  £.  (i  =  1,2,3)  are  the  LAGRANGIAN ,  convected  or  material 
coordinates  of  the  material  particle  and  may  be  taken  as  the  coordinates 


30 


of  the  particle  in  either  S  or  S  at  some  specific  instant  of  time  pro¬ 
vided  a  is  known.  Any  tensor  associated  with  the  material  particle  is 
also  said  to  be  objective  if  in  any  two  objectively  equivalent  motions 
it  obeys  its  appropriate  tensor  transformation  law  for  all  times. 

From  equation  (2.1.1a),  one  obtains 


V^t) 


Dx.(t) 

Dt 


aij(t) 


Vj(t)  +  X j ( t ) 


Dag(t)  Db.(t) 
“Dt  +  “Dt 


(2.1.2) 


where  and  ^  denote  differentiation  with  respect  to  time  keeping  the 
convected  coordinates  £ .  constant.  Equation  (2.1.2)  shows  that  velocity 
is  not  objective  and  it  follows  likewise  using  equation  (2.1.2)  that 
acceleration  is  not  objective.  Differentiation  of  equation  (2.1.2)  with 
respect  to  >T. (t)  gives 

J 


vi,j(x)  =  aik(t)  vk,m(t)  -= 


3xm(t)  Da i k ( t )  3xk(t) 


3x.(x) 


Dt 


3x.(t) 


(2.1.3) 


Since 


xk(t)  =  aj>k(t)  x^(t)  -  b.{t)  a  j  k  ( t ) , 


equation  (2.1.3)  is  written  as 


Vi>j(T)  =  a.k(t)  ajm(t)  vk>m(t)  +  a.k(t) 


Daik(t) 

Dt 


(2.1.4) 


from  which 


(.*)  ,X6 
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Da ..  (t) 


vj,i(T)  ■  #1kW  ajn1(t)  Vk(t)  +  aik(t>  ~ M  •  (2-1-5) 


Moreover,  from  equation  (2.1.1b), 


Da.,  (t) 


=  -  a . 


Da  -j  k  ( t ) 


jk  Dt 


(2.1.7) 


and  on  defining 


Da ..  (t) 

“ij^  s  aik^  Dt 


equations  (2.1.4)  and  (2.1.5)  are  written  as 


vi,j(T)  =  aik(t)  ajm(t)  vk,m(t>  ' 


and 


(2.1 ,7a, b) 


vj,i<T)  =  aik(t)  ajm(t)  vm,k(t)  + 

respectively.  From  equations  (2.1.7a,b),  one  then  obtains  the  ZORAWSKI 
relations  [28] 


d1j<T>  =  aik(t)  ajm(t)  dkm(t) 


“i ‘  aik^  ajm(t)  “km^  '  “ij^  ’ 


and 


(2.1 ,8a,b) 


6  nsrt  n a  .  (c .  f  .S'  b<  6  (£.  f .  a  if  ipe 
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where  are  components  of  the  vorticity  tensor  defined  by 


l  3vi  3vi 

“ij  "  2  vax.  "  9x .  • 

J  1 


(2.1.9) 


Thus  the  strain  rate  tensor  is  objective  while  the  vorticity  tensor 
is  not. 

With  respect  to  the  S  system,  the  stress  components  at  a 
material  point  P  of  the  stressed  medium  at  time  t  are  denoted  by  a. .(t) 

'  J 

and  with  respect  to  the  J  system  they  are  denoted  by  a.  .(t)  ,  t  =  t  -  a. 

*  J 

The  spatial  coordinate  dependence  of  the  stress  tensors  is  understood 
in  the  notation.  Due  to  the  objective  tensor  character  of  stress, 


a .  .  (t)  =  a.  (t)  a.  (t)  cr  (t)  . 
ij  '  lm  jn  mn 


(2.1.10) 


Taking  the  material  derivative  of  both  sides  of  equation  (2.1.10), 
one  obtains 


Do. . (t)  Do  (t)  Da.  (t) 

— Dr  =  aim*t^  ajn^^  — Dt  +  aim^  ~ Dt  amn^' 


lm 


Da.  (t) 


+  ajn^^  Dt  amn^^ 


(2.1.11) 


which  shows  that  the  material  derivative  of  the  stress  tensor  is  not 
objective.  However,  from  equation  (2. 1 . 1 1 ), objective  quantities  can  be 
obtained.  For  from  equations  (2.1.4)  and  (2.1.5), 
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=  -  a.n(t)  v..(t)  +  ad.(t)  vn>.(t). 
Addition  of  equations  (2-1 .1 2a ,b)  and  use  of  equation  (2.1. 

Da.  (t) 

Hit  =  ain^  '  aji^  “in^^ 

Substitution  of  equations  (2.1.12a),  (2.1.12b)  and  (2.1.13) 
into  equation  (2.1.11)  gives 


ft  °ij(x)  =  air(t)  ajs(t)  ft  °rs(t)  ’ 
=  a.r(t)  a.s(t)  ars(t)  . 


and 


Vo.  .(t) 
J - 

Vt 


-  a.r(t)  ajs 


Vo  At) 

(t)  - - - 

[Z)  Vt 


where 


-  Da.  •  (t) 

4ra..(t)  =  — H -  -  v,  m(t)  a,m(t)  -  v,  m(t)  a„ 


6t  ij 


Dt  J  »m  inr  #  i  ,nr  '  mj 


Da,  ,(t) 

/  +  \ - LJ_ — +  v  .(t)  a.  (t)  +  v  .(t)  a  .(t), 

aijlw  =  Dt  m,jv  inr  '  m,iv  mjv 


(2.1.12b) 

9)  yields 

(2.1.13) 

alternately 

(2.1 .14a,b,c) 


(t), 

(2.1 .15a,b,c) 


and 


( 3 ) .  oG 
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Vo.At)  Da..(t) 

~Pt  =  Dt  +  CTim(t>  “mj(t)  +  CTjm(t)  • 

6aj .  „  Vo. . 

Thus  gy-  ,  a^.  and  are  components  of  objective  tensors  called 
respectively  the  GREEN-OLDROYD  [29],  the  COTTER-RI VLIN  [30]  and  the 
JAUMANN  [31]  stress  rate  tensors.  There  is  also  the  TRUESDELL  [32] 
stress  rate  tensor  with  components 


- +  v  (t)  o..(t)  -  v.  (t)  a.  (t)  -  v.  (t)  a  .(t) 

Dt  m,mv  '  '  j  ,mv  '  imv  '  i  ,rrr  '  mjv  ' 

(2.1.16) 


formed  from  equation  (2.1.15a)  by  addition  of  expression  v  (t)  a..(t). 
In  fact  from  any  of  the  stress  rates  defined  above  others  can  be  formed 
by  addition  of  dimensionally  correct  terms  which  involve  components  of 
the  strain  rate  tensor.  If  the  stressed  medium  and  applied  forces 
undergo  a  rigid  body  motion  (d.  .(t)  =  0),  all  the  stress  rates  tensors 

*  vJ 

above  reduce  to  the  JAUMANN  stress  rate  tensor.  Moreover,  if  S  is  a 
co-rotational  reference  frame  (o>.  .(t)  =  0,  v.(t)  =  0)  and  the  stressed 

•  J  I 

medium  and  applied  forces  undergo  a  rigid  body  motion  (d. . ( t )  =  0)  and 

*  J 

if  the  stress  field  a.,  is  independent  of  time  with  respect  to  S,  then 

■  J 

the  JAUMANN  stress  rate  tensor  vanishes.  Being  objective  tensors,  all 
the  stress  rate  tensors  mentioned  vanish  with  respect  to  objectively 
equivalent  reference  frames  for  rigid  body  motions  of  the  stressed  ma¬ 
terial  whose  stress  field  is  independent  of  time  when  referred  to  a 
co-rotational  reference  frame. 

In  view  of  the  lack  of  a  unique  objective  stress  rate  tensor. 
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PRAGER  [33]  considered  the  implication  of  using  the  various  derivatives 
defined  by  equations  (2.1 ,15a, b,c)  and  equation  (2.1.16)  in  the  con¬ 
stitutive  equations  of  the  theory  of  plasticity.  PRAGER  concluded  that 
the  JAUMANN  derivative  defined  by  equation  (2.1.15c)  is  the  most  suitable 
for  use  in  the  constitutive  equations  of  plasticity  since  for  this  de¬ 
rivative  only  does  zero  stress  rate  imply  stationary  stress  invariants. 


2.2  CONVECTED  TIME  DIFFERENTIATION 

The  stress  rate  tensors  (2. 1 . 1 5a ,b ,c)  can  be  derived  by  a 
method  attributed  to  CAUCHY  [34]  and  generalized  by  OLDROYD  [22].  Con- 
vected  differentiation  with  respect  to  time  t  of  a  tensor  intrinsi¬ 
cally  associated  with  the  deforming  material  is  introduced.  This 
operation  involves  no  dependence  on  a  fixed  frame  of  reference  or  on 
the  motion  of  the  material  in  space. 

Restricting  discussion  to  a  tensor  T  with  components  T™1] 

referred  to  a  fixed  curvilinear  coordinate  system  x1  (i  =  1,2,3),  the 

(ST 

convected  time  derivative  of  T  is  defined  as  that  tensor  under  all 
coordinate  transformations  which  reduces  to  jjjr  (or  — )  in  any  convected 
coordinate  system  £.(i  =  1,2,3)  where  7  denotes  the  tensor  T  in  the 
convected  system  and  T ^  denotes  its  components  with  respect  to  this 
system.  T  is  assumed  to  be  a  tensor  under  all  coordinate  transformations, 
including  transformations  to  convected  coordinates.  Thus 


=<x6  3x^  3x^  =  Tmn  9x1  9x1 

y6  g £-ot  lj  ^ 


(2.2.1) 


Differentiating  equation  (2.2.1)  with  respect  to  time  t  holding  con¬ 
vected  coordinates  constant,  i ndicated  by  ,  and  using 
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one  obtains 


D_  , 9x\  =  3v^  9x^ 
Dt  3£a  3xn  9£a 


DTUV 


_  .  .  .  nTmn 

y6  _  9£  9£  Sx1  9xJ  r  l  j  Tmn  9yP 

Dt  „  m  „  n  1  Dt  1  D.i  I  i 


9xm  3xn  9£Y  9£6  Dt 


PJ 


9x 


+  Tmn  2^  _  jPn 


nP  9xJ 


1J  9xP 


Tmp  9vn , 
ij  9xP  ’ 


(2.2.2) 


Since  are  functions  of  x1  and  t, 

■  J 


DTmn 


9Tmn 


2 1  =  _2 _L  + 
Dt  9t 


9T 


mn 

ii  vp 


3x' 


Thus,  from  equation  (2.2.2), it  is  seen  that  the  components 


6  Tmn 
6t  1  j 


3Tmn 


9t 


LI  + 


9Tmn 


22  VP 


+  T 


mn  9v 


9x' 


PJ  9X1 


-  +  T 


mn  9v* 


^  3xJ 


-  T 


pn  9vm  _  -j-mp  9v_n 


^  9xP 


1J  9xP 


d?™1? 


reduce  to  — under  transformations  from  the  fixed  spatial  coordinate 

6  -rmn 


Dt 


system  to  a  convected  coordinate  system.  The  components  T..  .  are 
components  of  an  absolute  tensor  since  they  are  identically  equal  to  the 
tensor  components 
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9Tmn 

_ tj_  ,  Tmn 

3t  i j  ,p 


+  T 


mn 

PJ 


+  Tmn 
IP 


v 


P 

J 


-  Tpn  vm  -  Tmp 
1J  ,P  TJ 


V 


n 

»P 


where  subscript  commas  denote  covariant  differentiation. 

For  the  contravariant  stress  components  omn, 


6  mn  9amn  mn  p  pn  m  mp  n 
fit  Q  9t  a,p  V  "  a  v,p  '  a  vs 


(2.2.4) 


which  with  respect  to  a  rectangular  CARTESIAN  coordinate  system  are 


6  mn  ^amn 

o  j_  CT  pk  j.  ™  O  V  “  O  V 

ot  Dt  pn  m,p  mp  n,p, 


the  components  of  the  6REEN-0LDR0YD  objective  stress  rate  tensor.  Simi¬ 
larly,  for  the  covariant  stress  components  amn  , 


fit  amn 


9a 


mn 


9t 


a  vp  +  a  vp  +  a  vp  (2.2.5) 
mn,p  pn  ,m  mp  ,n 


which  with  respect  to  a  rectangular  CARTESIAN  coordinate  system  are 


fit  amn 


Dt 


v  +  a 
p  ,m  mp 


5 


the  components  of  the  COTTER-RIVLIN  objective  stress  rate  tensor. 


j 


38 

*"et  9mn  and  9mn  be  resPectively  the  covariant  and  contravariant 
components  of  the  metric  tensor  for  a  fixed  spatial  curvilinear  co¬ 
ordinate  system.  Then  since  Dg  /Dt  =  0  and  Dqmn/Dt  =  0, 

mn  3  '  * 


tt  g  =  2d 
6t  3mn  mn 


(2.2.6a,b,c) 


and 


5  npv  _  np  6  „ 

6t  ^9mn  0  6t  9mn 


„  6  np 

9mn  6t  0 


From  equations  (2.2.6b,c)  it  follows  that 


9mn  It  °nP  =  It  °"P>  ’  ^  > 


mn 


showing  that  the  raising  and  lowering  of  indices  is  not  commutative 
with  the  convected  time  derivative  operator  ^  except  when  the  stressed 
material  executes  only  a  rigid  body  motion  (d.  .  =  0)  at  P.  On  trans- 
forming  equations  (2.2.4)  and  (2.2.5)  to  a  fixed  spatial  curvilinear 
coordinate  system  such  that  d..  =  0  instantaneously  with  respect  to  this 

*  vJ 

system,  one  gets 


6  mn  3  mn  ,  mn  ,  p  pn  m  mp  n  , 
tt  o  =  ttO  +  a  vk  -  cr  co  -  a  K  co 
6t  3t  ,p  .p  .p 


and 


(2.2. 7a, b) 
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where 


and  the  subscript  commas  denote  covariant  differentiation.  Relation¬ 
ships  (2.2.7a,b)  define  components  of  absolute  tensors  which  reduce  to 
the  components  of  the  JAUMANN  stress  rate  tensor  given  by  definition 


(2.1.15c)  on  choosing  the  spatial  coordinate  system  as  rectangular 


CARTESIAN.  The  JAUMANN  derivative  thus  measures  the  time  rate  of  change 
with  respect  to  a  set  of  coordinate  axes  instantaneously  rotating  about 
itself  with  the  angular  velocity  of  the  material.  For  this  reason, 
the  JAUMANN  derivative  is  referred  to  as  the  co-rotational  time  derivative. 
Historically,  the  co-rotational  time  derivative  was  first  introduced  by 
ZAREMBA  [35].  In  this  thesis,  however,  it  will  be  called  the  JAUMANN 
derivative  as  it  is  in  current  plasticity  literature.  Since  the  JAUMANN 
derivative  does  not  depend  upon  the  components  of  the  strain  rate  tensor, 
it  is  the  simplest  of  all  the  objective  time  derivatives.  Moreover, 
it  is  the  only  objective  time  derivative  which  is  commutative  with  the 
operation  of  raising  and  lowering  indices. 


For  purposes  of  reference,  relationship  (2.2.7a)  is  rewritten 


as 


i  k  i  k  k  i 
+  v  +  ak  u>j  -  o.  , 


(2.2.8) 


Vt  9t 


which  upon  expanding  the  covariant  derivatives  becomes 
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with 
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1  =  1 
jk  "  2 


/9gmj  9gkm 
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the  CHRISTOFFEL  symbols  of  the  second  kind. 


2.3  PRANDTL.-REUSS  EQUATIONS  AT  FINITE  STRAIN 

Following  THOMAS  [36],  the  PRANDTL-REUSS  equations  for  plastic 
regions  valid  for  finite  strain  are  obtained  by  insuring  that  the  classi¬ 
cal  PRANDTL-REUSS  equations  do  not  violate  the  'principle  of  material 
objectivity1.  This  is  accomplished  simply  by  replacing  the  non-objective 
material  derivative  of  stress  deviation  with  the  JAUMANN  derivative  of 
stress  deviation.  The  strain  rate  tensor  is  not  replaced  since  it  is 
a  measure  of  the  instantaneous  rate  of  change  of  lengths  and  angles  of 
material  elements  in  the  deforming  material  [37]  and  is  an  objective 
tensor.  Thus  from  equations  (1.5.4a,b),  together  with 


Pp  Dp 
Vt  Dt  ’ 


one  obtains 
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(2.3. la, b) 
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the  constitutive  equations  for  a  compressible  elastic-perfectly  plastic 
MISES  solid  in  finite  strain.  These  reduce  to 


J__iL 

2  y  Vt 


Vs. 


,  A  = 


A" 

n  m 


3 


and 


(2.3.2a ,b) 


d 


k 

k 


=  0 


for  an  incompressible  elastic-perfectly  plastic  MISES  solid.  Similarly 
from  equations  (1.5.8a,b),  one  obtains 


I  Vs .  .  in  3f 

J _ LI  .  J_  6  +  y  i  _r_ 

2y  Vt  3K  Dt  °ij  I  Ar  3s.(  .  ’ 


and 


(2.3.3a,b) 


k  =  1  D£ 

k  "  K  Dt  » 


the  constitutive  equations  for  a  compressible  elastic-perfectly  plastic 
TRESCA  solid  in  the  plastic  range.  Equations  (2.3.3a ,b)  are  also  valid 
for  an  elastic-perfectly  plastic  solid  when  the  yield  surface  has  any 
finite  number  of  surfaces  with  piecewise  continuously  turning  normals. 

For  elastic  regions,  the  constitutive  equations  valid  for 
finite  strain  are  obtained  by  employing  proper  objective  tensors  in 
equations  (1.2.3a,b).  This  gives 
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Vs .  . 

-pi1  ■  fij 


a 


and 


(2.3.4a,b) 


the  PRANDTL-REUSS  equations  for  an  elastic  solid. 

By  adopting  the  constitutive  equations  for  a  hypo-elastic 
solid,  A.E.  GREEN  [38]  also  established  the  PRANDTL-REUSS  equations  at 
finite  strain  for  a  compressible  elastic-perfectly  plastic  MISES  solid. 
GREEN  employed  the  GREEN-OLDROYD  stress  rate  tensor  and  derived  the 
constitutive  equations 


6_ 

6t 


im 


m 


» 


and 


(2.3.5a,b) 
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i  2 

for  the  non-plastic  region  where  s.  <  2k  ,  and 


f  s 

6  eiJ  +  f1  c 1  m  fJ  =  fij  _  n  m  ciJ 
6t  5  S  m  "  S  Tm  T  2k2  S 


(2.3.6a) 


and 
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ft  P  =  -  K  d£  ,  (2.3.6b) 

O  ©  9  • 

for  the  plastic  region  where  s1  s'?  =  2k^  and  d1s'?  >  0,  If  the  MISES 

3  "i  J 

solid  is  incompressible,  equations  (2.3.5a,b)  and  (2.3.6a,b)  reduce 

to  equations  (2.3. 4a, b)  and  (2. 3,2a, b)  respectively.  A  discussion  of 
the  relationship  between  a  hypo-elastic  material  of  grade  two  and  an 
elastic-perfectly  plastic  MISES  solid  is  found  in  CHAPTER  V. 

To  eliminate  the  use  of  separate  sets  of  constitutive  equations 
in  both  the  elastic  and  plastic  regions,  THOMAS  [18]  proposed  a  single 
set  of  constitutive  equations  valid  for  both  regions.  With  the  as¬ 
sumption  that  the  constitutive  equations  have  the  form 


Vo.  . 


_  +  Bmn  o  , 

ij  mn  ij  mn 


with  A™1?  and  B11.".1  isotropic  tensor  functions  of  the  stress  deviation  in- 
1 3  i  J 

variants  and  strain  rate  deviation  invariants,  together  with  plastic  in¬ 
compressibility  and  the  VON  MISES  yield  criterion,  THOMAS  established  the 
constitutive  equations 


,mn 
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(2. 3. 7a ,b) 


h  d 


smsn 

1  -  (A21) 


+  b  a 


2k 


44 


In  general  y,  h,  b  are  invariants  dependent  on  position  and  time  but 
for  homogeneous  material,  they  must  be  material  constants.  For  plastic 
flow,  equations  (2.3.7a,b)  are  such  as  to  reduce  to  equations  (2.3.2a,b). 


2.4  PRANDTL-REUSS  CONSTITUTIVE  EQUATIONS  AND  FIELD  EQUATIONS  IN 

CURVILINEAR  COORDINATES  FOR  AN  INCOMPRESSIBLE  ELASTIC- 

PERFECTLY  PLASTIC  MISES  SOLID  IN  FINITE  STRAIN 

Let  x^i  =  1,2,3)  be  a  fixed  curvilinear  coordinate  system 


of  the  metric  tensor.  Then  the  governing  equations  for  the  plastic  flow 
of  an  incompressible  elastic-perfectly  plastic  MISES  solid  with  respect 
to  this  coordinate  system  with  commas  denoting  covariant  differentiation 
with  respect  to  the  metric  are: 

(i)  the  PRANDTL-REUSS  equations  (2.3.2a),  namely. 


Vs 


(2.4.1a) 


with 


and 


Vt 


Vs\ 


as  obtained  from  equation  (2.2.10)  on  replacing  oj  by  sj .  and  the 


equation  of  incompressibility 


d|  =  0  , 


(2.4.1b) 


where 


d1  =  —  gi v.  .  +  v .  ,  ) 
j  2  y  k,j  j,k' 


and  v1  are  the  contravariant  components  of  the  velocity  vector; 

(ii)  the  equations  of  motion 


•  •  •  • 


•  • 


a1^  +  pg1J  F.  =  p(|p  +  v'j  vJ) 


(2.4.2) 


where 


«1J  -  S1J  -  pg1J,  P  - 


1  /  rwu 
‘  3  ®mn  a  '  ’ 


p  is  the  density,  and  F.  the  covariant  components  of  the  body  force  per 
unit  mass; 

(iii)  the  equation  of  continuity 


H+<«Ai 


=  o 


(2.4.3) 


(iv)  the  VON  MISES  yield  criterion 


f  rdf  i  on  to  .  f:  ip  > 


and 
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(v)  the  identity 


s!  s  0  . 


(2.4.5) 


In  general,  the  governing  equations  must  be  used  to  deter¬ 
mine  eleven  unknowns;  namely,  the  six  components  s1J  of  the  symmetric 
stress  deviation  tensor,  the  three  velocity  components  v1 ,  the  density 
p,  and  the  hydrostatic  pressure  -p.  As  listed,  the  governing  equations 
consist  of  eleven  quasi-linear  partial  differential  equations  (2.4.1a), 
(2.4.1b),  (2.4.2)  and  (2.4.3)  and  two  finite  equations  (2.4.4)  and 
(2.4.5).  However,  equations  (2.4.1a)  are  not  independent  of  the  two 
finite  equations  since  on  taking  the  JAUMANN  derivative  of  both  sides 
of  equations  (2.4.4)  and  (2.4.5),  one  obtains 


Vs] 


sj  Vt  0  and  Vt  °* 


showing  that  there  are  only  four  independent  PRANDTL-REUSS  equations. 
Thus  there  are  available  nine  partial  differential  equations  and  two 
finite  equations  for  the  determination  of  the  eleven  unknowns. 

For  quasi-static,  incompressible  elastic-perfectly  plastic 
problems  in  the  absence  of  body  forces,  the  stress  and  velocity  field 
can  be  determined  with  ten  of  the  governing  equations;  namely,  equations 
(2.4.1a,b) ,  (2.4.4),  (2.4.5)  and  the  equations  of  equilibrium 


a1'?  =  0. 
9  J 
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The  density  p,  if  required,  is  determined  by  equation  (2.4.3).  The 
discussion  also  applies  if  the  flow  should  be  steady. 

If  in  the  general  case  the  unknowns  are  taken  as  the  six  in- 
•  • 

dependent  components  a1J  of  the  symmetric  stress  tensor,  the  three 
velocity  components  v1  and  the  density  p,  then  there  are  available  four 
independent  PRANDTL-REUSS  equations,  three  equilibrium  equations,  one 
incompressibility  equation,  the  continuity  equation  and  the  yield  condi¬ 
tion  for  the  determination  of  these  ten  unknowns. 

For  the  general  cases  considered  above,  THOMAS  [39]  remarks 
that,  as  an  alternative,  five  PRANDTL-REUSS  equations  may  be  used  pro¬ 
vided  the  VON  MISES  yield  criterion  is  not  used  directly  but  only  as  a 
boundary  condition  on  the  surface  S  enveloping  the  plastic  deforming 
region. 

In  APPENDIX  A  of  this  thesis,  the  PRANDTL-REUSS  equations 
(2.4.1a)  are  expressed  in  terms  of  the  physical  components  of  the  stress 
deviation  tensor  and  the  velocity  vector  with  respect  to  a  spherical 
polar  coordinate  system  (r,6,<j>)  and  a  cylindrical  polar  coordinate 
system  (r ,0 ,  z) . 
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CHAPTER  III 


THE  PLANE  RADIAL  FLOW  OF  AN  INCOMPRESSIBLE 

ELASTIC-PERFECTLY  PLASTIC  MISES  SOLID 


A  stress  field  for  plane  flow  of  an  isotropic  rigid-perfectly 
plastic  solid  through  a  converging  channel  with  rough  sides  was  ob¬ 
tained  by  NADAI  [40].  Later  HILL  [41]  obtained  a  kinematically  ad¬ 
missible  velocity  field  associated  with  this  stress  field.  The  stream¬ 
lines  of  this  velocity  field  are  radii  directed  through  the  virtual 
apex  of  the  channel  and  the  components  of  the  stress  deviation  tensor 
of  the  stress  field  are  independent  of  the  radial  distance  from  the 
virtual  apex.  In  this  chapter  an  investigation  is  undertaken  of  the 
plane  radial  flow  of  an  incompressible  elastic-perfectly  plastic  MISES 
solid  through  a  converging  channel  with  perfectly  rough  sides,  that  is, 
sides  on  which  the  frictional  stress  is  equal  to  the  shear  yield  stress. 
The  case  when  the  assumed  constant  frictional  stress  developed  on  the 
channel  walls  is  less  than  the  shear  yield  stress  is  also  considered. 

The  solution  of  NADAI  and  HILL  is  valid  for  a  rigid-perfectly 
plastic  solid  with  either  the  TRESCA  or  the  VON  MISES  yield  criterion 
provided  the  appropriate  value  of  the  shear  yield  stress  is  used.  How¬ 
ever,  this  investigation  is  restricted  to  a  MISES  solid,  that  is,  one 
whose  yield  is  governed  by  the  VON  MISES  yield  criterion. 

3.1  PLANE  RADIAL  FLOW  THROUGH  CONVERGING  CHANNEL 

Referring  to  Figure  2,  let  (r,  0,  z)  be  cylindrical  polar 
coordinates  and  let  the  channel  occupy  the  region  -  a  <_  9  <_  a  with 
the  flow  parallel  to  the  plane  z  =  0.  The  flow  is  steady  and  quasi- 
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static  with  stream-lines  that  are  radii  passing  through  the  virtual 
apex  of  the  channel.  It  is  assumed  that  the  VON  MISES  yield  criterion 
is  satisfied  throughout  the  field  and  no  unloading  occurs. 

With  reference  to  APPENDIX  A,  SECTION  A. 4,  the  governing 
equations  in  terms  of  the  non-zero  physical  components  of  the  stress 
deviation  sr,  sQ,  sr0,  sz,  velocity  vr  and  stress  ar,  aQ,  Tr0,  az  are 
the  PRANDTL-REUSS  equations 


(3.1 .la,b,c,d) 
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2y  Asz  ; 
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r  3r 


the  equation  of  incompressibility 


(3.1.2) 


the  equilibrium  equations 


?3.1.3a,b) 
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the  yield  criterion 


s 


2 

r 


(3.1.4) 


and  the  identity 


sr  +  se  +  sz  E  0  • 
Equation  (3.1.2)  requires  that 

v  =jaM 

r  r 


(3.1.5) 


and  this  gives  the  following  non-zero  components  of  strain  rate: 


d 


r 


1  g'(e) 


(3.1.6) 


Equations  (3.1.6)  show  that  ratios  of  the  strain  rate  components  do 
not  vary  along  a  radius.  The  considered  solid  being  elastically  and 
plastically  incompressible,  both  the  elastic  and  plastic  components  of 
velocity  satisfy  the  incompressibility  equation  (3.1.2)  and  the  uni¬ 
quely  determined  physical  components  of  the  elastic  and  plastic  strain 
rate  tensor  are  of  the  form  given  by  equations  (3.1.6).  Since 

/ikdi.(p) 

5 1  =  - J. - - 

J  />(?)  dn(p> 

n  m 

by  equation  (1.3.9b),  s ]  must  be  functions  of  0  alone.  It  follows 

J 

then,  from  equation  (3.1. Id),  that  Asz  =  0  and  since  X  is  non-zero 
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where  flow  occurs 
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sz  =  -  (sr+  se)  =  °- 


Substitution  of  s0  =  -  sr  in  the  yield  criterion  (3.1.4)  gives 


(3.1.7) 


Equations  (3.1. 3a , b )  and  (3.1.7)  are  also  valid  for  the  corresponding 
flow  of  a  rigid-perfectly  plastic  solid  and  the  stress  boundary 
conditions  Tr0  =  k  at  0  =  a  and  Tr0  =  k  at  0  =  -  a  are  the  same. 
Consequently. the  stress  field  is  identical  to  that  obtained  by  NADAI. 


3.2  DETERMINATION  OF  STRESS  FIELD 
■  ■  .  ■  .  1  ■  — 

NADAI ' S  approach  [40]  for  determining  the  stress  field  is  as 
follows.  The  stress  field  is  expressed  in  terms  of  the  angle  ip  between 
the  direction  of  the  algebraically  greatest  principal  stress  and  a  radius. 
The  angle  ip  is  a  function  of  0  only  and  it  has  the  same  sign  as  0  since 
the  friction  acts  so  as  to  oppose  the  relative  motion  or  tendency  for 
relative  motion  between  the  solid  and  the  channel  walls;  ip  is  zero  on  the 
axis  0  =  0  and  ranges  between  -  tt/4  and  tt/4.  The  yield  criterion  (3.1.7) 
is  parametrized  in  terms  of  \p  by  letting 


srQ  =  k  sin  2ip 


and 


(3.2. la, b) 


s  =  k  cos  2\p  . 
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From  the  equation  of  equilibrium  (3.1.3b), 


3a 
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=  -  2k  sin  2ip 


from  which 


a0  =  -  2k  J0  sin  2ip  de  +  h(r) 


and, hence  from  equation  (3.2.1b), 


ar  =  2k  cos  2ip  -  2k  /0  sin  2ip  de  +  h(r). 

o 


(3.2.2) 


Substitution  of  equations  (3.2.1a,b)  and  (3.2.2)  into  the  equation  of 
equilibrium  (3.1.3a)  gives 


2k  cos  2ip  +  2k  cos  2ip  =  -  r  ^  . 


Setting  each  side  of  this  equation  equal  to  a  constant  2ck  yields 


dh 

dr 


2ck 


and 


(3.2. 3a , b ) 


d9  _  cos  2\p 
dip  ”  c  -  cos  2ip 

Integration  of  these  differential  equations  yields 
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h(r)  =  -  2ck  £nr  +  Aq  , 


and 


9  =  -  t|<  +  c  tan'1  (/f-±-4-  tan  <p)  ,  (3.2.4) 

/7T7  cT 

where  Aq  is  an  arbitrary  constant  and,  since  ip  =  tt/4  on  6  =  a,  c  is 
given  by 


f7~- 7 


tan'1  (/|44 


\  IT  . 

)  =  4  +  <*• 


(3.2.5) 


Also,  c  varies  from  °°  to  1.1922  as  a  varies  from  0  to  ir/2.  Now 


\  sin  2t  d6  =  f  sin  2»  cC-°-S-Js-  ^  d* 


=  1  cos  2*  -  ±  +  §  An 


and,  hence 


p  =  -  0Z  =  -  }  (ar  +  oe) 

=  2k  [c£n  r  +  j  In  (c  -  cos  2\p)  +  A]  (3.2.6) 

where  A  is  an  arbitrary  constant.  Equations  (3.2. la, b)  and  (3.2.6) 
completely  determine  the  stress  field  to  within  an  arbitrary  hydro¬ 
static  pressure. 
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3.3  DETERMINATION  OF  VELOCITY  FIELD 

Elimination  of  X  between  equations  (3.1.1b)  and  (3.1.1c)  and 
substitution  of  equation  (3.1.5)  into  the  result  gives 


2ysr0 


2  2  ‘ 

s  .  -  us  +  s 

r0  M  r  r 


(3.3.1) 


Substitution  of  the  components  of  the  stress  deviation  given  by 
equations  (3 . 2 . 1  a ,b)  into  equation  (3.3.1)  gives 


g 1  ( 9 )  _  2y  sin  2i jj 
g(0)  "  k  -  y  cos 


(3.3.2) 


With  the  use  of  equation  (3.2.4b),  this  becomes 


dq  ( Q )  _  2  sin  2U>  cos  2\I>  diy 

g(6)  (n  -  cos  2ip)(c  -  cos  2\pj 


where  n  e  -  >  0  and  hence,  provided  cos  2\jj  t  n,  integration  gives 

I" 


-nr  cos2\j;  -  n|n  -jl/c-n 
9(e)  "  D  [(c  -  cos"2*F] 


(3.3.3) 


where  D  is  the  constant  of  integration.  For  n  =  0, 


"  c  -  cos 


which  is  the  expression  obtained  by  HILL  [41]  for  the  rigid-perfectly 
plastic  solid.  From  equations  (3.1.5)  and  (3.3.3), 


r  ;  :e  iz  .  t  *  '  )3  nor  •  r  V2C  u’ 
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n 


1 


..  _  D  r  cos  2\p  -  n  iC-n  _  _ o , 

v  =  -  LJ - 1 - SrJ  ,  n  /  cos  2\p  , 

(c  -  cos  2\p)c 


(3.3.4) 


where  c  is  given  by  equation  (3.2.6)  and  D  is  obtained  from  the  volume 
flow  per  unit  thickness  and  is  negative  since  v^  is  negative. 

The  scalar  invariant  X,  expressed  in  terms  of  the  physical 
components  of  the  strain  rate  tensor  and  stress  deviation  tensor,  is 


X  =  ^  (drsr  +  2  dresre  +  d9s6)  ' 


Substituting  equations  (3.1.6)  and  (3.1.8a)  gives 


x  =  _ L  9 ( Q.)..  [s  -  s  +  9 1  ( Q ).  s  ] 

A  2k2  Ls0  sr  gTeT  SreJ  ’ 


and  substituting  equations  (3.2.1a,b)  and  (3.3.2)  gives 


r  /I  -  n  cos  2^x 
A  "  kr  '  n  -  cos  2\p 


(3.3.5) 


For  plastic  flow,  X  is  positive. 

If  n  e  k/y  <  1,  then  according  to  equation  (3.3.4),  lim  v  =  0 

_1  6^±eo 

where  6  is  the  value  of  6  corresponding  to  2\p  =  cos  n.  In  the  region 
o 

-  0  <  0  <  6  ,  the  determined  velocity  field  and  the  stress  field  satisfy 

o  o 

the  governing  equations  and  X  >  0.  However,  in  the  regions  -  a  £  0  <  -  0Q 

and  0  <  0  <  a  the  velocity  field  and  the  stress  field  result  in  X  <  0. 

o  — 

To  obtain  a  positive  X  it  is  necessary  to  take  D  >  0,  but  then  the  friction 
shearing  stresses  =  k  at  0  =  a  and  =-kat@=-ado  not  oppose 


:  . 22  JLr-  1 )  :■  -  .  ; 
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the  flow.  Consequently,  the  velocity  field  given  by  equation  (3.3.4) 
is  not  admissible  in  the  regions  -  a  £  0  <  -  0Q  and  0Q  <  0  £  a  unless 
D  is  zero.  From  equations  (3.3.2)  and  (3.3.4),  one  obtains  that 

2n-c 

g'(e)  =  ~  0  sin  2^e  (cos  2f  -  n)c'n  ,  D  <  0  , 

(c  -  cos2i|;)c’n 

in  the  region  -  0Q  <  0  c  0Q  where  n  <  cos  2\p  £  1 .  Thus  for  c  <  2n, 

lim  g ' ( 0 )  =  0  and  the  velocity  profile  is  normal  to  the  radii 

0+±0„ 

o 

0  =  ±  0  .  For  c  =  2n,  lim  g '  ( 0 )  =  -  2D/n2  (1  -  n2)^2  >  0.  For 

0-*-±0 

o 

c  >  2n,  lim  g 1 ( 0 )  =  +  »  and  the  velocity  profile  is  tangential  to  the 
e+±6o 

radii  0  =  ±  0  .  Also,  infinite  values  of  d  „  result  as  the  radii  0  =  ±  0 
o  ru 

are  approached  from  the  region  -  0Q  <  0  <  0Q;  this  is  permissible  since 
the  material  considered  is  not  viscous.  Moreover,  the  parameter  X 
undergoes  an  infinite  discontinuity  at  0  =  ±  0Q  but  this  is  not  incon¬ 
sistent  with  the  constitutive  equations  (3.1.1). 

If  n  =  1,  equation  (3.3.5)  reduces  to 


from  which  it  follows  that  vr  =  0  for  -  a  £  0  £  a  is  the  only  admissible 
velocity  field. 

If  n  >  1,  then  for  X  >  0  the  velocity  field  is  such  that  the 

radial  velocity  decreases  as  0  varies  laterally  along  a  0  -  line  from 

0  =  ±  0  to  0  =  0;  0  is  the  value  of  0  corresponding  to  2\p  =  cos"1  -  . 

0  0  n 

This  is  physically  unlikely;  moreover,  since  all  common  elastic-plastic 
materials  have  small  values  of  n  <  1 ,  the  phenomenon  is  of  no  great 
practical  importance. 


ftt  Jertt  rtoua  z\  b[9rt  ^taolsv  9ftt  0  <  A  n9rlJ  ,  <  n  ii 

■  iuo  1  •  r.*- o  nq  . 


3.4  CONCLUSION 
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For  the  region  -  0Q  <  6  c  e  ,  the  solution  to  the  converging 
radial  flow  problem  involving  a  channel  with  perfectly  rough  sides  is 
provided  by  the  velocity  equation  (3.3.4)  and  by  equations  (3.2.1a,b) 
and  equation  (3.2.7)  which  determine  the  stress  field  to  within  an 
arbitrary  hydrostatic  pressure.  For  the  regions  -  a  d  <_  -  0Q  and 
0Q  <  6  <  a,  the  solution  is  provided  by  vr  =  0  and  the  same  stress 
equations.  These  solutions  indicate  that  the  velocity  v^  is  zero  at 
the  channel  walls  whereas  the  velocity  field  obtained  by  HILL  for  the 
corresponding  problem  involving  a  rigid-perfectly  plastic  solid 
indicates  that  there  is  slip  at  the  channel  walls.  If  n  «  1,  the 
angle  a  -  0Q  is  small  compared  with  a  as  shown  by  the  values  in  the 
following  table.  These  values  are  for  the  particular  case  of 
a  =  24° 1 7 1  (c  =  2). 


TABLE  1 

VARIATION  OF  a  -  0  WITH  n  =  —• 


[a  =  24°1 7 ' 

i — i 

C\i 

II 

a 

n 

a  -  0 

0 

0.001 

0°0'4.511 

0.01 

0°0' 12.6" 

0.05 

O 

o 

O 

0.1 

0°4'40" 

Figure  3  shows  g(0)  for  an  elastic-perfectly  plastic  solid 
with  n  =  0.1  and  D  =  -1  and  also  for  an  rigid  plastic  solid  for 
a  =  24°  1 7 '  (c  =  2)  and  D  =  -1.  For  smaller  values  of  n,  of  the  order 
of  10-3,  which  are  realistic  for  metals,  the  difference  between  the 


[(S  =  a)  'U°*S  = 
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CD 


FIGURE  3  -  GRAPH  OF  g(e)  FOR 
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velocity  profiles  for  the  elastic-perfectly  plastic  and  rigid-perfectly 
plastic  solid  is  negligible  except  near  the  sides. 


From  the  present  solution  involving  the  elastic-perfectly 


plastic  solid,  it  is  possible  to  find  the  solution  when  the  constant 
frictional  stress  on  the  channel  walls  is  mk  (0  <_m  <_  1).  The  stress 
distribution  within  a  sector  ±  6(6  _<  a)  is  the  stress  field  for  a 
channel  of  angle  20  and  a  constant  frictional  stress  k  sin  2$.  Thus 
to  find  the  solution  for  a  channel  of  angle  20  and  frictional  stress 
mk,  it  is  required  to  determine  c  from  equation  (3.2.4)  with 

=  y  sin'^m.  Then  the  stress  field  is  given  by  equations  (3.2.1a,b) 
and  (3.2.6)  and  the  velocity  field  by  equation  (3.3.4)  provided  m  is 
such  that  0  <  m  <  A  -  n2,  n  <  1,  in  which  case  there  is  slip  along  the 


there  is  developed  a 


channel  walls. 


non-deforming  region  and  there  is  no  slip  along  the  channel  walls.  In 
the  deforming  region,  the  velocity  field  is  given  by  equation  (3.3.4). 


■ 

8f  m  b9bfvo-.q  (».e.e)  norJ6UP9  ti  bf»W  'iJfsofov  srtJ  bns  (8.  bns 

(K£.£)  noiJsupa  *d  nsvrp  et  bfsrt  y?  sofsv  sril  .norgsi  '«  '°'ot'  9i’ 


CHAPTER  IV 


ELASTIC-PERFECTLY  PLASTIC  FLOW  THROUGH 

A  CONVERGING  CONICAL  CHANNEL 

SHIELD  [42]  considered  the  axial ly-symmetric  flow  of  a  rigid- 
perfectly  plastic  material  forced  through  a  rigid  conical  channel,  out¬ 
lined  the  method  of  solution  for  a  general  yield  criterion  and  gave 
solutions  using  the  VON  MISES  and  TRESCA  yield  criteria.  These  solu¬ 
tions  involve  streamlines  that  are  radii  passing  through  the  virtual 
apex  of  the  cone  and  a  constant  friction  stress  between  the  channel 
wall  and  the  material  is  assumed.  The  corresponding  problem  with  an 
incompressible  elastic-perfectly  plastic  MISES  solid  is  considered  in 
this  chapter. 

Axially  symmetric  flow  through  the  converging  channel  is  con¬ 
sidered  and  the  channel  is  assumed  sufficiently  long  that  the  inlet  and 
exit  effects  can  be  neglected.  It  is  further  assumed  that  the  flow  is 
steady  and  quasi-static  and  that  all  points  in  the  stress  field  satisfy 
the  VON  MISES  yield  criterion  and  that  no  unloading  is  taking  place. 

A  solution  with  streamlines  that  are  radii  directed  through  the  virtual 
apex  of  the  channel  is  sought. 

4.1  GOVERNING  EQUATIONS  FOR  RADIAL  FLOW 

Let  (r,6,<|>)  be  spherical  polar  coordinates  with  the  axis  of 
the  channel  given  by  0  =  0  and  the  conical  surface  by  6  =  a  as  illustrated 
in  Figure  4.  It  is  assumed  that  the  friction  stress  acting  on  the  ma¬ 
terial  at  the  surface  of  the  channel  6  =  a  is  TrQ  =  mk  where  m  is  a 
constant.  Since  k  is  the  shear  yield  stress  and  Tre  is  non-negative 
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FIGURE  4  -  SPHERICAL  POLAR  COORDINATE  SYSTEM  FOR 

CONVERGING  FLOW  THROUGH  A  CONICAL  CHANNEL 

« 
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for  converging  flow,  m  must  satisfy  0  £  m  <  1 . 

With  reference  to  APPENDIX  A,  SECTION  A. 3,  the  governing 
equations  for  axially  symmetric  flow,  with  the  assumption  that  the 
velocity  component  v.  is  zero,  are 
the  PRANDTL-REUSS  equations 

v  !ir  fre!V  .  .  r!V  le , 

vr  3r  r  30  2y  ^3r  “  Xsr^  * 


8s0  s  0  3v  v 

vr  3r~  +  ~r~  30“  2y  "  As0], 


(4.1 . la ,b,c,d) 


asre  ,  sr  -  s0  9vr 


■1  9vr 


V  3  r  +  2r  36  2y  ^2r  30  ,wr0 


-  As„J  , 


3s . 

vr  37"  =  2yE“  "  Xs^  » 


where  sr>  sQ,  s^  and  srQ  are  the  non-zero  physical  components  of  the 
stress  deviation  and  vr  is  the  radial  component  of  velocity; 

the  equation  of  incompressibility 


3v 


r 


3r 


+ 


(4.1.2) 


the  equilibrium  equations 


3a._  i  3t 


— -  +  —  +  -  (2a  -  a  -  a  +  t  cot  0)  =  0, 

r  30  r  '  r  6  <j>  r0  ' 


(4.1.3a ,b) 


3x  .  i 


+  1  — +  -  [(a  -  a  )  cot  0  +  3t  fl]  =  0; 
8r  r  30  r  LV  6  r  r0J 


and  the  VON  MISES  condition 
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(ar-a0)2  +  (<V<V2  +  (°<t>  '  °r)2  +  6Tr02  =  6I<2’  (4-1'4) 


where  c^,  cQ,  and  ir0  are  the  non-zero  physical  components  of  the 
stress  tensor.  Only  three  of  the  PRANDTL-REUSS  equations  are  indepen¬ 
dent  since  the  sum  of  equations  (4.1.1a),  (4.Llb)  and  (4.1. Id)  is 
identically  zero.  Consequently,  there  are  seven  independent  equations 
for  the  six  unknowns  sr>  sQ,  s^,  srQ,  p  and  vr-  The  prior  assumption 
that  vn  =  0  is  the  reason  for  the  additional  equation.  If  a  solution 

with  v  =  0  can  be  found  which  satisfies  the  governing  equations 

u 

(4.1.1)  to  (4. 1.4), the  prior  assumption  is  justified. 

Equation  (4.1.2)  requires  that 


(4.1.5) 


where  g ( 0 )  is  a  function  of  0  only  and  B  is  a  constant.  The  components 
of  the  strain  rate  tensor  are  then 


d 


r 


B 


g(e) 

3 

r 


re 


=  ii 
2 


(6) 


The  ratios  of  the  strain  rate  components  obtained  from  equation  (4.1.5) 
do  not  vary  along  a  radius.  Therefore,  following  the  same  argument 
given  in  the  two  dimensional  plane  problem,  the  components  of  the 
stress  deviation  tensor  can  be  shown  to  be  functions  of  6  only. 

Consequently,  eliminating  X  from  equations  (4.1.1a)  and 


(4.1.1c)  gives 


J-ne  nt  m?  .f  f  o  (ijio  2uiu6i  of*  or  ob 

3/1J  *o  ?• '!.  9  Idcnq  msTq  »“iio  .i:  rb  o  J  srl*  ni*  r  r v r 3 

pr,  rJ  6n  r  nr  Ts  ,  i Ht ns  pa  2noO 


65 


9' (e)  - 

iler  ■ 


4p  s 


r6 


[sr0  - 2  (2y  -  sr  +  s0)] 


and  from  equations  (4.1.1a)  and  (4.1.1b),  together  with  use 
(4.1.2), 


9l (Q)  =  rfr _ ^jh 

gTeT  sr0  Lsr  +  s0  J* 


Equating  the  right  hand  sides  of  equations  (4.1.6)  and  (4.1 
simplifying  gives 


sre2  =isrse  +  !sr2  -  s02  -  usr  '  2>JV 


Equation  (4.1.8)  may  be  rewritten  as 


6s  2  =  3 ( s r  +  sQ)2  -  3(sr  +  3s0)  s0  -  6y  (sr  +  2< 


Combining  equation  (4.1.9)  with  the  yield  criterion  (4.1.4) 


(sp  -  Sj2  +  (s  +  2se)2  +(2sr  +  se)2  +  3(sr  +  s0) 


r  "6 


-  3srs0  -  9s02  -  6u(sr  +  2s0)  =  6k2 


if  the  relationship 


s  +  sQ  +  s.  =  0 
r  0  4) 


(4.1.6) 

of  equation 

(4.1.7) 

7)  and 

(4.1.8) 

0) .  (4. 1 .9) 
gives 

(4.1.10) 


(4.1.11) 
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is  used.  Upon  expansion  of  each  term,  equation  (4.1.10)  simplifies 
down  to 


3sr2  -  2usr  -  2k2  =  (4y  -  3sr)  s0. 


Consequently, 


s 


6 


3sr2  -  2ysr  -  2k2 
-  3sr 


for  sr  / 


(4.1.12) 


and  from  this  and  equation  (4.1.11), 

2k2  ‘  2ysr  4u 

S(|)  =  4y  -  3sr  for  sr  ^  3  *  (4.1.13) 

Substituting  equation  (4.1.12)  in  equation  (4.1.9)  and  simplifying 
gives 


sre2  =  {-  9sr4  +  18ysr3  -  12y2sr2  +  15k2sr2  -  24yk2sr 

-  4k4  +  16u2k2H4u  jTs-}2  .  (4.1.14) 

The  right  hand  side  of  equation  (4.1.14)  can  be  factored  further  to 
give 

sr62  ={(3sr2  -  4k2) (k2  -  4y2  +  6ysr  - 


and  hence 


M  jt  yar  +  - 
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1/2 

sre  =  4i T-  3s  {(4k2  -  3sr2)(3sr2  -  6usr  +  4y2  -  k2)}  (4.1.15) 

where  the  positive  root  is  taken  since  srQ  is  non-negative  throughout 
the  field.  For  sr0  to  be  defined  it  is  necessary  that 

(4  -  3n2)(3n2  -  63n  +  432  -  1)  >  0 

Sr  u  2  2 

where  n  =  •£—  and  3  =  .  The  second  factor  3q  -  63n  +43  -  1  has 

discriminant  12  -  1232  <  0  if  3  >  1.  It  is  assumed  that  3  >  1  as  it 

appears  that  no  solution  to  the  flow  problem  is  possible  if  3  £  1. 

Consequently,  the  second  factor  has  no  zeros  and  since  it  has  positive 

2 

concavity  it  is  positive  for  all  n.  Thus  srQ  is  defined  if  4  -  3n  >_  0, 
that  is,  if 


a  condition  that  is  clearly  satisfied  because  of  the  yield  criterion. 

From  the  requirement  that  vr  <0,  equation  (4.1.2)  indicates  that 
9v 

— -  >  0.  Then  from  equation  (4.1.1a)  together  with  the  further  re- 
3r  -  3v 

quirements  that  sr0  >  0  and  <0,  one  deduces  that  sr  >  0.  Con- 

p 

sequently,  0  <  sr  £  —  k  and  because  of  the  assumption  of  axial  sym¬ 


s  =  —  k,  s  =  s  = 

r  /3  e  't’ 


"  /3  °n  ^  =  O’ 


metry , 


,r  >  a  Vr  9fd*raaoq  ar  mardonq  rtor^9f#  o*  nol-Ju  o  on  61  .  •  9CCl[ 


,0  <  Sn£  -  >  bsnrtsb  ei  ^  aurtT  ,n  ff«  irf  svrJfJ  <q  at  Ji  tfhwanoa 


jsrtJ  29i63fbrif  (S.f.P)  noUeupe  ,0  »  JertJ  Jnaiiwitupst  9(1}  mo-' 
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or  for  all  6  if  the  channel  is  perfectly  smooth,  resulting  in  spherically 
symmetric  flow. 

From  the  requirement  0  £  sr  <_ —  k,  it  follows  from  equation 
(4.1 .15)  that 


0  1  sr6  1  0  "  ~V)1/2  k  =  m'k. 
rw  43 


Thus  a  solution  cannot  be  obtained  if  the  prescribed  shearing  stress 
at  the  surface  is  greater  than  m'k.  It  is  only  for  the  special  case 
of  a  rigid-perfectly  plastic  material  (m '  =  1)  that  the  prescribed 
shearing  stress  at  the  surface  can  be  k. 


4.2  DETERMINATION  OF  STRESS  FIELD 

Integration  of  the  equation  of  equilibrium  (4.1.3b)  gives 


=  f(r)  ■  /6  {(Sfi  "  SJ  COt  0  +  3srA*de  (4.2.1) 


0 


0 


re- 


since  the  components  of  the  stress  deviation  are  functions  of  6  only. 
From  equation  (4.2.1)  an  expression  for  the  hydrostatic  pressure, 

1  k 

“  P  E  TQk  *  is  obtai’ned  and  is 


-  P 


=  f(r)  -  sn  -  ;9{(s„  -  s.)  cot  6  +  3s  }d6.  (4.2.2) 


■0  Jo  "  e  * 


re- 


The  equation  of  equilibrium  (4.1.3a)  gives 


r  =  .  (!iri 

r  9 r  V  90 


+  3sr  +  sr0  cote) 


ae 


(4.2.3) 


vu»i  ’  K  n\2 


29vh  (dE.r.fr)  rnuridrUup9  to  nol'J6up9  9rtJ  to  noHeipotfil 


( f .S.£) 
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since  Qr  =  sr  "  P  and  2 -  a0  -  =  3sr.  The  left  hand  side  of 

equation  (4.2.3)  is  a  function  of  r  alone  and  the  right  hand  side  is 
a  function  of  0  alone;  consequently, 


9R  = 
3 


ck 


(4.2.4) 


and 


ds 


re 


de 


+  3s„  +  s  .  cot  e  =  ck 
r  ry 


(4.2.5) 


where  c  is  constant.  Since  for  converging  flow  p  must  increase  as  r 
decreases,  c  is  positive. 

Integration  of  equation  (4.2.4)  gives 


p  =  -  ck  £n  (£-)  -  h(6) 

o 


where  rQ  is  a  constant  which  depends  on  the  specified  stress  at  some 
point  ( r ,6 )  in  the  field  and  from  equation  (4.2.2) 

h(6)  =  -  s0  -  /6  {(sQ  -  s  )  cot  G  +  3sr0}d6.  (4.2.6) 

o  ^ 

In  order  to  obtain  s  Q  as  a  function  of  0,  the  non-linear 

ro 

differential  equation  (4.2.5)  with  sr  related  to  sr0  by  equation 

(4.1.15)  must  be  solved  subject  to  the  conditions  sr0  =  0  on  6  =  0 

and  s  =  m'k  on  0  =  a.  The  constant  c  in  equation  (4.2.5)  is  deter- 
r0 

mined  by  a,  3  and  m' . 


-  !  >b  zi  .S  >)  norJsups  :ir  o  i  6*2100  -vti  .jo  =  6  .10  J*  -  ,^2-bfi 


Equation  (4.2.5)  is  written  in  the  non-dimensional  form 

T '  (6)  +t(6)  cot  0  +  3n(e)  =  c  (4.2.7) 

s  s 

where  t(0)  =  — and  n ( © )  =  jr-  .  For  the  limiting  case  y  +  «>  for 

fixed  k,  it  follows  from  equation  (4.1.5)  that  n(6)  —  [1  -  x2(0)]1//2. 

/3 

Thus  for  the  special  case  of  a  rigid-perfectly  plastic  material, 
equation  (4.2.7)  becomes 

T  1  (e )  +  x(0)  cot  0  +  2/3  [1  -  I2(0)]^2  =  c 

which  in  non-dimensional  form  is  the  equation  obtained  by  SHIELD  for 

the  rigid-perfectly  plastic  case. 

To  solve  equation  (4.2.7)  a  value  of  c  was  chosen  and  the 

equation  was  integrated  numerically  to  give  0  as  a  function  of  t  in 

the  interval  0  <_  x  £  m1  for  a  specified  3.  The  solution  was  started 

at  0  =  0,  x  =  0  and  continued  until  x  =  m'  was  attained;  the  value  of 

0  corresponding  to  the  prescribed  x  =  m'  at  the  surface  being  the 

cIt 

semi-angle  a  of  the  conical  channel.  This  required  that  be  known  on 
0=0.  From  equation  (4.2.7), 

(^-)  =1  i m  (c  -  3n  -  xcot0) 

dy0=O  0+0 

=  c  -  2/3  -  lim  xcot0. 

0+0 


(4.2.8) 
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But 


1  im  Tcote  =  (-ttt) 

e+o  de 


dr 

d6 


e=o 


and  thus,  from  equation  (4.2.8),  it  follows  that  (~) 


2 


Since  t  is  non-negative  in  the  field,  c  >_  2/3  with  equality  holding 
only  for  spherically  symmetric  flow.  The  results  of  the  numerical 
integration  are  shown  graphically  in  Figure  5  for  the  case  3=^-  =  10. 


The  numerical  solution  of  equation  (4.2.7)  involves  the 


numerical  determination  of  n  from  x  using  equation  (4.1.15)  and  sub¬ 
stitution  of  the  values  of  n  so  obtained  in  equations  (4.1.12)  and 


(4.1.13)  gives  the  corresponding  values  of  sQ  and  s^.  The  results  so 


obtained  are  shown  graphically  in  Figures  6  and  7  for  the  cases  3  =  10 
and  c  =  10  and  7  respectively. 

The  function  h(e)  defined  by  equation  (4.2.6)  is  also  shown 
graphically  in  Figure  7  for  3  =  10  and  c  =  10  and  7. 

4.3  DETERMINATION  OF  THE  VELOCITY  FIELD 

Substituting  equation  (4.1.12)  in  equation  (4.1.7)  gives 


v  =BaM  =  ^exp  t./ 

r  r  r  o 


Hence 


A  bne  OT  *  o  bms  Of  =  a  icrt  ^  nt 


FIGURE  5  -  FUNCTION  t(Q)  FOR  VARIOUS  VALUES  OF  C 
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■  ^  exp  C-  /9{— 2 

r  o  43 


4  -  3n2 
-  6Bn  +  3ri 


1/2  4B2  -  38n 

>  {^T-}de 


(4.3.1) 


=  \  exP  [-  A(e)], 
r 

where  B  is  a  constant  depending  on  the  volume  flow  and  is  negative 
since  v^,  is  negative.  For  the  special  case  of  the  rigid-perfectly 
plastic  MISES  material,  equation  (4.3.1)  becomes 


vr  =  K exp  2>/3  /e  ^  "  x2(e))^2  d0] 

r  o 


which  is  the  result  obtained  by  SHIELD.  The  integrand  in  equation 
(4.3.1)  has  a  singularity  at  n  =  ^  .  For  0  <  0  <  0Q,  where  0Q  is 
the  value  of  9  corresponding  to  n  =  1  .  the  integral 


1  0  9  l/?  9n  46  -  3 

A(6n)  >  1  [ ( 4S2  -  3) (B2  -  1)] '  J  (-  3  + 


6 


n  -  e 


-)de 


■  y(B)  9q  +  6(B)  / 


9 


o  dG 


0  n  -  J 


where  y(3)  and  6(3)  are  positive  numbers  for  a  specified  3  >  1.  Now 


0 


/ 


o  d0 


o  n  - 


3 


,  .  rV£  1  ,d9N 

/  - T  W  dr| 

e+0  o 


n  ‘  B 


i  •  r  /d0  \  f60  e  dll  1 

>  lim  J  T> 

"  e+0+  dn  0o-e  o  n  -  j 


9,00  !  {f  .'£.*)  .tol-j  iiif  ,  f»riS*Mi  232IM  3tm!q 


Consequently 


lim  /  0  ■■■d9  =  +  oo 

e^eo°  n  -  j 

dQ 

since  ^  is  a  decreasing  function  and  is  finite  and  negative  at 
n  =  j  .  Hence  the  integral  A ( 6 )  does  not  converge  when  the  upper  limit 
of  integration  is  0Q  and  g ( 0o )  =  0. 

The  scalar  invariant  X,  expressed  in  terms  of  the  physical 
components  of  the  strain  rate  tensor  and  the  stress  deviation  tensor 
is 


X  = 


2k 


!!r  +  ire  !V+(S  ,3)^ 

'r  9r  r  96  ls0  V  r  J 


and  must  be  positive  during  plastic  deformation.  Using  equations 
(4.1.5)  and  (4.1.7),  the  form  of  X  can  be  expressed  as 


r3n  -  4g 

Ln3  -  1 


■]  for  v 


r 


/  0  . 


Also  3n  -  43  is  negative  for  0  <  n  <  >  3  >  1.  Thus  in  the  region 

v3 

0  <  0  <  0Q  with  vr  <  0  and  as  given  by  equation  (4.3.1),  X  is  positive. 
However. in  the  region  0Q  <  0  <  a  ,  where  a  is  the  value  of  0  corres¬ 
ponding  to  t  =  m',  the  velocity  field  given  by  equation  (4.3.1)  and  the 
stress  field  result  in  negative  X.  If  a  positive  constant  B  is  assumed 
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in  the  region  0Q  <  0  £  a,  then  >  0  and  X  is  positive,  but  then  the 

positive  shearing  stress  at  0  =  a  does  not  oppose  the  flow  and  does 

work  on  the  material.  Consequently,  the  velocity  field  is  not  admissible 

in  the  region  9 Q  <  0  £  a  unless  B  is  zero.  Thus  the  region  0Q  £  ©  £  a 

does  not  deform  during  the  flow.  Similar  non-deforming  regions  were 

indicated  by  the  solution  for  the  plane  radial  flow  problem  considered 

in  CHAPTER  III.  Again  the  parameter  X  has  an  infinite  discontinuity 

at  0  =  0Q  but  this  is  not  inconsistent  with  the  equations  (4.1.1).  Of 

2 

interest  is  the  fact  that  this  discontinuity  in  X  occurs  when  sr  =  k  /y 
in  both  the  plane  converging  radial  flow  problem  and  the  axially  sym¬ 
metric  converging  radial  flow  problem  just  considered. 

The  function  g(0)  in  equation  (4.3.1)  was  evaluated  numerically 
for  the  cases  where  3  =  10  and  c  =  10  and  7  respectively  for  the  region 

0  <  0  <  6  and  where  the  maximum  s  -  =  m'k  is  attained  on  the  channel 
—  —  o  rt> 

walls.  Results  so  obtained  are  shown  graphically  in  Figure  8.  For 
practical  values  of  $  =  y/k,  the  angle  a  -  0Q  is  very  small.  The  re¬ 
lationship  between  a  -  0Q  and  a  is  shown  graphically  in  Figure  9  for 
8  =  10  and  100  respectively. 


78 


FIGURE  8  -  GRAPHS  OF  g(e)  FOR  a  =  15°06.7'  (8  =  10,  c  =  10) 


a  =  25°22 . 2 1  (@  =  10,  c  =  7) 


-  25°22_2 


79 


FIGURE  9  -  VARIATION  OF  a  -  0o  WITH  SEMI-ANGLE  a 
FOR  3  =  100  AND  10  RESPECTIVELY 


(RADIANS) 
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CHAPTER  V 


PLANE  RADIAL  FLOW  PROBLEMS  OF 
A  HYPO-ELASTIC  MATERIAL  OF  GRADE  TWO 

In  classical  elasticity,  deformations  of  a  material  are 
measured  from  a  'natural'  state  of  the  material,  taken  to  be  an  equi¬ 
librium  state  without  external  forces  and  without  stress.  ZAREMBA  [43] 
and  JAUMANN  [44]  suggested  dropping  the  idea  of  a  'natural'  state  and 
proposed,  instead  of  the  usual  stress-strain  relations,  relationships 
between  the  rate  of  stress  and  the  rate  of  strain.  Generalizing  these 
works,  TRUESDELL  [45]  proposed  the  theory  of  hypo-elasticity.  The  re¬ 
searches  of  JAUMANN  are  applicable  to  a  'material'  classified  in  current 
terminology  as  hypo-elastic  of  grade  zero.  TRUESDELL ' S  work  is  appli¬ 
cable  to  a  wider  class  of  materials  since  the  theory  of  hypo-elasticity 
reduces  to  the  classical  theory  of  elasticity  for  infinitesimal  defor¬ 
mations  of  the  material  measured  from  the  'natural'  state  configuration 
and  may  even  be  regarded  in  special  cases  as  an  enlargement  of  the  theory 
of  large  elastic  strains.  Moreover,  hypo-elasticity  theory  predicts 
in  special  cases  a  yield-like  phenomenon  called  hypo-elastic  yield  with¬ 
out  the  assumption  of  a  yield  condition. 

In  this  chapter,  the  comparison  as  exhibited  by  GREEN  [46]  and 
TRUESDELL  [47]  between  a  hypo-elastic  material  of  grade  two  and  an 
elastic-perfectly  plastic  MISES  material  is  given.  Also,  two  flow  problems 
are  considered^  namely,  the  steady  state,  quasi-static,  plane  radial  flow 
of  an  incompressible  hypo-elastic  material  of  grade  two  through  a  con¬ 
verging  infinite  channel  with  smooth  walls  and  a  similar  problem  with 
dynamic  effects  considered. 
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5.1  HYPO-ELASTICITY  AND  PLASTICITY 
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A  hypo-elastic  material  is  one  whose  constitutive  equations 
may  be  written  in  the  form 


(5.1.1) 


•  •  •  •  •  • 

where  H1Jmn  =  HJlmn  =  Hljnm  are  the  components  of  a  hypo-elastic  re¬ 
sponse  function  -  an  isotropic  tensor  function  of  the  stress  components 
a1J*.  TRUESDELL  [48],  NOLL  [23],  and  RIVLIN  and  ERICKSEN  [24]  have 
shown  that  when  the  response  function  is  a  polynomial  function,  the 
components  Hljmn  can  be  represented  as 


where  C-j , 
nomials  in 


Ljijmn  _  lm  ,r  in  r  in  r  lr  nv 
H  =  g  C C-j g°  +  C20  +  C3a  ar  ' 


•  0 


im/r  in  ,  r  in  r  jr  n% 
+  a  C CAg°  +  C^a  +  C6a°  ar  ) 


IS  m/r  in  r  j  r  i  ,  r  j  r  1 1  \ 
+  a  as  (Cig  +  C8a  +  Cga  ar  ) 


Jn  X  r  11 


+  ^C10(giJgmn  +  gmVn) 


+  lc11(giVm+  gmnaji  +  gJ'Vn  +  gjiamn) 


Ic^g'VV  +  ♦  gjVpap"  ♦  gJ1ompcp"), 


C  are  called  the  response  coefficients  and  are  poly- 
the  stress  invariants  J-|,  J2»  and  J3  de^ned  as 


aJMnoqmoo  »=,•••  5  k  sriJ  *o  no  i*nui  i®**  '  '  Ml  m*  K!,oq* 


26  bft^i52'3'  S'* 


„  ft*.  m  m,„  =  n  »H 


Moreover,  a  hypo-elastic  material  is  said  to  be  of  grade  n  if  Hn^mn 
are  polynomials  of  degree  n  in  the  stress  components  a1J.  For  a  hypo- 
elastic  material  of  grade  n,  n  >_  1 ,  the  constitutive  equations  (5.1.1) 
are  form  invariant  under  a  change  of  objective  stress  rate.  Not  being 
so  for  n  =  0,  hypo-elastic  materials  of  grade  zero  are  considered  in 
TRUESDELL  and  TOUPIN  [10]  as  being  physically  irrelevant. 

Since 


«_  Jj  ,  Paij 

st  a  '  vt 


(o 


ingmj 


+  anjgim)d 


mn 


9 


where  is  the  JAUMANN  derivative  of  stress  given  by  equation  (2.2.9), 
the  constitutive  equations  (5.1.1)  may  be  rewritten  as 


Vail  =  (hU™  +  aingmj  +  anjgim}  d 
P  t 


mn 


or,  on  using  the  mixed  components  o»,  as 


Vo j 

Vt  Jpj 


=  g  .H1 Pmnd  +  o1d?  +  oPdl 
9PJ  mn  p  J  J  P 


(5.1.2) 


Expanding  out  the  right  hand  side,  equation  (5.1.2)  becomes 


ad  (S  T.i  )  norJscpa  «jb  ?  bm>(  .p i  v  :  \  o  ocx3 
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+  +  Cg¥aj  +  Cyd^^a!?  +  Cg^cr^a!? 


+  Cn'i'a^aP  +  C  -,  n  d  1  +  C-ntd^aP  +  a^dP) 
9  p  J  10  j  llpj  pj 


(5.1.3) 


+  C-|  2  ( dJ.CTpCT  j  +  <^>  +  oJdP  +  cjdj. 
where  $  =  aid'?  and  ¥  =  gla^d.  are  joint  invariants  of  the  stress 

J  1  J  K 

tensor  and  the  rate  of  strain  tensor. 

For  the  following  particular  choice  of  the  response  coefficients 
C.  (i  =  1 . 12)  as  proposed  by  GREEN  [26]: 

C  =  .  iusd  (c1)2  c  =  C  --  -P--  q1 

L1  l-2v  27  'ai '  >  L2  ^4  9  ai  ’ 


4ya‘ 


C5  3 


»  C  i  o  ~  2y,  C-j  i  1 » 


=  0  otherwise, 

and  with  v  the  POISSON  ratio  and  a  a  material  constant,  one  obtains  from 
equation  (5.1.3)  the  constitutive  equations  for  a  class  of  compressible 
hypo-elastic  materials  of  grade  two,  namely. 


]_?±. 

2p  Vt 


v  Akxi 
T2v  dk6j  ' 


2a 


(s"fn)sl 
n  m  j 


(5.1.4) 


Contraction  of  equation  (5.1.4)  yields 


1  Vok  1+v  ,k 
2 vVT  "  l-2v  ak 


(5.1.5a) 
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or  simply 


D£=  . 

Dt 


K  d 


k 

k  • 


Thus  from  equations  (5.1.4)  and  (5.1.5a),  one  gets 


2y  Vt 


2a 


(sVls1 
'  n  m'  j 


(5.1.5b) 


which  together  with  equations  (5.1.5b)  are  the  constitutive  equations 
of  a  compressible  hypo-elastic  material  of  grade  two.  For  an  incom¬ 
pressible  hypo-elastic  material  of  grade  two,  the  constitutive  equations 
are 


and 


Vs 


2y  Vt 


=  d1  - 
J 


2a‘ 


(s7)sl 
'  n  nr  j 


(5. 1 ,6a,b) 


which  upon  identifying  ^5-  with  -U  are  seen  to  be  equivalent  to  those 

J  2k 

for  the  incompressible  elastic-perfectly  plastic  MISES  solid  in  a  plastic 

state.  For  the  hypo-elastic  solid  the  second  term  on  the  right  hand 

side  of  equation  (5.1.6a)  is  present  for  any  value  of  s^.j  that  occurs 

durina  deformation  whereas  for  the  elastic-perfectly  plastic  solid  it 

•  *  2  i  VsJi 

is  zero  unless  sis'?  =  2k  and  s.  0* 

By  considering  simple  shear  flow  of  the  hypo-elastic  solid 

with  constitutive  equations  (5.1.6a,b),  TRUESDELL  [47]  has  shown  that 


*ri  ibr  noqu  ri^Mw 


2'1UD30  fbftt  ..2  ^0  9U  T  6V  \nf>  ^0  '  1  *j  MG  3  r  <  i.  3 .  r  )  Wfi  M  0  9b  ^  2 
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as  the  shearing  strain  is  increased  the  corresponding  shearing  stress 
first  increases  to  a  maximum  value  called  the  hypo-elastic  yield  and 
then  decreases.  Also,  if  ^  >  1,  this  shearing  stress  approaches, 
asymptotically,  a  value  lower  than  the  hypo-elastic  yield  which,  along 
with  the  normal  stresses  required  to  maintain  the  shear  flow,  satis¬ 
fies  the  VON  MISES  yield  condition.  For  large  values  of  ,  which 
are  typical  for  metals,  the  hypo-elastic  yield  in  simple  shear  flow, 
the  asymptotic  value  and  the  shearing  stress  predicted  by  the  VON  MISES 
yield  condition  are  very  nearly  equal  and  the  normal  stresses  required 
to  maintain  shear  are  negligible.  The  solutions  to  the  radial  flow 
problems  considered  in  CHAPTER  III  and  CHAPTER  IV  are  not  valid  if 
<  1  and  the  reason  for  this  may  be  the  non-existence  of  a  VON  MISES 
type  yield  for  the  hypo-elastic  solid  if  ^  c  1.  This  is  not  of  great 
practical  importance  since  all  common  elastic-plastic  solids  have  large 
values  of  . 

5.2  HYPO-ELASTIC  FLOW  IN  AN  INFINITE  CONVERGING  CHANNEL 

The  steady  state,  quasi-static  plane  flow  of  an  incompressible 
hypo-elastic  material  of  grade  two  in  a  converging  infinite  channel  with 
smooth  walls  is  considered.  The  geometry  of  the  problem  is  the  same  as 
for  the  plane  flow  problem  of  an  elastic-perfectly  plastic  MISES  solid 
considered  in  CHAPTER  III. 

Since  plane  incompressible  flow  is  considered,  the  only  non¬ 
zero  physical  components  of  stress  deviation  and  velocity  are  s^,  Sq 
and  v  .  These  components  are  functions  of  r  alone  with  s^,  -  ~  Sg  since 

the  flow  is  axially  symmetric.  The  three  unknowns  sf,  vp  and  the  hydro¬ 
static  pressure-p  are  determined  by  three  governing  equations;  namely, 
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the  constitutive  equation 


ds^  dvM  0  2  dv  v  o 
v  _ L  =  our _ L  2o_  / _ r  rw  2-, 

^  dr  3  dr  "  r  sr  * 


r  dr 


(5.2.1) 


the  equation  of  incompressibility 


dv  v 
_J1  +  JL  =  o 
dr  r  u’ 


(5.2.2) 


and  the  non-trivial  equation  of  equilibrium 


r  Hr  ^sr"p^  +  2sr  =  °* 


(5.2.3) 


together  with  an  arbitrarily  assigned  initial  value  for  s .  This 
condition  on  s^  completely  determines  the  stress  field  to  within  an 
arbitrary  hydrostatic  pressure. 

Equation  (5.2.2)  requires  that 

vr  =  -  r  >  0  (5.2.4) 

with  the  constant  A  >  0  for  converging  flow.  Equation  (5.2.1)  then 
gives 


^r  =  2m 

dr  kr 


(5.2.5) 


where  kE^,a>0.  Integration  of  equation  (5.2.5)  yields 
2a 


.  wo  r 


■  -.■> 


ups  to  nort6,tg9tnI 
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and 


k  -1  r 

-  ^  tanh  —  =  log  r  +  log  C  if  |s  |  <  k, 

k  -1  Sr 

-  2jJ  coth  —  =  log  r  +  log  C  if  |sr|  >  k, 


(5.2.6a,b,c) 


sr  =  ±  k  if  |sr|  =  k. 

C  is  a  positive  constant  of  integration.  If  sr  =  0  at  r  =  rQ  >  0,  then 
from  equation  (5.2.6a), 

r 

sr  =  k  tanh  log  —) .  (5.2.7) 

s 

Figure  10  shows  the  graph  of  tt-  plotted  against  —  and  shows  that 

K  ro 

s  changes  algebraic  sign  at  r  =  r  and  that  as  y  -*  «>  for  fixed  k  the 

Sr  sr 

curve  of  approaches  pointwise  the  curve  —  =  1  in  the  interval 

0  <  —  <1  and  the  curve  -r~  =  -  1  in  the  interval  1  <  zr  <  °°*  If 
—  r  k  r 

o  o 

initially  sf  =  k^ ,  |k-j|  >  k  at  r  =  rQ,  then  from  equation  (5.2.6b)  it 

follows  that  sf  is  infinitely  discontinuous  at  r  =  rQ.  Such  behaviours 

in  s  are  considered  unrealistic  and  the  solutions  (5.2.7)  and  (5.2.6b) 
r 

are  considered  inadmissible. 

The  initial  condition  sf  =  k  at  r  =  rQ  >  0  gives  the  solution 

=  k  for  all  r  >  0,  (5.2.8) 


and  from  equation  (5.2.3)  it  follows  that 


I 
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p  2k  log  r  +  log  C-j ,  C-j  a  positive  constant. 
If  P  =  P0  r  =  rQ,  then 

p  =  2k  log  -p-  +  p 
o 

and  hence 


or  =  k  -  2k  log  y-  -  P0> 

0 


-k  -  2k  log  y-  -  p  , 
0 


and 


(5.2.9a ,b,c) 


az  =  -  2k  l09  r  -  p0- 

o 


The  stress  field  for  the  corresponding  problem  involving  an  elastic- 
perfectly  plastic  MISES  solid  is  obtained  from  equations  (3.2. la, b  and  3.1.3a) 
on  setting  ip  =  0.  It  is  identical  to  the  stress  field  defined  by 
equations  (5.2.9a,b,c)  if  in  these  equations  the  material  constant  k 
is  interpreted  as  the  asymptotic  value  of  the  shearing  stress  of  the 
hypo-elastic  material  of  grade  two  in  simple  shear. 

If  the  above  flow  problem  is  not  quasi-static  and  dynamic 
effects  are  taken  into  consideration,  then  the  equation  of  motion 


P)  + 


2s. 


=  pv 


dv, 

r  dr 


r 


(5.2.10) 


. 
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must  be  used  rather  than  the  equation  of  equilibrium  (5.2.3).  Since  the 
flow  is  axially  symmetric,  the  velocity  vf  is  a  function  of  r  only  and 
hence  from  the  equation  of  incompressibility  (5.2.2), 


(5.2.11) 


with  the  constant  A  >  0  for  converging  flow.  The  stress  deviation 
component  sr  is  then  as  given  by  equations  (5. 2.6a ,b ,c)  depending  upon 
the  specified  initial  value  of  sr  at  r  =  rQ>  0.  For  the  same  reasons 
as  discussed  previously,  the  only  admissible  value  taken  for  sr  is 

sr  =  +  k  (5.2.12) 

for  all  r  >  0.  Substitution  of  equations  (5.2.11  and  5.2.12)  into 
equation  (5.2.10)  and  integration  yields 

1  A  2 

p  =  2k  log  r  -  j  p  +  log  C-, ,  C-j  a  positive  constant. 

r 

If  p  =  pQ  at  r  =  rQ,  then 

p  =  2k  log  -p—  +  (—^2  "  +  P0* 

0  ro  r 

The  stress  components  are  then 

A2  1  1 

a  =  k  -  2k  log  f" +  2  - $  '  po  * 

r  0  r  ro 

afl  -  -  k  -  2k  log  f-  +  ^  -  -7]  -  P0> 

0  0  r  r 
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and 


a 


z 


P 


=  -  2k  log  —  + 
ro 


CHAPTER  VI 


CHARACTERISTIC  STUDY  FOR  PLANE  STRAIN  FLOW  OF  AN 

INCOMPRESSIBLE  ELASTIC-PERFECTLY  PLASTIC  MISES  SOLID 

The  governing  equations  for  the  plane  plastic  flow  of  a  rigid 
perfectly  plastic  MISES  or  TRESCA  solid  form  a  hyperbolic  system  and 
the  characteristics  of  the  stress  and  velocity  fields  coincide  [49]. 

The  characteristics,  called  a  -  lines  and  3  -  lines,  are  mutually 
orthogonal  and  are  also  called  slip  lines  since  their  directions  at 
every  point  coincide  with  those  of  the  maximum  shear  strain  rate.  By 
convention,  if  an  a  -  line  and  3  -  line  are  regarded  as  a  pair  of  right- 
handed  curvilinear  reference  axes,  the  line  of  action  of  the  algebraically 
greatest  principal  stress  lies  in  the  first  and  third  quadrant.  The 
state  of  stress  at  a  point  in  the  deforming  material  is  completely  speci¬ 
fied  by  the  mean  compressive  stress-p  and  the  angle  <J>  which  is  taken 
to  be  the  anti-clockwise  angular  rotation  of  the  a  -  line  from  the  x  - 
axis  of  a  CARTESIAN  coordinate  reference  system.  The  two  families  of 
characteristics  are  defined  by 

=  tan  $  for  a  -  lines, 


and 


dy. 

=-cot  d>  for  3  -  lines, 
dx 

Moreover,  for  the  stress  field,  the  compatibility  relations  are  the 
HENCKY  equations  [50], 


92 


nor  Sob  9nK  bM  ,  29**  teanr  rrv'tu^  bsbnfirt 

.  irf  -  jo  9d^  norJ&Jcn  ^Tugne  92rw>ioofo-r}n6  erl 3  9d  o3 


;  1 3  906  zno r36f9o  \3  f  \  fdfJeqmoa  ertd  ,bf9*r^  229^*2  9dd  10^  «‘19V09'ioM 


93 


dp  +  2k  d<J>  =  0  along  an  a  -  line, 


and 


dp  -  2k  d4>  =  0  along  a  3  -  line, 

with  k  the  appropriate  yield  stress  in  pure  shear.  For  the  velocity 
field,  the  compatibility  relations  are  the  GEI RINGER  equations  [51], 

du  -  v  d4>  =  0  along  an  a  -  line, 


and 


dv  +  u  d<f>  =0  along  a  3  -  line, 

where  u  and  v  are  velocity  components  in  the  positive  a-  and  3- 
di recti ons  respectively. 

In  this  chapter,  a  corresponding  study  is  made  of  the  governing 
equations  for  the  quasi-static,  steady  state  plastic  deformation  of  an 
incompressible  elastic-perfectly  plastic  MISES  solid  in  plane  strain. 

i  k  i 

It  is  shown  that  if  |-|  <1  this  system  of  equations  admits  four  distinct 

y 

families  of  real  characteristics.  Also  compatibility  relations  along 
these  characteristics  are  derived.  In  the  limiting  case  of  the  rigid- 
perfectly  plastic  MISES  solid  (y  ■+  °°) ,  the  GEI  RINGER  equations  are  re¬ 
covered  and  the  stress  and  velocity  characteristics  coincide.  In 
this  chapter, cylindrical  polar  coordinates  (r,6,z)  are  used.  In 


,  enr  r 
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APPENDIX  B  ,  the  study  is  undertaken  in  a  different  manner  and  rec¬ 
tangular  CARTESIAN  coordinates  are  used. 
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For  axially  symmetric  flow  of  a  rigid-perfectly  plastic  MISES 
solid,  PARSONS  [12]  has  shown  that  the  governing  equations  admit  no 
real  characteristics  except  possibly  a  curve  in  a  meridional  plane 
along  which  the  radial  velocity  is  zero.  No  corresponding  study  is 
made  in  this  thesis  of  the  system  of  governing  equations  for  an  in¬ 
compressible  el astic-perfectly  plastic  MISES  solid  deforming  in  the 
plastic  state  under  axially  symmetric  conditions.  In  contrast,  for  a 
rigid-perfectly  plastic  TRESCA  solid  for  which  the  HAAR-KARMAN  hypo¬ 
thesis  has  been  adopted,  it  is  known  that  the  governing  equations  for 
stress  and  velocity  are  hyperbolic  with  characteristics  which  coincide 
with  the  slip  lines  in  a  meridional  plane  [52].  As  yet  there  are  no 
published  researches  pertaining  to  the  plastic  flow  of  an  el astic- 
perfectly  plastic  TRESCA  solid  in  finite  strain. 

6.1  CHARACTERISTICS  OF  THE  STRESS  FIELD  AND  VELOCITY  FIELD  FOR  PLANE 

STRAIN  DEFORMATION  OF  AN  INCOMPRESSIBLE  ELASTIC-PERFECTLY  PLASTIC 

MISES  SOLID 

The  governing  equations,  in  curvilinear  coordinates,  for  the 
plastic  flow  of  an  incompressible  elastic-perfectly  plastic  MISES  solid 
are  given  in  CHAPTER  II,  SECTION  2.4  of  this  thesis.  In  APPENDIX  A, 
SECTION  A. 4  these  equations  are  expressed  in  cylindrical  polar  coordi¬ 
nates  (r ,0 ,z) .  The  number  of  governing  equations  is  greatly  reduced, 
however,  for  steady  state,  quasi -static  plane  flow  with  the  deformation 
independent  of  z  and  parallel  to  the  (r,0)  -plane.  Since  the  solid  is 
incompressible, each  incremental  distortion  in  a  state  of  plane  strain 


-09i  bn6  ^9nn6m  itmsb*  tb  6  nr  rt9>U-T9bnu  ar  ybuJa  art*  ,  - 

< 
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consists  only  of  a  pure  shear.  The  stress  component  az  normal  to  the 

plane  of  flow  is  thus  equal  to  the  mean  hydrostatic  pressure  -p  and 

the  stress  deviation  component  sz  vanishes  identically.  Also  the 

velocity  component  v  and  the  stress  components  a  and  a_  are  zero. 

^  rz  wz 

It  then  follows  from  equation  (A. 4. 3)  that 

sr  =  -  sQ  (6.1.1) 


and  from  equation  (A. 4. 2)  that 


sre  ■  /k2  -  sr2 


(6.1.2) 


With  appropriately  specified  boundary  conditions,  the  unknowns  sr>  vr, 
v0  and  p  are  determined  by  the  constitutive  equation 


3sr  *  Ve  3sr  +  ,sr  3vr  .  1  /J7TT  ,wsr  3v, 

vr  3r  +  r  30  +u^2  ^  3r  +  r  ^k  sr  ^2  ^  30 


+  JF~- 


2  ,lJSr 


r  'k2 


3v0  y  sr2  3v0  _  1  /2 


2  US, 


s-  1)  #  -  =  7/k  -  V  D  ve  t6-1 


ys  v 
h  r  r 

,  2  r 


obtained  from  the  PRANDTL-REUSS  equation  (A. 4. 4a)  with  appropriate 
substitutions,  the  equations  of  equilibrium 


!ii 

9r 


r  /k2  -  s  2 


9sr  9£  =  _  ffr 
96  3r  r 


(6.1.3b) 


.3a) 
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1  !!r  .  1  9£  =  2sr9 
r  30  "  r  ae  r 


and  the  equation  of  incompressibility, 


(6.1.3c) 


+  l^e 

3r  r  ae 


(6.1.3d) 


For  the  determination  of  the  characteristics  and  compatibility  re¬ 
lations  along  any  existing  characteristics,  the  above  system  of  partial 
differential  equations  are  supplemented  with  the  relationships 


9si 

8r 


as. 


dr  + 


ae 


de  =  ds. 


9 


+  §§■ de  =  dp’ 


dr  + 


de  =  dv 


r  * 


(6.1 .3e,f ,g,h) 


dr  dr  +  ae  de  dV 

The  non-zero  physical  components  of  the  stress  deviation  tensor  are 
parametrized  in  terms  of  the  angle  ip  between  the  line  of  action  of  the 
algebraically  greater  principal  stress  and  a  radius  by  setting 

sr  =  k  cos  2\p 


(6.1.4) 
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and  hence* by  equation  (6.1.2), 

srQ  =  k  sin  2\p. 

The  characteristics  of  the  stress  field  are  determined  from 
equations  (6.1 .3,b,c,e,f) .  Substituting  sr  given  by  equation  (6.1.4) 
in  these  equations: 

sin  ^  cos  If  +  If  =  +  "r  cos  2^’ 

cos  2ip  +  y  sin  2ty  ||  -  7  |^  =  -  7  sin  2ip,  (6.1 .5a,b,c,d) 

|idr  +  f|d9  =  <#, 


and 


If  dr  +  i de  ■ dP 

where  P  e  The  coefficient  determinant  A  of  equations  (6.1 .5a,b,c,d) 

2k  5 

is 

sin  2 ip  -  1  cos  2<p 

„  _  cos  2^  7  sin  z<l> 

As  =  r 

dr  de 

0  0 


1  0 

0 

0  0 

dr  de 


or;  sup®  nr 


,  a  zoo  -  +  =  S|  -  :  S  *oo  ^  *  f[  *  "»» 


and  the  determinant  equation  A$  =  0  yields,  after  simplification,  the 
two  ordinary  differential  equations 
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and 


(6.1 ,6a, b) 


=  tan  (\p  +  J) . 


Equations  (6.1. 6a ,b)  define  two  families  of  orthogonal  characteristics 
curves  of  the  stress  field.  The  system  of  equations  (6.1. 5a ,b ,c ,d)  are 
field  equations  governing  the  stress  field  and  apply  to  any  incompressible 


material  deforming  in  plane  strain  and  for  which  the  relationship  (6.1.2) 


is  valid.  Since  this  relationship  holds  for  a  rigid-perfectly  plastic 
MISES  or  TRESCA  solid  in  plane  strain,  curves  defined  by  equations 
(6.1 . 6a , b )  are  also  characteristics  of  the  stress  field  for  these  solids. 
Figure  11  illustrates  the  geometrical  relationship  of  the  characteristic 
curves  with  respect  to  the  cylindrical  polar  coordinate  system  and  the 
CARTESIAN  coordinate  system.  It  is  seen  that 


*  +  f  =  e  +  <p 


(6.1.7) 


and  hence  the  family  of  characteristics  defined  by  equation  (6.1.6a) 
are  a  -  lines  and  the  family  defined  by  equation  (6.1.6b)  are  8  -  lines. 

To  obtain  the  compatibility  relations  along  the  a  -  lines  and 
8  -  lines  of  the  stress  field,  the  terms  in  any  one  column  of  the  coef¬ 
ficient  determinant  are  replaced  by  the  corresponding  terms  appearing 


4 
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FIGURE  11  -  DIAGRAM  SHOWING  RELATIONSHIPS  BETWEEN  THE 
CHARACTERISTICS  OF  THE  STRESS  AND  VELOCITY  FIELDS 
AND  ALSO  CONVENTION  USED 


) 


o 


on  the  right  hand  side  of  equations  (6. 1 .5a ,b,c,d) .  Setting  the  re¬ 
sulting  determinant  equal  to  zero  and  simplifying  yields 

dP  -  (sin  2i|<  +  cos  2^  ^-)(de  +  dtp)  =  0.  (6.1.8) 

The  compatibility  relationship  along  an  a  -  line  is  then  obtained  by 
1  d  r 

replacing  "pr  ^q"  gi ven  by  equation  (6.1.6a)  in  equation  (6.1.8)  to  give 

dP  +  d(6  +  ijj)  =  0 


or 


dp  +  2k  dcj)  =  0  ,  (6.1.8a) 

if  equation  (6.1.7)  is  used.  This  is  the  familar  HENCKY  relation  along 

1  dr 

the  a  -  line.  Similarly,  replacing  -  ^  given  by  equation  (6.1.6b)  in 
equation  (6.1.8)  yields 


dp  -  2k  deb  -  0  ,  (6.1 ,8b) 

the  HENCKY  relation  along  a  8  -  line. 

The  characteristics  of  the  velocity  field  are  determined 

from  equations  (6.1 .3a,d,g,h) .  Using  the  expression  for  sr  given  by 

equation  (6.1.4),  equation  (6.1.3a)  can  be  written  as 

0  3v  i  k  3v 

(cos2  2ip  -  2 )  jp1  +  7  sin  2^(cos  — 

+  sin  2ip  (cos  2<p  +  -)  -  7  cos  2i|i  39“ 


■ 
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r  sin  2^  (cos  2i|>  +  -jj)  Vg  +  Jr  cos^  2 ip  v 


+  2^sin  2i|>  v  +  2  -  s  i  n  2\b  — 

y  r  9r  r  y  ^  96  r  ' 


(6.1.9) 


The  coefficient  determinant  A  of  equations  (6.1.9)  and  (6.1.3d,g,h)  is 


\  E 


cos  2\p-2 
1 

dr 

0 


•Jr  sin2i^(cos2^--) 

y 

0 


de 

o 


sin2iJ;(cos2^+^-) 

0 


0 

dr 


-  Jr  cos ^2ip 

1 

r 

0 


de 


which  set  equal  to  zero  yields 


r2  (cos  +  T?)  "  r  2  sin  2\p 


y 


dr 


(cos  2\|>  - 


=  0  .  (6.1.10) 


d0 

This  equation  is  a  quadratic  equation  in  r  ^r  and  has  discriminant 
J\ -  (—)^  .  Thus  the  system  of  partial  differential  equations  (6.1.9) 

l  k  i  i  k  i 

and  (6.1.3d,g,h)  is  hyperbolic  if  |— |  <  1,  parabolic  if  |-|  =  1  and 

y  y 

i  k  i  k 

elliptic  if  |jj|  >1.  Considering  only  the  case  where  0  c  -  <  1,  which 
is  the  most  important  from  the  practical  standpoint,  and  on  defining 


-  =  sin  2y, 
P 


equation  (6.1.10)  yields  the  two  differential  equations: 


r  ^  =  tan  (ip  +  y  -  f) 


(6.1.11a) 


,r$  n(2 
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and 


r  ^  =  tan  -  y  +  J).  (6.1.11b) 

Equations  (6.1.11a,b)  determine  two  non-orthogonal  families  of  charac¬ 
teristics  for  the  velocity  field.  Figure  11  shows  their  geometrical 
relationship  with  respect  to  the  a  -  and  3  -  line  of  the  stress  field. 

The  family  of  characteristics  defined  by  equations  (6.1.11a)  are  called 
a-j-  lines  and  the  family  defined  by  equation  (6.1.11b)  as  3-|  -  lines. 

For  the  limiting  case  of  a  rigid-perfectly  plastic  MISES  solid  (y  ■*  0), 
the  characteristics  of  the  velocity  field  are  seen  to  coincide  with 
the  characteristics  of  the  stress  field. 

The  compatibility  relationships  along  the  -  and  3-j  -  lines 
are  found  by  firstly  replacing  the  terms  in  any  column  of  the  coefficient 
determinant  with  the  corresponding  terms  on  the  right  hand  side  of 
equations  (6.1.9)  and  (6.1 .3d,g,h) .  Setting  the  resulting  determinant 
equal  to  zero  gives 

2  sin  2y  l>r  If  +  F"  !f]  +  (cos  ^  +  sin  2y) 

|\!i  .  ds.  .  4§.]  .  1  (cos  2ip  -  sin  2y)  ^  =  0  (6.1.12) 

Lr  ds  dr  r  drJ  r  '  v  de 

where  s  is  the  arclength  along  a  characteristic.  On  defining  £  -  ^  'Is 
equation  (6.1.11a)  is  written 

r  =  tan  (y  -  C)- 


bs r r eo  9*16  (firt.r.d)  anotteups  bonrtsb  aorJaHsdDfiifiJ  o  to  orT 

9bra  bnsri  Jrlgii  si  *  no  zwsi  pnrbnoq29moo  9di  f  w  v,A  Jnwr  n9Jsb 


:  pn  ii  jb  iQ 


.aiJef  i o ) D6Ti]>flo  e  £iols  d;  .nsfoib  si-  ar  a  9i9riw 


Hence  along  an  a-j  -  line, 
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and 


de 

ds 


=  -  -  sin 


(?  -  y) 


(6.1.13a) 


=  cos  (c  -  y) •  (6.1.13b) 

Equation  (6.1.12)  then  becomes 

sin  2y  [vr  |£  +  -p-  ||]  +  1  sin  (t  +  y)[v0  cos  (?  -  y) 

(6.1.14) 

dvr  dvfl 

+  vr  sin  (?  -  y)  1  +  cos  (?  +  y)^-  -  sin  (c  +  y)^  =  0, 

if  the  equations  (6. 1.1 3a, b)  and  the  trigometric  relations 

sin  2y  +  cos  2\p  =  2  sin  (c  +  y)  cos  (c  -  y) , 

sin  2y  -  cos  2ip  =  -  2  sin  (c  +  y)  sin  (c  -  y) , 

are  used.  The  velocity  is  now  resolved  into  the  two  components  u  and 
v  where  u  is  the  component  along  the  a-j  -  line  and  v  the  component  along 
the  3-j  -  line.  Relationships  between  these  components  u  and  v  and  the 
components  vr  and  vQ  are 

v  =  u  cos  (c  -  y)  +  v  sin  (5  +  y)  (6.1.15a) 


'  c  rtf  o  rr  >  ,*  .  •  <  vvcn  *  yjr)o'*-’V  ••  ■T 

.'Jbccqi  oo  3.  n  iris  ?c  'terottSM.. 
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and 


v0  =  -  u  sin  (c  -  y)  +  v  cos  (c  +  y)  .  (6.1.15b) 

With  the  use  of  these  relationships,  it  follows  that 

7  sin  (c  +  y)[vq  cos  (c  -  y)  +  vr  sin  (c  -  y)] 

=  ~  sin  (5  +  y)  cos  2y  (6.1.16) 


and 


cos  (c  +  y) 


dv^,  dvQ  , 

-  sin  (c  +  y)  -  cos  2y 


ds 


ds 


+  (u  sin  2y  +  v) 


(6.1.17) 


Now  an  expression  for  along  an  a-j  -  line  is  presently  derived 

involving  the  variable  p  =  0  ^*11S  ^eads  t0  a  more  concise  form 

for  the  compatibility  relation  along  an  a,  -  line.  Since  ?  =  tt/4  -  ip,  then 


=  _ 

9r  ~  9r  * 

(6. 1 . 18a ,b) 

-  & 

96  "  "  90  * 


and  the  equations  (6.1.5a,b)  become 


cos  2c  If- 1  sin  25  If  =  If  -  1  sin  2c. 
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and 


sin  2?  fir  +  l  cos  2s  ||  =  -  1 1|  +  1  cos  2? 


3 r  r  *  80 


Solving  for  ff  and  §§■  gives 


30 


*£.  =  cos  2C  2E. 
3r  os  ^  3r 


_  r\  3P 

r  S1n  2?  80 


and 


(6.1.19a,b) 


=  1  -  (r  sin  2?  +  cos  2?  §£■)  . 

Hence  along  an  a-j  -  line, 

3s  =  cos  U  ■+  y)  -  7  sin  (c  +  y)  |f  -  1  sin  (c  -  y) . 

Equation  (6.1.17)  then  becomes 

dv  dvfi  , 

cos  (c  +  y)  -  sin  (c  +  y)  =  cos  2y  ^  -  (u  sin  2c  +  v) 

(6.1.20) 

1  sin  (c  -  y)  +  (u  sin  2y  +  v)[cos(c  +  y)  §f  -  f  si"  (S  +  y)  |^]. 

Also  from  equations  (6.1.15a,b)  and  (6.1.1 9a ,b) , 

sin2y[vr  +  ^r-ff]  =sin2Y[u{-cos(c  +  y)  ff  +  sin(?  +  y)  f  §£} 

+  v{sin(c  -  Y )  |f  +  cos  (?  '  y)  f  |f> 
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+  “  (u  s i n ( c  -  y)  -  v  cos(c  -  y) } ] .  (6.1.21) 

Substitution  of  equations  (6.1.20)  and  (6.1.21)  into  equation  (6.1.14), 
after  simplification,  results  in 

cos  2y  ^cis  +  v  ^fr  cos  ^  "  y)  “  ft  r  sin  ^  -  y)]  =  0 


or,  with  use  of  equations  (6.1.1 3a ,b) , 


cos  2y  [ 


du 

ds 


+  v( 


3P  dr  3P  d6 n  _  n 

3r  ds  30  ds'-1 


(6.1.22) 


Thus  the  compatabi 1 ity  relation  along  the  -  line  of  the  velocity 
field  is 


du  +  v_  =  0 

ds  2k  ds  u 


or 


du  +  dp  =  0  .  (6.1.23) 

As  observed  earlier  for  the  limiting  case  for  the  rigid-plastic  MISES 
solid  (y  -*  0),  the  characteristics  of  the  velocity  field  coincide  with 
the  characteristics  of  the  stress  field.  Along  the  ot  -  lines,  the  HENCKY 
relation  (6.1.18a)  holds  and  hence,  in  the  limiting  case  y  +  0,  the 
relation  (6.1.23)  becomes 


du  -  vd<|>  =  0, 
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the  GE I  RINGER  equation  for  an  a  -  line  of  the  velocity 
forming  rigid-perfectly  plastic  MISES  solid. 

The  compatibility  relation  along  a  3-|  -  line 
in  an  exactly  analogous  manner  to  that  given  above  for 
It  follows  from  equation  (6.1.11b),  which  must  be  used 
equation  (6.1.11a),  that 


field  of  a  de¬ 
may  be  obtained 
the  a-j  -  line, 
rather  than 


H?  =  F  cos  <?  +  y)> 
=  sin  (c  +  y)  , 


(6.1 .24) 


with  s  the  arclength  along  a  8-j  -  characteristic.  Hence  from  equations 
(6.1 . 5a ,b ) ,  it  follows  that 


^  |  costs  +  y)  -  Sin(c  -  y)  f~  -  J  costs  -  y)  ||  .  (6.1.25) 


Substitution  of  equations  (6.1 .15a ,b) ,  (6.1.24)  and  (6.1.25)  into  equation 
(6.1.14)  yields,  after  simplification,  the  relation 


cos 


u 


dP )  . 

ds ' 


0 


9 


or 


along  a  3-|  -  line. 


dv  +  Ik"  dp  =  0. 


(6.1.26) 


For  the  limiting  case  y  -*•  °°»  equation  (6.1.26)  becomes 


(es.r.3)  .  |§  <r  -  5)a°>  :  -  3f  (r  -  ?)"»*  -  (r  +  ;)*03  ::  ■  it 
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dv  +  udc}>  =  0  , 

the  GEIRINGER  equation  for  the  8-j  -  line  of  the  velocity  field  of  a 
rigid  perfectly  plastic  MISES  solid. 

Compatibility  relations  (6.1.23)  and  (6.1.26)  are  observed 
to  have  exactly  the  same  form  as  the  GEIRINGER  equations  if  the  latter 
equations  are  written  in  terms  of  u,  v  and  p  rather  than  u,  v  and  0 
as  is  customary  in  plasticity  literature, 

6.2  DETERMINATION  OF  THE  CHARACTERISTICS  OF  THE  STRESS  AND  VELOCITY 

FIELDS  FOR  THE  CONVERGING  PLANE  FLOW  PROBLEM 

In  this  section,  the  characteristics  of  the  stress  and  ve¬ 
locity  fields  are  determined  for  the  plane  strain  flow  problem  con¬ 
sidered  in  CHAPTER  III  and  which  involved  a  converging  infinite  channel 
with  perfectly  rough  sides.  Since  the  stress  and  velocity  fields  have 
already  been  determined,  no  use  is  made  of  the  method  of  characteristics 
in  determining  these  fields. 

In  SECTION  6.1,  the  differential  equations  for  the  character¬ 
istics  of  the  stress  field  for  the  plastic  flow  of  an  elastic-perfectly 
plastic  MISES  solid  in  plane  strain  are  given  by  equations  (6.1 . 6a ,b ) . 

For  any  specific  problem,  the  differential  equations  can  be  solved  if 
either  ^  is  a  known  function  of  r  and  e  or  ^  is  a  known  function  of  \p. 

ds 

For  the  problem  considered  in  CHAPTER  III,  is  given  by  equation 
(3.2.3b).  Hence  equations  (6.1.6a ,b)  may  be  written  as 

—  =  -  (LjLliO-JljL)  for  an  a  -  line,  (6.2.1a) 

r  c  -  cos  2\p' 


-f  >0  n9i  io'tq  wott  nr  -  e  i^lc;  r  •  w  *h  *t«-  *b  is  abTstt  tftoof 
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and 


dr  _  1  -  sin  2\j>  , 
r  c  -  cos  2\p  ^ 


for  a  3  -  line 


(6.2.1b) 


where 


c  >  1.1922. 


The  constant  c  appearing  in  equations  (6.2. la, b)  is,  of  course,  a 
function  of  the  semi-angle  a  of  the  converging  infinite  channel  as 
given  by  equation  (3.2.6).  With  a  change  of  variable  from  ip  to  z  = 
tan  ip ,  equations  (6.2. la, b)  become,  respecti vely , 


dr  r  z 
r  "  _2 


1 


z  +1 


c+1 


z2  + 


(£ii) 

[c+V 


z2  + 


y— ]  dz  for  an  a  -  li 


ne 


and 


dr  =  [ 

r 


z2+l 


1 

c+1 


1 


z2  + 


(£zl) 

c+1 


-y-]  dz  for  a  3  -  line, 


<5r> 


where  the  right  hand  side  of  each  equation  is  expressed  in  partial 
fraction  form.  Integration  of  these  two  differential  equations  give 
the  following  equations  for  the  characteristic  curves: 


r  =  b 


/ 


c+1 


c-cos2i|j 


exp  [- 


-  tan 


-1 


/-2 


c  -1 


tan  i^)]  for  an  a  -  line, 

(6.2.2a) 


and 


,  ?f  1:  <o  norJt  g£4fll 


no 


tamjO]  for  a  3  -  line, 


(6.2.2b) 


The  constant  b  is  the  constant  of  integration  and  the  angle  ifj  varies 


between  -  ^  to  ^  .  If  rQ  is  the  value  of  r  at  which  an  a  -  line  or 


o 


3  -  line  intersects  the  axis  6  =  ip  =  0,  then  b  =  r  /  — h-  . 

o  c+ 1 

The  equations  (6.2.2a,b)  are  also  the  equations  for  the  slip 
lines  of  stress  field.  NADAI  [40]  obtained  equations  (6.2.2a,b)  as  the 
equations  of  the  slip  lines  in  the  flow  problem  of  a  rigid  plastic  solid 
in  a  perfectly  rough  converging  channel  and  plotted  curves  of  the  slip 
lines  for  the  cases  where  a  =  24° 1 7 '  (c=2)  and  a  =  90°  (c=1.1922).  It 
follows  from  equations  (6.2.2a,b)  that  for  rQ  >  0  the  characteristics 
of  the  stress  field  terminate  on  the  channel  walls  at  finite  distances 
from  the  virtual  apex  of  the  converging  channel. 

The  characteristics  of  the  velocity  field  for  the  plastic  flow 
problem  of  an  elastic-perfectly  plastic  MISES  solid  in  plane  strain  and 
for  which  0  <  —  <  1  are  determined  from  equations  (6.1.11a,b).  For  the 

p* 

flow  problem  under  consideration,  these  differential  equations,  with  the 
use  of  equation  (3.2.3b),  may  be  expressed  as 


dr 


r 


/l-n2  +  sin  2\p  cos  2i jj  .. 

n  -  cos  2\p  c  -  cos  2\p  ^ 


for  an  -  line, 


and 


(6.2. 3a , b ) 


for  an  3i  -  line. 


q f  r 2  slid*  To'i  anorJ6up9  9uj  oaf  r  e  d«cr.S.  snord  e  art 


Ill 


Again  c  >  1.1922  and  is  related  to  a  by  equation  (3.2.5).  Also  n  =  —  . 

r* 

With  a  change  of  variable  from  to  z  =  tan  ,  equations  (6.2.3a,b) 
become 


dr 

r 


1  /Un 
c+1  / 1+n 


(1  - 


2z 


[- 


Av 


+  z2)(1  -  z2) 


1-n 

1+n 


z2)(|rf  +  z2)(1  +  z"> 


“]  dz  for  an  a-j  -  line, 


and 


dr 

r 


1  /Tn 
c+1  / 1+n 


(1  +  -T==+  Z2)(l  -  Z2) 


[- 


+  z‘> 


— ]  dz  for  a  $-|  -  line. 


Integration  of  these  equations  give  the  following  characteristic  curves 
of  the  velocity  field  valid  provided  cos  2\p  f  n: 


n 


(|tam|i  +  /]+"|)C  "  rc  /  1-n2  t  -1 

r  =  b1  - 3 - exp  —  tan 


(|^|  +  tan2^)2(c'n)cosij- 


c  -1 


and 


c+1 


tamp)]  for  a-j  -  lines, 


(6.2. 4a , b ) 


n 


r  = 


( |  tamp  -  /j^-l  )C_n  c 

b-  - - L-LL— - exp  [-  ft-t  /  -o- 


(cTT  +  tan2^)2^C"n^C0S^ 


c“n  c2-l 


tan 


-1 


(/^7y  tamp)] 


for  3-|  -  lines. 


The  constant  of  integration  b-j  is  related  to  r  ,  the  value  of  r  at 
which  the  a-j  -  or  3-|  -  line  intersects  the  axis  0  =  \p  =  0,  by 


b 


1 


1 _ 

r  (£d) c  (hR\  “  n  -|2  ( c- n ) 
Llc+1;  M+n;  J 


For  n  -*  0,  it  is  noted  that  equations  (6.2.4a,b),  in  the  limit,  be¬ 
come  equations  (6.2.2a ,b) .  This  follows  since 


_ ] _ _  /  c+1 

(cos*)(|=f+  tan2*)1'2  c-cos2^ 

and  lim  b-,  =  b.  This  result  is  expected,  since  for  plane  strain  flow 
n+0  1 

of  a  rigid  perfectly  plastic  MISES  solid,  the  characteristics  of  stress 
and  velocity  coincide. 

Figure  12  shows  the  characteristics  of  the  velocity  field  of 

the  flow  problem  considered  in  CHAPTER  III  for  an  elastic-perfectly 

plastic  MISES  with  —  =  0.10  and  a  perfectly  rough  channel  with  semi- 

y 

angle  a  =  24° 1 7 ' ( c=2 ) .  In  the  region  -0  <_  0  £  a,  corresponding  to 

-  i  cos’1  -  <  ifc  <  ic  >  the  ou  -  lines  are  tangent  to  the  radial  line 

2  y  —  —  4  I 

e  =  -  e  at  the  virtual  apex  r  =  0  and  terminate  on  the  radial  line 
o 

0  =  a.  In  the  region  -  a  <  0  <  -  0Q,  corresponding  to  -  j  <  ip  < 

-  1  cos’1  -  ,  the  an  -  lines  are  again  tangent  to  0  =  -  0  at  r  =  0 

2  u  I  0 

but  terminate  on  9  =  -  a.  The  S,  -  lines  are  reflections  of  a,  -  lines 
about  the  axis  0=0.  The  radial  lines  6  =  ±  9Q  are  limiting  lines  of 
the  velocity  characteristic  field.  Since  ±9Q  are  values  of  0  corres- 

k 

ponding  respectively  to  values  of  ip  such  that  as  cos2i|i  =  — ,  these 
limiting  lines  coincide  with  the  inner  boundaries  of  the  non-deforming 


sonfa  awoffcrt  aMT  .(d,6$.S.d)  anoUeups  amoo 


t'+O 


{^nut  *  '-••';)■  aoo) 


^0  bfoi^  VOfOO?  V  9,-ii  %0  .  O' J  :r1!)jDfc1ur!D  9f1J  aW'J><!  S  't'  |  i 

-rm92  rttrw  fannsrio  riguo'i  y,rj:>9Vi9q  fi  bn6  Or.'.  =  " rw  ri.  J'  jr.^biq 

oi  JnggnsJ  91 6  zvnU 

-29'nOO  8  ^0  29UffiV  9*t6  0±  90fir2  .biSTT  D  *f$2  i "IdfOG'iS  1-  ^IfDO’^V  srii 
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FIGURE  12  -  CHARACTERISTICS  OF  THE  VELOCITY  FIELD  FOR 

CONVERGING  FLOW  IN  A  PERFECTLY  ROUGH  CHANNEL 

a  =  24° 1 7 1  (c  =  2),  n  e  -  =  0.10 
'  y 
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regions  exhibited  in  CHAPTER  III.  If  the  frictional  shearing  stress 
on  the  channel  walls  is  equal  to  / 1  ■  nL  k,  n  e  -  ,  then  the  limiting 
lines  coincide  with  the  channel  walls;  if  less  then  they  fall  outside 
the  channel  walls  and  the  end  points  of  the  characteristic  curves  lie 
on  the  channel  walls  at  finite  distances  from  the  virtual  apex. 

6.3  CONCLUSIONS 

The  characteristic  study  in  this  chapter  has  involved  a 
frame  invariant  form  of  the  PRANDTL-REUSS  equations.  Provided  that 
the  incompressible  elastic-perfectly  plastic  MISES  solid  has  a  -  value 
less  than  unity,  the  governing  equations  for  plane  flow  are  hyperbolic 
and  hence  the  method  of  characteristics  can  be  used  for  solutions  of 
problems,  either  rotational  or  irrotational ,  involving  finite  defor¬ 
mations.  Extensive  search  of  the  literature  indicates  that  much  of  the 
researches  concerning  the  application  of  characteristics  to  elasto- 
plasticity  has  been  restricted  to  problems  involving  only  small  de¬ 
formations.  It  is  hoped  that  the  results  developed  here  can  be  used 
both  in  extending  our  knowledge  of  elasto-plasticity  and  in  solving 
many  technological  problems. 
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APPENDIX  A 


GOVERNING  EQUATIONS  IN  A  SPHERICAL  POLAR  AND  A  CYLINDRICAL 
POLAR  COORDINATE  SYSTEM  FOR  AN  ELASTIC-PERFECTLY  PLASTIC 
MISES  SOLID  IN  FINITE  STRAIN 

In  CHAPTER  II,  SECTION  2.4  of  this  thesis,  the  governing 
equations  for  an  elastic-perfectly  plastic  MISES  solid  in  finite  strain 
are  given  with  respect  to  a  general  curvilinear  coordinate  system.  In 
this  APPENDIX  A,  these  governing  equations  are  expressed  in  both  a 
spherical  polar  and  a  cylindrical  polar  coordinate  system  using  the 
physical  components  of  the  vector  and  tensor  quantities  involved.  The 
detailed  transformation  of  the  PRANDTL-REUSS  constitutive  equations  is 
included.  .However,  since  the  field  equations  expressed  in  these  two 
coordinate  systems  are  contained  in  various  texts,  notably  McCONNELL  [53] 
and  LANGLOIS  [54],  these  equations  are  only  entered  here  for  easy  reference. 

A. 1  CHRISTOFFEL  SYMBOLS  AND  PHYSICAL  COMPONENTS  FOR  A  SPHERICAL  POLAR 

COORDINATE  SYSTEM 

For  a  spherical  polar  coordinate  system  r,  6,  <j>  as  illustrated 
in  Figure  4,  the  line  element  is  given  by 

ds2  =  dr2  +  r2d02  +  r2sin20d<J>2  . 

The  covariant  and  contravariant  components  of  the  metric  tensor  are 
gn  -  1.  g22  =  r2>  g33  -  r2sin2e,  9ij  =  0(1«)  , 


and 
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11  i  22  1  33  1  ii  „ 

9  =  1 »  9  =  7»  9  =  -o - o-,  g  J  =  0  (l^j) ; 

r  r  sin  6 


and  the  non-zero  CHRISTOFFEL  symbols  of  the  second  kind  are 


^ 22  ~  ~  ^  ^ 


33 


2 

-  rsin  0, 


r12  =  r21 


1  r2 

r*  33 


=  -  sinecose, 


rl3  =  r31 


1  r3 
r’  l23 


3 

r32=cote. 


In  terms  of  the  corresponding  covariant  and  contravariant 
components,  the  physical  components  of  velocity  are 


vr  =  V1  =  y1 >  ve  =  T  *  ry2>  v*  =  rsfne  =  rsin9v''- 


and  the  physical  components  of  the  symmetric  stress  tensor  are 


11  .22  _  ,_2 _22  °33  2  .  2a  33 

a  -  a-,-.  =  a  »  ex*  -  ^  r  0  »  a4>  2  2  r  sin  0a  , 

r  11  0  r  v  r  sin^0 


r0  r 


12  =  ro12,  x 


r<j>  rsin0 


^  =  rsin0a^,  t 


04>  2  . 

Y  r  sin0 


23  2  .  Q  23 

•  =  r  sin0a 


The  physical  components  of  the  stress  deviation  tensor,  sr>  sQ,  s^, 

s  s  and  s  are  similarly  related  to  the  covariant  and  contra- 

r0  *  r4>  04> 

variant  components  of  the  stress  deviation  tensor.  From  the  relation¬ 


ships 


fe1  3  tr 


V  '  » rt  'J.  29  '--too  I  to  •'  1 1 

e^6  \t  o  v  to  >o»onc  ,  .;oo  iaotey  q  orlt  ,  3+ i9no<  n«  o 


,  venrzi  =  -gf;- -  *  ,  »i  =  *  Qv  ,fv  =  ,  -  V 


9"<6  -10209*  229-1*2  ontamm^  9ri*  to  aii  snoqmoo  T t>:o t ay,rtq  At  ona 


it  follows  that 


sr  "  ar  +  p>  s«  =  afi  +  P>  s*  =  +  P. 


e  “e 


(j)  <(> 


re 


Tr0  *  srcf>  "  Tr(j)*  s0<t>  "  T0<j)» 


with 


P  =  ■  I  (or  +  °e  +  V- 


The  physical  components  of  the  symmetric  strain  rate  tensor  are 


.  H  =i!!i  +  !r  j  .  V  .  cote 

ar  9r  *  a0  r  90  r  *  <j>  r  r  v0  rsin0  94> 


.  .  1,1  +  h  .i/_L 

ar0  2'r  90  9r  r  '*  r^  2vrsin0  94>  9r  r 


9vr  9v^  v 


,  -  if_L_  !!i  +  1  !!i  cote  , 
a0(j)  2Vsin0  9(j)  r  90  '  r  V  * 


and  the  physical  components  of  the  skew- symmetric  vorticity  tensor 


03 


r0 


3  V 

=  Ir^-aiT  "  fr  (rve^’  “r<t>  =  2^F  9lr^rv4i^  "  rsTne  3^’ 


“e*“  2?ik  Sr  -  !e  Vne)L 


n  :  "fL  ciis;  i  51  .iMiJ?  j  ie  9.1i  zrns  o  ,  too  f  •  r<i\(  C'T 


ev€  r 


_L_  +  ^ _  _L  =  h 

$C  •*  t  '  €  1 


V  kv6  v€  r  r  v  av6  iV€  r v  r 


v6 
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A. 2  GOVERNING  EQUATIONS  IN  SPHERICAL  POLAR  COORDINATES 

The  equations  of  motion  (2.4.2)  in  terms  of  the  physical  com¬ 
ponents  of  stress  and  velocity  are 


!!r  + 

1  3xr0  1 

3x  , 
r4> 

3r 

r  3r  rsine 

34> 

3xr0  + 

l  3oe  l 

3x0+ 

3r 

r  30  rsine 

34) 

3t,"0  + 

1  t04>  .  1 

3a , 

9 

3r 

r  30  rsine 

34> 

where 

> 

<ro 

> 

Q 

3  v, 

r  r 

Dt  =  3t  +  V 

r  3r 

Dv. 


Dv 


0 


V 


e 


r4>  c  L  0<f> 


Dv. 

_  f 

Dt  ^4> 


r  v0  3V  v*  3vr  1,  2,  2, 

r  38  rsine  3<J>  '  r^v6  +  v4>  )  • 


0ve 

Dt 


3t  r  3r'rv0'  r  30  rsin©  3<j> 


cote  2 

r  v4>  * 


Dv,  3v,  a 

_ i  =  _ i  +  jl  i_/rv  )  + 

Dt  3t  r  3r'rV 


Q  9  (cinflv  1  +  & _ 

rsine  36  IS1  V  rsine  3<}> 


p  is  the  density  and  Fr>  F0,  are  the  physical  components  of  the  body 
force  per  unit  mass.  For  quasi-static  problems  in  the  absence  of  body 
forces,  the  right  hand  sides  of  equations  (A.2.1a,b,c)  are  zero.  These 
equations  are  then  called  the  equations  of  equilibrium. 

The  equation  of  continuity  (2.4.3)  in  terms  of  the  density  and 
the  physical  components  of  velocity  is 


v6 


*  V 


.  6  1 

*  +  U  13 


■‘Sfen?H+  (e>vy"fj)  el'TT  T  *  i  t  7c 


ybod  9ftt  aJngnoqmoo  Uora^lq  aril  9T6  3  ,e^  bne  ^rjnsb  9ri>  2  q 
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ft  +  h  f7(r2pvr>  +  TiTne  |^(s1nepve)  +  f^)  =  0. 


If  the  solid  is  incompressible,  there  is  also  the  equation  of  incom¬ 
pressibility  (2.4.1b), 


!!n  +  !!n  +  l!!i+ cote  +^L_!!i=0 

3r  3r  r  30  r  0  rsine  3(f) 


The  identity  (2.4.5)  is 


sr  +  s@  +  s4>  *  °- 


(A.2.2) 


and  the  VON  MISES  yield  condition  (2.4.4)  is 


V2  ♦  s92  ^  ♦  2(sre2  ♦  s0/  +  .„«)  -  2k2  . 


which  with  use  of  equation  (A. 2. 2)  may  be  written  as 


or,  alternatively,  as 

(or-o0)2  +  (vV*  +  (Var)2  +  6(Tre2+Te<t>2+Tr4,  )  =  6k  • 

The  finite  PRANDTL-REUSS  constitutive  equations  (2.4.1a) 
in  terms  of  physical  components  of  stress  deviation  and  velocity  are 


-moanr  *o  (lorieups  srt*  oate  at  msitt  .-fdfaasiqrnowt  at  bi  oa  9(1*  '• 


‘ 
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(a) 


9s. 


9s-  vQ  9s^  v .  3s 

T  {W  -  2sre)  +  TiiHe  (sT  -  2sine  sre> 


r  .  ..  ,  '9 


9t  +  vr  3r  +  r  v30 


-  fer-  hr  (rva)]^-  [f-  (rvj  - 


0'  r  L3r  '  <Jr  ~  sin0  3<|>  r 


1  ^Vrj  sr4) 


(b) 


3v 

2y  ^3r“  '  Asr^  * 


(A. 2. 3a) 


9S0  9S0  VQ  9Sfl  V,  dsa 

3t  +  vr  aF”  +  F~  +  2sr0^  +  rsin0  ^3^  2cos0  sQ^) 


1  9Vfi  Vr 

2y  V  30”  +  F~  "  AS0^  * 


(A. 2. 3b) 


(c) 


9s  ,  9s  .  vq  3s,  v,  i  3s , 

_£  +  v  _£+_!_£+_£  fJ _ £  +  ?s  +  2cot0  s  1 

3t  vr  3r  r  30  r  ^sin@  3<j>  ^sr4>  ^cot0  sQ(p) 


1 


3v. 


9v, 


+  w  '  It  (rV]  -?*■ +  [3T  ■  fe  (sin6  V3  Fsine 


■  2tJ  (ZZ7 


!!i  +  !r  +  cote  v  .  Xs  ) 

rsin0  34>  r  r  e  AV’ 


1 


(A. 2. 3c) 


(d) 


3s  .  3s  fi  v„  3s  „ 

_J Ci+v  _J C2.+  _i  [_IH.  +  (s  -  SJ] 
3t  r  3r  r  L  30  v  r  0/J 


v  9  s 

+  - $-5-  ^  -  sin0  Sfl.  -  COS0S  ,  ) 

rsin0  9<j>  09 


(dt;S.A) 


Ws+V  V  -*-v  +  " 


26  ,  \ 


♦e2  e"^  '  «c 
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,3vr  3vfl 

+  r  W~  "  V  ~2r 


S  y,”  S  A  3v ,  ,  sfl. 

r  ®  +  (v  +  r  — -  — - - — )  -Q$- 

V(j>  r  9r  sine  d<p  '  2r 


(e) 


^<f.8  ‘  36  (sln9%)]  2rsine 


1  ,1  3vr  3vfl  Vfi 

'  2u  ^2  39~  +  W  '  *  Xsr6^ 


9s  ,  9s  ,  vn  9s 

— ri  +  v  +  _1  (_ii  .  s  ) 

9t  r  9r  r  v  96  be<j>' 


(A, 


v,  l 

+  ^  (7777  ^7^  +  S„  -  SQ  +  cote  Sr0) 


r  xsin6  9<f)  r  6 


+  &  <"♦> 


j_  fV,  Vfr  .  r!!i  i_  ,c,nflu  n  !re 
sine  9(j)  ^  2r  ^  96  (sin6v<j)^  2 r 


-  [!!l.  i_  (rv  )]  !§S 

Lae  9r  v  6;J  2r 


9v  9v ,  v 


=  2»  [I  (rsTne  W  +  3^  "  ^  '  *Sr<t>]’ 


(A, 


(f) 


9s Q,  9sA ,  vA  9sft , 

Q£  +  v  ii  +  _i  (_ ®£  +  s  ) 

at  r  9r  r  1  96  br<K 


+  (sfl-  sjcose  +  s  sine] 

rsine  L  9<J>  0  <p 


_9vv,  a  ,  x-.  srd>  r 9  /  \  1  9VTt  Sr6 

+  fer1  -  I?  <rvfl)]  2T  '  [3F  (rV  '  sine  3d>  ] 


2r 


+  tiT  -  le  (sin%)]  27?T7l 


2.3d) 


2.3e) 


ev!  <vf  f ,  r,  - 


<*,*  6,03  +  O2  '  -4  +*,-W*)  ^  + 


(  ven'ra)  -  -  -  . .  +  ni*.  1  [2  .-!-•-(  vi)  4] 


fcS  t(4venta)  ef  -  ^3  + 
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=  2y  [1  ( 


_ +  cote  ' 

rsine  3<J>  r  30  "  r  <jr 


1 


ASg^].  (A.2.3f) 


1  •  • 

In  these  equations,  the  joint  invariant  A  =  — d1.  s'?  ,  expressed  in 

2  k  j  1 

terms  of  the  physical  components  of  the  strain  rate  and  stress  deviation 
tensor,  is 


x  =  ^7  [sr  af~  +  T  (vr  +  +  ^  (vr  +  cot6ve  +  iTne 


3v. 


3v 


+  (Hr.  +  r  "  0 

r  v30 


3v _  3v 


.  v  )  +  _r±  (4 —  44  +  r  Hi  .  v  ) 

3r  V0 ^  r  'sm0  34>  3r  4>y 


3v^  3v 


+  (77^^+  ^  -  cotev.)]. 


r  vsin0  3<j)  30 


4>‘ 


Details  of  the  derivation  of  the  above  PRANDTL-REUSS  equations 
are  as  follows.  Rewritten  here,  equations  (2.4.1a),  the  PRANDTL-REUSS 
constitutive  equations  with  respect  to  a  general  curvilinear  coordinate 
system,  are 


Vs 


ij 


Vt 


=  2u(d1J  -  As1J) 


with 


(smjgki  +  si*gkjKm  , 


Vt  ~  Dt 


km 


=  ii^+  smn4ksim)vk 


at  '3xk  mk- 


mk' 


.  (s%ki  +  simg«Km. 


A2 

-mi  +  eve:t03  +  iv>  V  1+  ~e  S1  *  x 
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1  2  3 

With  x  =  r,  x  =  0,  x  =  <{>  and 

0)  i  =  1,  j  =  1, 


Dt  at  '  v9xk  '  Amk' 


1fr  +  (a7~+  2rmlsml)v1  +  (^7“+  2rm2sml)v2 


+  (*£  +  2r^sml)v3 
ax3  m3 


3S11  +  l41  V1  +  v2  (is!].  .  2x1s12)  +  v3  (l41  _  2xlsin2x2sl3) 


3t  '  3X1 


9x 


9xv 


9s  9s  v  9s  v  9s 

9t“  +  vr  9r“  +  7~  "  2srG^  +  rsine  ^9$~  “  2sin0  srV  * 


,ml  kl  .  1  m  kl  x  _  -2111  9c13  11 

(s  g  +  s  g  )ojkm  =  -  2s  g  u>12  -  2s  g  w 


13 


ira  _  9_  (rv  }]  .  fn  f9_  ( rv  )  .  J_  !!n 
r  l96  9r  v  0/J  r  L9r  v  V  sin0  9(j) 


], 


n  n  9v 

2u(d  -  Xs  )  =  -  Xsr), 


combine  to  yield  equation  (A. 2. 3a); 

(2)  1  -  2,  j  =  2, 

-  It $22  +  +  24s2V  +  +  2rm/>2 


* 


V  lf?  '7  t  ’",l  ‘l  +  -  — )  *  -  -  -  - 

S)lm  T~~ 


(£t2SxSnrarXS  -  -  M)  %  +  <Sf2fxS  -  M)  *v  +  f,  +  l'ji 


#S  ■  t  «S  •  1“  (s) 
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+  (^-  +  2r2,s2m)v3 

3x3  m3 


3  .22  ,  ,3s22  .  2  _22\..l  .  ,3s22  ,  2  _21\.,2 

IT"  7  )v  T2-  7  )v 


+  —  2  sinx^cosx^3)v0 

9x 


]_  !!i  +  r9_  (!e )  +  !!0]  v  +  r i_  (fe J  +  f!r0]  !e 

r2  at  +  Lar  \2>  r3  J  vr  Lae  \2>  r2  J  r 


+  rl_  _  2cos8  -|  cp 
Ld<p  '  2'  2  se<j)J  rsine  * 


/  m2  k2  2mnk2\  _  9  m2  22 

-  (s  g  +  s  g  )wkm  =  -  2s  g  a> 


2m 


(x1) 


2  /  12  .  .32  x 

1  2  (s  Wgi  s  ^23 


“T  {sre^ar  (rve)  ae  ^  +  sTne  ^a<j>  “  ae  ^sin0  V^* 
r 


2p  (d22  -  As22)  =  H  7 4r  +  7  -  *sfl> 


combine  to  form  equation  (A. 2. 3b); 
(3)  1  =  3,  j  =  3, 


Ds33  .  3s7+  (li7  +  2r31sm3)v1  +  (^2~+  2rm2sm3)v2 


Dt  3t 


3X 


3x 
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+  (^i+  2r33sm3)v3 


9x 


8s33  .  / 8s33  .  2  33\  1  ,  ,3s33 


2N_33 


at  +  )v  +  2  cot(x^)s 

3x  x  3x 


/3s33  ,2  _13  ,  <L6\  . 

+  (— T"  +  t  s  +  2  cot(x  )s  )v 

3x  ,  x1 


,2\_23\  3 


i  r9sd>  .  3si  +  ^e!fi 


r— i  + 


r2sin2e  3t  'r  3r  '  r  30 


+  f4,  (rnb-  a/+  2sr*  +  2  cote  seo>)]’ 


/  m3  k3  3m  k3x  _  0  m3  33 

-  (s  g  +  s  g  )ajkm  =  -  2s  g  oj 


3m 


9J3  33  q  23  33 
-  -  2s  g  go^i  -  2s  g  w 


32 


2  rSr4>  r  1  9vr  13/ 

r2sin20  1  2  Lrsine  94)  "  r  3r  [r\n 


+  2rsin0  ^  "  36 


2y  (d33  -  As33)  =  2y  (-5-^ 


1 


2  2  ?  ? 
r  sin  0  r  sin  0 


•) 


2v  (_] —  +  j:  +  cote  , 

2  .  2„  vrsin0  3<f>  r  r  v0  AV 
r  sin  0 


combine  to  form  equation  (A. 2. 3c); 


4>  e1  i  i  ,$*=  gSnt  ,y 


(4)  i  =  1,  j  =  2, 
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Ds12  3s12  "-12 


12 


Dt  9t 


+  r^lslm)vl  +  (5“+  r^m2  +  rm2$lm)v 


+  (2^+r^m2  +  4sl>3 


9s 


12  ,12 


12 


9t 


+  (as-  + 1  s12)vl  +  (as^+  (_xlj  s22  +  1  S11)V2 

3x  x1  9x^  X1 


+  -  x1sin2x2)s32  +  (-sinx2cosx2s13)v3 

9x 


i_  (ili)  +  [9_  till 9)  +  !l0]  v  +  ri_  (!li)  .  fi  +  V,  !e 

9t  1  r  '  L9r  K  r  ;  2  J  vr  L99  1  r  '  r  r  J  r 

r 


+  %  (-F)  -  ^  -  cos9  rifne  • 


9  /  reN  sine 


/  m2  kl  ,  lm  k2\ 

-  (s  g  +  s  g  )o,km 


22 _kl  ,  _12_k2 


(s1 2g1 1  +  s11gk2)a),1  -  (sfcV‘  +  s'V^)^ 


,  32  kl  ^  13  k2s 

-  (s  g  +  s  g  )a)k3 


11  22  22  11  32  11  13  22 

sV  w21  -  s^gMa)12  -  s  g  o)13  -  s  g  o)23 


9v 


-  ^2  {(Sr'Se)  [3T“  9f  (rve)] 
2r 


/!  o  xnra-)  +  “  ’(  V  r?  x  -  -  r;  L  + 


S<:.,  LS,  >J. 


fT.i  ■  r>  $pi  % 
92+  92)“ 
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a  9v 

+  s0(f)  ^9r  ^ rv4)^  "  sine  9$“-' 


3vfl  a 

'  snJ,  sine  ^‘30  (sinev^)], 


2„  /J2  , .12,  .  2^  rl  ,1  3vr  .  3v0  \  ,c  -1 

2y  (d  -  As  )  -  p-  .  xsr0] 


combine  to  form  equation  (A. 2. 3d); 
(5)  i  =  1,  j  =  3, 


Ds13  _  9s13  .  /9s13  .  r3  _lm\  1  /9s13  ,  rl  „m3  3  Jnu  2 

“Dt - 3F-+(y+Imls  )v  +(^2-+rm2s  +  fm2s  )v 


+  <2? +  rI/3  +  4sl>3 


3s13  +  (lij!  +  ^  s13)yl  +  (3ili  .  xls23  +  cotx2s13]y2 


9t  \  1  1 

9x  x 


9x 


13 


+  (^1  .  x1sin2x2s33  +  Jy  s11  +  cotx2s12)v' 
Sx'3  x 


1  r!!ra  +  v  !!rt  +  la  (!iri  .  s  ) 

rsine  L  3t  vr  3r  r  1  30  s0*' 


+  Sr  -  s0  +  COt9  Sr9^> 


-  (sm3gkl  +  slmgk3)u>km  =  -  s1 1g33u>3i  -  s^g11^  -  s12g33o>32  -  s^g11^ 


Sv[£raSxioo  +  £Vx  -  ♦  rv(£fa  /  ♦  -;f 


.  1  fVsr  r  3_  ,  >  1 _ 3vrn 

2rsin0  1  r  L3r  1  rv4) J  '  sine  3<t>  J 


Sr9  r3ve  3  ,  .  „  n 

+  “T  [3T  '  30  <slneV] 


!ei  [!!r  _  3_  (  )L 
r  l30  3r  '  VJ’ 


o  /j!3  ^  _ 1 3  x  2 u  rl  ,  1 

2u  (d  -  Xs  )  =  —f-E  [7  ( — 


9v  9v,  v 


i  .  _!)  _ 


rsin9  L2  vrsin6  3<J)  9r  r 


Xs 


combine  to  form  equation  (A.2.3e);  and  for 
(6)  i  =  2,  j  =  3, 


Ds23  _  3s23  ,  ,3s23  ,  „2  .m3  x  „3  .2rn„l 
“Dt"  '  3t  mls  mls  /v 


23 


+  ( 05  +  r2  sm3  +  r3  s2m)v2 

+  1  9  +  m2s  m2s  ;v 


9x 


+  (JS__  +  ^  JIIO  +  ^  s"»)v 
+  V  “5  1 mV  ;v 


23 


2  m3  ,  J  2m\  3 


9x 


3  '  1m3J  ’  "m3' 


(4  +  1tsZ3  +  jtsZ3|v1 

9x  x  x 


+  +  iy  S13  +  COtX2S23)v2 

9x^  X 


+  (ls!i  -  sinx2cosx2s33  +  ^-  s21  +  cotx2s22)v3 
3X3  * 


„v6 


.€  «  t  c$  *  t 
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r  sine 


]_  [IfMt  v  !!§i+  (!Xs  )^. 

-2--*  ~  9t  r  9r  [  96  sr<p'  r 


3s, 


+  ("a*  +  ^ s a  "  SJ  cose  +  s  Q  sine)  — 4— r-1  , 

d(p  e  <p  re  '  rsineJ  * 


•m2~k3  ^  _3m_k2\  _  31  22..  12  33 

km  "  "  s  9  21  “  s  9  W31 


/  me  ko  jm  k^n 

■  (s  g  +  s  g  )oo 


22  33  33  22 

s  9  w32  ~  ^  9  ^23 


2r3sin0  '"3e  "  3r  "  Sr6^3r  ^rv4>^  '  sine  3^ 


and 


se's<i)  r3ve  3  ,  .  „  n 

+  HTie  [3T  '  36  (sin6V]’ 


2y  (d23  -  AS23)  =  —  [1  (-4 


9vQ  i  3v,  .n 
_ . 0  ,  i cote  v 

r2sine  L2  Vrsine  94>  r  96  "  r  V 


Xse^ 


combine  to  form  equation  (A.2.3f). 

A. 3  CHRISTOFFEL  SYMBOLS  AND  PHYSICAL  COMPONENTS  FOR  A  CYLINDRICAL 

POLAR  COORDINATE  SYSTEM 

For  a  cylindrical  polar  coordinate  system  r,  e,z  as  illus¬ 
trated  in  Figure  2,  the  line  element  is  given  by 

ds2  =  dr2  +  r2de2  +  dz2. 


it  ss 


v 


The  covariant  and  contravariant  components  of  the  metric  tensor  are 


9 1 1  9£2  r  *  933  ^  >  9..  ~  9 


and 


911  =  1.  g“  -  4,  -  1,  g1J  -  0  (i/j); 


.22  1  33 


1J 


and  the  non-zero  CHRISTOFFEL  symbols  of  the  second  kind  are 


r1  =  -  r2  -  — 

22  r,  i12  i21  -  r  • 


In  terms  of  the  corresponding  covariant  and  contravariant 
components,  the  physical  components  of  velocity  are 

1  v2  2  3 

vr  =  v,  =  v  ,  v0  =  -  =  rv  ,  vz  =  v3  =  v  , 

and  the  physical  components  of  the  symmetric  stress  tensor  are 


11 


°r  "  Q11  =  a  *  a0  “Z 


22  «  . 2  22 


33 


r  0  »  a3  a33  "  a  » 


12  12 


13 


r6  r 


-  ro  ,  cr -j  3  0  ,  t 


0z  r 


23  23 

=  ra  , 


with  similar  relationships  holding  between  the  covariant,  contravariant 
and  physical  components  sr»  Sg,  sz,  s^g,  srz,  sQz  of  the  stress  deviation 


tensor.  Also 


. 


r, 


Sf1  *  ss 


gi  ves 
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sr  =  0r  +  P*  se  =  CT0  +  P>  sz  =  °z  +  P. 

sr6  Tr0’  srz  "  Trz’  s0z  ~  T®z’ 
with  p  =  -  |  (ar  +  a0  +  az) . 

The  physical  components  of  the  symmetric  strain  rate  tensor  are 


I  =  lit  d  .  1  3v0  +  V  .  3vz 
r  3r  *  6  '  r  30  r  ’  z  dz~  ’ 


i  i  3vfi  ,  9v  3v 

dr0  =  2  (F  30“  +  sr  ■  7^)  *  =  7  (^T  +  7T1) 


r z  2  v3r  3z 


1  9vfl  l  9v7 

dn  =  -  ( — 2.  +  I  — L) 

0z  2  °z  r  36  ' * 


and  the  physical  components  of  the  skew- symmetric  vorticity  tensor  are 


1  8vr  3  1  9v> 

wr0  =  2r  ^30”  "  3 7  ^rv0^*  wrz  =  z  ^dz 


3v. 

3r" 


-I  A 

w0z  "  2  3z 


1  8vz 

r  30 


) 


A. 4  GOVERNING  EQUATIONS  IN  CYLINDRICAL  POLAR  COORDINATES 

The  equations  of  motion  (2.4.2)  in  terms  of  the  physical  com¬ 
ponents  of  stress  and  velocity  are 


1  3xr0 
r  36 


3t 


rz 


3z 


v°e 


,Dvr 
p  Dt 


r 
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lire  +  1  Hi  +  Hez  .  2  .  ,Dve  .  . 

3r  r  90  9z  r  Tr0  "  p^Dt  "  V  * 


9xrz  1  9tA7  9q7  Dv_ 

rz  x  1  oz  z  ,  rz  _  ,  z  r-  \ 

3r  r  90  97"  r  "  plDt”  "  Fz j » 


where 


DV  -  3Vr  +  v  III  +  *£  !V  +  3vr  V 


Dt  9t  r  9r  r  90  vz  9z  '  r  ’ 


^.!!e  +  V9  f  *e  3  V  ave 

Dt  9t  r  9r  r  90  vz  9z  * 


Dv  9v  9v  vQ  9v  9v 

_ i.  =  _ z  +  v  _z+_6__i.+  v  __z 

Dt  9t  vr  9r  r  90  z  9z  ’ 


p  is  the  density  and  F  ,  FQ,  Fz  are  the  physical  components  of  the  body 
force  per  unit  mass.  For  quasi-static  problems  in  the  absence  of  body 
forces,  the  right  hand  sides  of  equations  (A.4.1a,b,c)  are  zero  and  the 
equations  are  called  the  equations  of  equilibrium. 

The  equation  of  continuity  (2.4.3)  in  terms  of  the  physical 
components  of  velocity  is 


!§■  +  rtr  (rpvj  +  7  fg  (pvfl)  +  (pv7)  =  0. 


9t  r  9r 


0 ' 


9z 


Also  if  the  elastic-perfectly  plastic  solid  is  incompressible,  there  is 
the  equation  of  incompressibility 


9v, 


r  90 


y  ^6 


^bod  9Dn92ds  9rtt  n'r  an^f  joiq  a ;  J6.t2-f26up  10 '  .zzm  Jrnu  idq  so*m  t 


.0  =  (svq)t;  +  (8vqj^*  ,vqr  4^+‘t 
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The  VON  MISES  yield  condition  (2.4.4)  is 


+  s92  +  sz2  ^  2(sre2  ♦  srz2  +  sez2)  -  2k2  (A.4.2) 


and  the  identity  (2.4.5)  is 


sr  +  s0  +  sz  =  °- 


(A. 4. 3) 


The  PRANDTL-RtUSS  constitutive  equations  (2.4.1a)  in  terms 
of  the  components  of  stress  deviation  and  velocity  are 
(a) 

3s  3s  vfi  3s  3s  s  9v  j, 

3t  r  3r  r  '30  ^sr0;  vz  3z  r  L36  3r  vrVJ 


s  3v  3v 


(A. 4. 4a) 


(b) 


3se  !!e  +  !e  ,i!e  +  ,,  w  w  !! 

3t  V  3r  r  ^30  2sr0^  vz  3z 


0 


S  3v  rs  3VQ  -1  3v 

,  _r0.  r r  9_  frw  )-|  _  s  ( — §.  .  J - 1.) 

+  r  ^-30  3r  'rv0'J  s0z  '3z  r  30  ' 


i  9va 

<3T  +  O  '  Xss]’ 


(A. 4. 4b) 


(c) 

3sz+  !fz  +  !e!iz  +  v  !! 
3t“  vr  3r  r  30  z  3z 


3S,  v„  3s  3S  3vf  3v  3Vg  1  3vZ' 

Z  +  v..  rz—  +  —  +  v,  ~  +  srz  ' 3z  3r  ^  s9z  3z  ”  r  30  ^ 


3  V 

■  2^3T  -  XSZL 


(A. 4.4c) 


a*6  ,y 
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(d) 


8sre  !lre  (^fre  * 

9t  vr  9r  +  r  ^  96  sr  ‘  s0^ 


Sr"Sfl  9vr  a  sfl7  9vr  9v7 

+  _ L _ 2_  r _ £1  _  n _ 21  / _ L _ In 

2r  l96  9r  lrVJ  2  l9z  9r  ; 


sr7  9vr  1  9v7  i  l  3v„  vA 

~2~  ^Fz  F  9©“^  =  2y  ^  ^F  FF“  +  FF”  “  F~^  Xsr0-^  (A*4*4d) 


(e) 


9s . 


!fn+  v  !!iz  +  !e  ,^rz  -  .  v  _ 
3t  vr  3r  +  r  39  ~  s9z^  z  3z 


rz 


!ez  r!V  9_  /  n  +  !re  ,!!i  i!!zx 
2r  l90  "  9r  vrVJ  2  ^9z  "  r  90  ; 


s  -ST  9v^  9v7 

+  2  ^Fz”  "  Fr“^  =  2y^T  ^FF~  +  9z“^  "  Xsrz^  * 


•1  ,9vz  ,  9V 


(A.4.4e) 


(f) 


9sft  9sa  va  9sa  9sfl 

9z  +  ,91  +  JL  ( — §1  +  c  )  +  v  — — 

9t  vr  9r  r  90  srz'  vz  9z 


+  JCL  r^Jl  .  9_  (rv  )]  +  !li  (^1  -  ^1) 
+  2r  ^90  9r  ^rv0'J  2  '9z  9 r  ' 


s  ~s  9v^  i  9v  i  9vQ  -I  9v 

♦  -V^  -  1  +  7  *T>  ‘  Xsez].  (A.4.4f) 


2  v 9z  r  90 


1  ,  i  j  . 

where  the  joint  invariant  A  =  o  di  si  1S 

2  k  j 


n  9v  sA  9vfl 

^4  [sr  3T+F<3T+Vr) 

2k 


9v. 

r'  +  Sz  FT 


s  9v  9v0  9v  9v 

+  “IF  +  r  FF~  V0 ^  +  Srz  ^FF  +  Fz“^ 


,  *e*  0v+i9i  v.S0i: 

SV  +  („«  +-~)  —  +  -.' 


1  1$  1 
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9vfl  1  av 

+  S0z  ^az“  +  r  W^‘ 


Details  of  the  derivation  of  the  above  PRANDTL-REUSS  equations 
are  as  follows.  Proceeding  as  in  SECTION  A. 2  above  with  =  r,  x2  =  0 
and  x3  =  z  and 
(1)  1  =  1,  j  =  1, 


Ds 


11 


Dt 


9s11  ,  9s11  ..1  ,  ,9s11  .  J  .21  .  _1  .12\  ..2  ,  9s11  .3 
-Tr+—rv  (“2”  +  r22s  +  r22s  )v  +71“  V 


9x 


9x 


9x 


9Sr  9s  v  9s  9S 

9t”  +  Vr  +  7~  ^90“  "  2sr0^  +  vz  9z~"  * 


/  ml  kl  ,  lm  kK  _  0  1m  11  t 

-(5  9  +  S  g  -  -  2s  g  ulm 


2sre  ^2r  ^39  "  3r  ^rv9^  r  ^30  3r  ^rv0^’ 


2y(d11  -  As11)  =  2u(-37r-  -  As r ) 


9v. 


combine  to  form  equation  (A. 4. 4a); 
(2)  i  =  2,  j  =  2, 


22  22 
Ds  _  9s 


+  (^T“  +  2r?i)  v'  +  ( 


22 


22 


22 


Dt  9t  '  V.l  '  21 1  v  '  '—2 


2  '  -1  1  '9s  +  2T22s12)  v2  +  ^5-  v3 


9x 


9x 


9x 


i  9sa  v0  9se  N  9s0 

r  6  .  w  2.  +  -2.  ( — 2.  +  2s  J  +  v 

=  “2  [^T  Vr  9r  r  l90  ^sr0y 

r 


z  9z 


■]. 


r€ 


.llfi)  4  '  ^  ■  -  -  Et(9«>  4  -  V  4J  e-,8  -  * 
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/  m2  k2  ,  2m  k2v  9  2m  22 

-  (s  g  +  s  g  )u3km  =  -  2s  g  co 


2m 


2  a  sn7  3v  ,  9v 

±_  r _ Hi  ; _ L  _  L_  (ru  )\  +  ( _ 2.  _  J _ 

2  L  2r  l90  9r  ir  e'^  2r  v  90  r  99  ;J> 


9  ^22  ,  22 N  2u  /l  8v6  A  vr  ,  \ 

2u(d  -  xs  )  =  -^  (fae-  +  — -  *s0) 


combine  to  form  equation  (A. 4. 4b); 
(3)  i  =  3,  j  =  3, 


Ds33  3s33  3s33  k 

Dt  -  3t  3xk 


3b,  3s,  v„  3s  3s 

z  +  V  r4  +  —  +  V. 


3t  r  3r  r  36  z  Sz  ’ 


,  m3  k3  .  3m  k3s  _  ,  m3  33 

(s  g  +  s  g  )ukm  -  -  2s  g  u>3m 


9v  9v  9v  ,  9v 

e  ( _ r  _ Zn  .  ( _ 0  _  J _ z\ 

srz  '9z  9r  1  s0z  v9z  r  90  ' * 


2y(d33  -  Xs33)  =  2y(-^  -  Xsz) 


9v. 


combine  to  form  equation  (A. 4. 4c); 
(4)  1  =  1,  j  =  2, 


Ds^=  3s^.+  (iili+  TlsU)^ 
Dt  3t  V  21 


i  ats  ,  u*>  H 1 


Y> 


x  ; 


-  -  m„“<£V£*  *  EVma)  - 
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+  ( 


9s 


12 


+  r 


9x 


22 


9s 


12 


ax' 


1  (!!ri  +  v  !!re  +  ve  8sre  .  v  3sre  .  Vse  , 
r  Vr  3r  VlT+  *z  1T  + T"  V’ 


(sm2gkl  +  s%k2)co.  -  -  (smVVm  +  slmg22u 


1m 


2m 


=  1  [V!i  {!!n .  i_  (rv )}  +  fez  3v 

r  L  2r  *90  9 r  ^rV}  2  ^9z  "  W' 


sr7  9vfl  1  9V7 

+  —A  ( _ sL  _  J _ L\  i 

2  v3z  r  90  )U 


2u(d12  -  AS12)  -  2k  r1  f1  3Vr  +  3V0  \  1 


combine  to  form  equation  (A.4.4d); 
(5)  i  =  1,  j  =  3, 


13  13  13  13  n 

Ds1'3  9s  J  ,  ..1  9s1,3  ,  ,9sIJ  ,  .,1  23n  2  ,  9s1,3  3 

Tt — §r +  v  rr +  (rr +  r22s  >v  +  ~r v 


3X1  3x‘ 


9x' 


!frz  .  u  !!rz  .  ^  ,!frz  ,  ,A„ 

3t  r  3r  r  '  30  s0z^  3z  vz  ’ 


,  m3  kl  x  lm  k3>  23  11  33  11  12  33  11  22 

(s  g  +  s  g  )iokm  =  -  s  g  i»12  -  s  g  u>13  -  s  g  w32  -  s  g  oo31 


S„,  3vr  a  Srfl  3Vfl  1  3v7 

_§i  r _ L  .  £_  (rv  n  +  _E°  f — —  .  J _ 1.) 

2r  l30  3 r  '  0;J  2  l3z  r  30  1 


s  -s  3v  3v 

+  r  Z  ( -  — ) 

2  '  9z  9r 


' 


1  *  X  *1 


T 


i  o  i  q  i  9  V  9  v 

2U(d13  -  Xs13)  =  2y[l  (^+3^)-  Xsrz] 


combine  to  form  (A.4.4e);  and  for 
(6)  1  «  2,  j  -  3, 


Ds23  .  3s23  ,3s23  _2  .m3.  ,1  .  ,3s23  .  r2  .m3.  2 
Dt  3t  V  rmlS  )v  +(^2“  +  rm2s  )v 


3s23  3 

+  v3 

3x 


^Ir +  (^r +  r2i$23)vl  +  (is? +  4s13)v2  +  - t - 

3t  3x  31  3x3  U  3x3 


1  [!!ez  +  v  !!ez  +  !e  (!fez  +  frTfez)+  !fez] 

r  L  9t  r  9r  r  v  90  -  '  v->  ^  J  * 


z  9z 


/  m3  k2  ,  2m  k3\  _  m3  22  2m  33 . 

-  (s  g  +  s  g  )uk  -  -  S  g  <%  -  s  g  <o3m 


13  22  33  22  21  33  22  33 

s  g  co^ -j  “  s  g  o>23  ■  s  g  co^-j  -  s  g 


i  s  9v  *  s  9v  9v7 

1  r  rz  r _ r  9_  /  \  i  +  _r£  / _ L _ l 

7  [f  {96^  "  9r  lrv0j)  2  l9z  9r 


) 


Sfi-ST  9VA  1  9V7 

4.  _§ _ l  t _ 2.  .  I _ L)  1 

2  v9z  r  96 


9Va  i  3v. 


2y(d23  -  XS23)  =  (3^  +  7  90^  “  Xs0z^ 


combine  to  form  equation  (A.4.4f). 


r  . 


+  Sv(ereJl  +~r;)  ♦  rv(r>f^  ♦ 


sv  +(__  _  +  -ir. 


m£“  V  8  -  mS“  V  -  ‘  "  m.^1  3  2  +  6  •' 


t  b  _„2 


APPENDIX  B 


CHARACTERISTIC  STUDY  FOR  PLANE  STRAIN  FLOW  OF  AN 
INCOMPRESSIBLE  ELASTIC-PERFECTLY  PLASTIC  MISES  SOLID 

In  this  appendix,  a  supplementary  study  of  the  characteristics 
of  the  governing  equations  for  steady  state,  quasi-static  plane  plastic 
flow  of  an  incompressible  elastic-perfectly  plastic  MISES  solid  is  given 
using  rectangular  CARTESIAN  coordinates  (x,y,z).  The  procedure  is  dif¬ 
ferent  from  that  used  in  CHAPTER  VI,  SECTION  6.1  and  alternate  forms  for 
the  compatibility  relationships  are  derived. 

B.  1  GOVERNING  EQUATIONS  IN  RECTANGULAR  CARTESIAN  COORDINATES 

For  plane  strain  flow  independent  of  z  and  parallel  to  the  (x,y)- 

plane 


=  0  . 


Since  d  =  0,  s  =  0  and  o  is  equal  to  the  hydrostatic  pressure  -p. 
z  z  z 

Consequently, 


and  the  VON  MISES  yield  criterion  is 


where  k  is  the  yield  stress  in  pure  shear.  sx,  and  are  stress 
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deviation  components. 

The  governing  equations  used  with  appropriate  boundary  con¬ 
ditions  to  determine  sx,  -  p  and  the  velocity  components  vx  and  v 
are  the  PRANDTL-REUSS  equation 


9v. 


(vw  ^r-  +  v..  -  (2k2  -  s  2)  pL  -  -  s  )  s  -x 

x  9x  'y  V  sxy  9y 


9s, 

r  — -  T  V  - - 

x  9x  y  9y 


k2  9Vw 

+  ( —  +  s  )  s  — - 
vy  xJ  bxy  9x 


0  9v 

S  2  — t  =  o  • 

*  )j  u  1 


(B.1.3a) 


the  equations  of  equilibri 


urn 


8Sx  +  !!**.  .  =  o 

3x  9y  9y  * 

(B.l. 3b) 

9SX  9s 

x _ +  °SL  =  o* 

3y  9x  9y  u* 

(B.l. 3c) 

the  equation  of  incompressibility 


9v  9  v 

_*  +  -1  =  o- 

9x  9y  u’ 


(B. 1 .3d) 


and  the  VON  MISES  yield  criterion  ,(B.1.2).  Equation  (B.1.3a)  is  obtained 
from  equation  (2.4.1a)  by  putting  i  =  j  =  1  and  x^  =  x,  x^  =  y,  x^  =  z 
and  noting  that  rectangular  CARTESIAN  tensor  components  are  physical 
components. 

The  yield  criterion  (B.1.2)  is  satisfied  identically  on  defining 


s  =  -  k  sin  2<p 

A 


(B.l ,4a) 


mu  ridr  f  rups  znoM  -iup  j  9f 


ii  r  f  i  jfzat^qmoon  r  nci*J6up9  srlJ 


.(s.r.a)  noMeJMo  ir9\\  2beim  m\  s  1  bn 


=  F  <  # y.  =  sx  ,x  -  ,x  H  ■  1  =  i  -  r  P*,r:f;?uq  ( => ; .  A .  ,  noi;  6up3  mot'* 


t»S  nra  *  -  3  a 


and 
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s  e  k  cos  20 


(B. 1 .4b) 


where  0  +  J  is  the  angle  between  the  direction  of  the  algebraically 
greater  principal  stress  in  the  (x,y)  -  plane  and  the  positive  x-axis. 
Substitution  of  equations  (B.1.4a,b)  into  equations  (B.1.3a,b,c)  together 
with  use  of  equation  (B.1.3d)  yields  respectively 


9v. 


2  f  +  vy  +  2  cos  2*  aT 


,  9v  .  9v 

+  <£  +  2*>  5T-  si"  2*>  °> 


(B. 1 .5a,b,c) 


C05  2*  ||+  sin  2*  f*  +  §7"  0, 


and 

sin  24,  |f  -  cos  2<i>ff +  §f  =  0, 

k 

where  P  =  ^  .  Considering  only  the  case  where  0  c  -  <  1  and  defining 

—  =  sin  2y,  equation  (B.1.5a)  can  be  written  as 

y 

9v 

2  sin  2y  (vx  |f  +  vy  ffi I  +  2  cos  2<t> 

(B.1.6) 

9v  ^vv 

+  (sin  2y  +  sin  24>)  ^  -  (sin  2y  -  sin  2*)  =  0. 


The  four  governing  equations 


used  to  determine  the  unknowns  0,  P,  vx  and 


Yl  r63Ui4«ef8  -liJ  n:S  J  ■  ■>■'  ...  srtj  <  J  +  ♦  9  9  W 

^IsvfJDsqzsi  2bf9fx  (bt.f.8)  nofieupe  ’K)  92U  fIJN 


5  5  $  0  .  ...  di  v  +  ™  #)  ' 


V, 


antnrtsb  bnt  f  >  0  9i9rtw  9ebd  9rti  *fno  eoh«f«ic  .  h  '  9  9  :w 
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Vy  are  equations  (B.1.6),  (B.1.5b,c)  and  (B.1.3d).  For  the  determination 
of  any  existing  characteristics  these  governing  equations  are  supplemented 
with  the  relations 


f  dx  +  gdy  =  dp, 

If  dx  +  ff  dy  =  d*, 

3v  3v 

_  dx  +  _  dy  =  dv 

8v  8v 

3x  dx  +  aT  dy  ‘  dv 


X* 


y 


(B. 1 . 7a ,b,c,d) 


B.2  DETERMINATION  OF  CHARACTERISTICS  AND  COMPATIBILITY  RELATIONSHIPS 

Let  T  be  any  member  of  the  family  of  curves  in  the  (x,y)  -  plane 
governed  by  the  differential  equation 


^  =  tan  c,  C  =  c(x ,y) , 

where  dx  =  cos  c  ds,  dy  =  sin  £  ds  and  ds  is  the  elemental  arclength 
along  r.  Along  r  it  follows  that 

-  -  cos(2*  -  c)  §£  -  s1n(2*  -  c)  ff  (B.2.1) 
with  use  of  equations  (B.1.5b,c)  and  (B.1.7a)  and  that 


g.c«St|t-S1nt|l 


(B.2. 2) 


:  *■  ~ 7 

(Vb 

vb  «  Xu  ^ 


2clIH2MQITAJ3fl  YTIJISITA^MOD  CMA 


30  HOIT-*.i1IM513T3Q  S.8 


(^,x)  9f1J  nr  29VTUO  ^ o  yl  5  !  J  T’o  <9  19  i  yns  ad  .  3sJ 

nor^6up3  .rJna'isWi)  srl;  yd  bgmgvog 
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with  the  use  of  equation  (B.1.7b).  From  equations  (B.2.1)  and  (B.2.2), 


3x  '  sin  2(4>-C)  [sin(2<f>‘?)  3^  +  sin?^  > 


(B.2.3a,b) 


ay 


=  Tin  2(»-  g)  C-  cos(^-?)  al  -  c°s  c  f]  . 


and  hence  from  equations  (B.1.5b,c), 


!!f  =  sin  2U-C)  tsin  4  ^  +  *1"(2*-0  £]  , 
!y  =  sin  2(<j>-c)  [‘  cos  s  3s  '  cos  (2+‘?)  3?] 


(B.2.4a,b) 


It  follows  from  equations  (B.2.3a,b)  and  (B.2.4a,b)  that  |£,  |£- 

9  P  TT 

and  are  not  determinate  when  5  =  <J>  and  c  =  4>  +  j  if 


and 


dP  + .  d(J>  =  0  when  5  =  cf> 


dP  -  d<J)  =  0  when  5  =  <j>  + 


(B.2.5a,b) 


TT 


Hence  the  characteristics  of  the  stress  field  are  defined  by 


and 


dy_ 

dx 


dy. 

dx 


=  tan  <p 


(B. 2.6a  ,b) 


=  -  cot  (p 


Curves  defined  by  equation  (B.2.6a)  are  called  a-lines  and  those  defined 
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by  equation  (B.2.6b)  are  called  3-lines  and  the  compatibility  relation¬ 
ships  are  given  by  equation  (B.2.5a)  and  equation  (B.2.5b)  for  the  a- 
and  3-lines  respectively.  For  the  limiting  case  of  a  rigid-plastic 
solid,  the  characteristics  and  compatibility  relationships,  which  are 
known  as  the  HENCKY  equations,  are  the  same  as  those  above. 

For  any  other  case,  it  follows  from  equation  (B.1.6)  with 
use  of  equations  (B.1.7c,d)  that 


9v 


3x-  [sin(2?-2<j>)  -  sin  2y]  +  sin  2?  sin  2y  [vx  ff  +  vy  fy] 


(B.2.7) 


dv  dv 

+  2  cos  c  [sin  2y  +  sin  2<£]  -  2  sin  c  [sin  2y  -  sin  2(f)]  =  o. 


9v 


If  — -  is  determined, then  by  equations  (B.1.7c,d)  and  (B. 1.3d),  so  are 


9  x 


3v .  9v  .  9v. 


x  ,  and  -r-X-  .  In  equation  (B.2.7)  it  is 

\  8vx  9vx  9vy 

seen  that  when  c  =  4>  +  Yorc  =  <l>-Y  +  |-»  (and  hence  -  ,  3x 


the  partial  derivatives  ,  9x 


3v 


and  — ^)  is  not  determinate  if 
oy 


dv. 


;in  2c  sin  2y  [vx  fy]  +  cos  £  [sin  2<Y  +  sin  2(^  ds 


dv 


(B.2.8) 


-  sin  c  [sin  2y  -  sin  2<j>]  -  0. 


Consequently,  the  characteristi cs  of  the  velocity  field  are  defined  by 


^  =  tan  (4>  +  y) 


p-  =  -  cot  (<f>  -  y). 
dx 


and 


(B.2.9a  ,b) 


« 


wV6 
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Curves  defined  by  equation  (B.2.9a)  are  called  a-j-lines  and  those  de¬ 
fined  by  equation  (B.2.9b)  are  called  3-j-lines.  Also,  since  y  i  0, 
substitution  of  equations  (B.2.3a,b)  into  equation  (B.2.8)  is  permissible 
and  the  result  is 

sin2U-”j'Y  [{vx  sTn(24,-C)  -  vy  cos(2*-?)}  |f-  +  {vx  sin?-vy  cos?} 

dv  dv 

+  cosc[sin2y+sin(t>]  ^ - sinc[sin2y-sin2<j)]  =  0.  (B.2.10) 

For  an  a-j-line,  c  =  ^  +  Y  and  equation  (B.2.10)  reduces  after  simplification 
to 


[vx  cos  (4>-y )  +  vy  sin  (<J>-y)] 

HP 

+  [v  cos  (<j>+y)  -  vx  sin  ( ^  =  0.  (B.2.11) 

On  defining 


and 


U  =  v¥  cos  (4>-y)  +  v  sin  (<j>-y) 

X  y 

V  =  v  cos  (4>+y)  -  vY  sin  (<j)+y), 
y  * 


equation  (B.2.11)  becomes 


dU 

ds 


+ 


(B.2.12a,b) 


or  equivalently  » 


-  USnra-rSnrs^nra  -  [tnr8+YSri^2]3«03  + 
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dU  +  W  dp  =  °>  (B.2.13) 

the  compatibility  relation  along  an  c^-line.  Similarly  for  a  6,-line 
for  which  .?  =  4>  -  y  +  j  ,  equation  (B.2.10)  reduces  to 

3T  [vy  cos(<t>+y)  -  vx  sin  (<f>+y)] 

+  [vx  cos +  v  si n (<J>-y) ]  =  0. 

With  use  of  equations  (B.2.12a,b),  this  equation  is  written  simply  as 


or  equivalently 


dV  +  Ik  dp  '  0  *  (B.2.14) 

the  compatibility  relation  along  a  3-j-line. 

Equations  (B.2.13)  and  (B.2.14)  are  alternate  forms  of  the 
compatibility  relations  along  an  o^-line  and  a  3^ -line  derived  in  CHAPTER 
VI,  SECTION  6.1.  From  equations  (B.2.12a,b),  it  is  observed  that  V  is 

I 

the  component  of  the  velocity  perpendicular  to  the  a-|-line  and  U  is  the 
component  of  velocity  perpendicular  to  the  3-|-line.  If,  however,  the 
velocity  is  resolved  into  two  components  u  and  v  where  u  is  the  component 
along  the  a-j-line  and  v  the  component  along  the  Byline,  then 


r!  flu  1.10,1.  ftorJ  Js-t  i- r rd  )sqnoa  art 

.0  =  •  [(r-*)<iU  v  +  (y-*)2<»  ,v]  + 
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and 


U  =  u  cos  2y 


V  =  v  cos  2y . 


Use  of  these  relations  in  equations  (B.2.13)  and  (B.2.14)  yield  re¬ 
spectively 


and 


du  +  v  ^-  =  0 


along  an  a-j-line 


dv  +  u  ^-  =  0 


along  a  3-j-line  , 


thus  recovering  the  form  of  the  compatibility  relations  along  the  velocity 
characteristics  derived  originally  in  CHAPTER  VI. 
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